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ABSTRACT 

The  objective  of  this  report  is  to  summarize  water  quality 
constituents  in  the  AOSERP  study  area  and  to  examine  relationships 
between  these  constituents  and  changes  in  land  formation,  hydrology,  and 
development. 

Summaries  of  routine  parameters,  nutrients,  and  metals  are 
presented  for  sampling  sites  along  the  Athabasca  River  to  detect 
longitudinal  and  seasonal  changes  in  water  quality  and  to  determine  the 
effects  of  point  source  effluents  on  the  river.  Changes  in  ion 
concentrations  on  the  Athabasca  River  were  attributed  to  inputs  from  the 
Clearwater  River  rather  than  industrial  saline  discharge  via  Poplar 
Creek.  The  downstream  effect  of  the  Fort  McMurray  sewage  treatment 
plant  was  limited  to  the  mainstem  upstream  of  the  confluence  with  the 
Muskeg  River.  Elevated  concentrations  of  nickel  and  vanadium  occurred 
during  high  flow  periods,  indicating  that  concentrations  of  these  metals 
in  surface  waters  were  associated  with  weathering  of  natural  bitumen 
rather  than  with  industrial  effluents.  A principal  component  analysis 
(PCA)  was  used  for  the  simultaneous  examination  of  selected  water 
quality  parameters  on  the  Athabasca  River.  Sites  exhibiting  similar 
water  quality  characteristics  were  delineated  on  schematic  maps  of  the 
river. 

Baseline  data  and  relationships  among  parameters  also  are 
presented  for  east,  west,  and  south  drainages  entering  the  Athabasca 
River  between  Fort  McMurray  and  Embarras  Airport,  as  well  as  the 
Athabasca  Delta  drainage.  An  overall  analysis  of  the  four  regions  was 
conducted  using  PCA  to  delineate  those  sites  with  similar  water  quality 
characteristics.  Site  groupings  often  reflected  the  geological  type  of 
the  region. 
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1. 


INTRODUCTION 


1.1  ATHABASCA  OIL  SANDS  AREA 

The  Alberta  Oil  Sands  Environmental  Research  Program  (AOSERP) 
study  area  encompasses  about  28  400  km2  in  northeast  Alberta  (Figure  1). 
The  area  occurs  between  56°  and  59°  north  latitude  and  111°  and  113° 
west  longitude,  excluding  Wood  Buffalo  National  Park. 

The  AOSERP  study  area  includes  most  of  the  mineable  area  of 

the  Athabasca  Wabiskaw-McMurray  deposit  of  oil  sands.  Suncor  Inc.  and 

Syncrude  Canada  Ltd.  are  the  two  oil  sands  mining  and  processing  plants 
in  the  study  area.  Fort  McMurray,  located  in  the  southern  part  of  the 
study  area,  is  a major  population  centre.  Anzac,  Fort  Chipewyan,  and 
Fort  MacKay  are  small  settlements.  Increased  demands  on  surface  waters 
associated  with  mining  development  include  recreational , domestic,  and 
commercial  uses.  Despite  these  demands  there  is  a concern  to  minimize 

impacts  of  development  on  water  quality  within  the  Athabasca  Basin. 

This  report  summarizes  regional  water  quality  data  and,  when 
applicable,  impacts  of  development  on  surface  water  quality. 

1.1.1  Regional  Surface  Water  Quality  Monitoring  Program 

Several  specific  water  quality  studies  related  to  industrial 
development  were  initiated  in  1976  throughout  the  AOSERP  area.  In 
1977,  a three-year  co-ordinated  regional  sampling  program  was  conducted 
by  the  Water  Quality  Control  Branch,  Pollution  Control  Division,  Alberta 
Environment.  Sampling  sites  were  selected  in  four  major  regions  of  the 
study  area  including: 

1.  Mainstem  Athabasca  River 

2.  Tributaries 

3.  Lakes 

4.  Athabasca  Delta 

All  recorded  data  were  stored  in  the  National  Water  Quality 
Data  Bank  (NAQUADAT).  Water  quality  monitoring  activities  have 
continued  from  1981  to  the  present,  focussing  on  the  Athabasca  mainstem 
and  the  Athabasca  Delta. 
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1.1.2  Reports 

Several  water  quality  AOSERP  reports  have  been  published  which 
dwell  on  specific  water  quality  aspects  within  the  study  area  (Akena 
1979;  Costerton  and  Gessy  1979;  Nix  et  al . 1979).  Seidner  (1980) 
summarized  regional  water  quality  in  the  AOSERP  study  area  based  on  1976 
and  1977  data  sets.  Description  of  sampling  sites  and  analytical 
procedures  used  in  the  regional  sampling  program  have  been  prepared  by 
Akena  (1980).  Non-AOSERP  surface  water  quality  data  for  the  region  have 
been  compiled  by  Akena  and  Christian  (1981).  Recently,  a compilation  of 
non-NAQUADAT  and  NAQUADAT  water  quality  data  (historical  to  1982)  has 
been  prepared  for  the  Athabasca  River  Basin  (Alberta  Environment, 
Planning  Division  1983a,  1983b).  In  these  reports,  the  AOSERP  study 
area  data  are  presented  for  Reach  7 of  the  Athabasca  Basin.  A detailed 
water  quality  summary  for  the  mainstem  Athabasca  within  the  AOSERP  study 
area  is  being  prepared  (Akena  in  prep.).  The  present  report  will 
emphasize  regional  characteristics  of  the  study  area. 

1.1.3  Report  Objectives 

The  objective  of  the  co-ordinated  regional  water  quality 
program  was  to  establish  baseline  data  for  aquatic  resources  in  the 
AOSERP  study  area,  thus  providing  a standard  against  which  any  impact 
could  be  measured.  Oil  sands  development,  however,  had  already 
occurred:  Suncor  Ltd.  (originally,  Great  Canadian  Oil  Sands  Ltd.)  began 

operations  in  1967;  Syncrude  Canada  Ltd.  began  operations  in  1977. 
Because  additional  oil  sands  development  was  planned  for  several 
watersheds  in  the  area,  a broad  regional  data  base  was  required.  This 
regional  report  is  a summary  of  water  quality  constituents 
characteristic  of  watersheds  associated  with  the  Athabasca  River  in  the 
north,  east,  west,  and  south. 

Specific  objectives  of  the  report  include: 

1.  descriptions  of  regional  baseline  states; 

2.  identification  of  physical,  chemical,  and  biological 
processes  affecting  surface  water  quality; 
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3.  identification  of  significant  correlations  between  water 
quality  constituents  and  hydrology; 

4.  evaluation  of  the  significance  of  natural  levels  of 
physical  and  chemical  parameters  to  biological 
communities;  and 

5.  when  possible,  an  assessment  of  impacts  of  development  on 
surface  water  quality. 

1.2  DESCRIPTION  OF  THE  ATHABASCA  OIL  SANDS  AREA 

1.2.1  Ecoregions 

Ecoregions  are  land  areas  "characterized  by  a distinctive 
regional  climate  as  expressed  by  vegetation"  (Subcommittee  on 
Biophysical  Land  Classification  1969).  Strong  and  Leggat  (1981) 
delineated  ecoregions  for  Alberta  using  map  units  identified  by 
vegetation,  climate  and,  to  a lesser  extent,  soils.  The  approach  used 
was  similar  to  that  of  the  Canadian  Committee  on  Ecological  Land 
Classification  (Hills  1976). 

Turchenek  and  Lindsay  (1982)  used  the  ecoregions  defined  by 
Strong  and  Leggat  (1981)  as  a framework  for  studying  soils  in  the  AOSERP 
area.  Vegetation  (a  reflection  of  climate)  is  often  used  to  describe 
ecoregions  when  climatological  data  are  insufficient.  Turchenek  and 
Lindsay  (1982)  subdivided  ecoregions  into  ecodistricts  and  physiographic 
units  (adapted  from  Pettapiece  1981).  Ecodistricts  are  land  areas  with 
distinct  patterns  of  "relief,  geology,  geomorphology  and  associated 
regional  vegetation"  (Turchenek  and  Lindsay  1982).  The  Boreal  Mixedwood 
and  the  Northern  Boreal  are  the  two  ecoregions  represented  in  the  AOSERP 
study  area. 

Boreal  Mixedwood,  which  constitutes  the  largest  land  area 
(43%)  of  Alberta,  is  characterized  by  deciduous  forests  underlain  by 
gray  luvisolic  soils.  Jack  pine  communities  are  found  on  sandy 
material;  aspen  poplar/white  spruce  occur  on  moderately  drained  soils. 
More  than  70%  of  the  annual  precipitation  in  this  region  occurs  in 
summer.  July  is  the  wettest  month. 
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The  Northern  Boreal  Ecoregion  is  characterized  by  deciduous 
(aspen)  and  coniferous  (white  spruce)  forests  underlain  by  gray 
luvisolic  soils.  The  Northern  Boreal  region  is  distinguished  from  the 
Mixedwood  by  an  increase  in  abundance  of  white  spruce  and  a colder, 
drier  climate. 

1.2.2  Cl imate 

Longley  and  Janz  (1978)  present  summary  climatological  data 
for  the  AOSERP  study  area.  A climatic  summary  for  the  Northern  Boreal 
Forest  and  the  Boreal  Mixedwood  Forest  is  presented  in  Table  1. 

Climatologists  rely  on  a continuous  record  of  30  years  to 
develop  sets  of  normal  values  for  weather  stations.  Such  long-term 
climatic  information  is  available  for  Fort  McMurray  (1944  to  present). 
Other  data  are  available  for  Embarras  Airport  (February  1943  to 
September  1962)  and  Fort  Chipewyan  (October  1962  to  present),  plus 
several  forestry  lookout  stations  (summer  data  only).  The  climate  of 
the  AOSERP  study  area  is  continental,  characterized  by  cold  winters  with 
relatively  little  snow,  and  cool,  wet  summers.  Most  precipitation 
occurs  in  summer  as  major  rainfalls. 

Mean  daily  temperatures  for  both  Fort  McMurray  and  Fort 
Chipewyan  are  16.3°C  in  July  and  -21.5°  and  -26.2°C  for  the  two  locales 
in  January.  Fort  McMurray  and  Fort  Chipewyan  receive  223  and  226  days 
of  frost  each  year,  respecti vely.  At  Fort  McMurray,  the  average  annual 
precipitation  is  435  mean  (305  mm  rainfall  and  140  cm  snowfall).  The 
amount  of  precipitation  decreases  northward  from  Fort  McMurray  to  Fort 
Chipewyan. 

1.2.3  Vegetation 

Thompson  et  al . (1978)  present  a classification  scheme  for  the 
AOSERP  study  area,  based  on  major  vegetation  types  described  by  Stringer 
(1976).  Turchenek  and  Lindsay  (1982)  identified  the  vegetation  and  soil 
types  associated  with  each  drainage  system  and  presented  descriptions  of 
ecological  units. 


Table  1.  Summary  of  climatic  data  for  Alberta  ecoregions  (adapted  from  Strong  and  Leggat  1981). 
Values  expressed  as  (minimum  to  maximum)  mean. 
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Northern  Alberta  is  located  within  the  Boreal  Forest  region 
(Rowe  1972).  The  four  subdivisions  of  the  Boreal  Forest  which  occur 
within  the  AOSERP  study  area  are  Boreal  Mixedwood,  Athabasca  South, 
Northwest  Transition  Zone,  and  Upper  Mackenzie. 

The  Boreal  Mixedwood  Forest  region  occupies  most  of  the  study 
area.  Here,  aspen  poplar  (Populus  tremuloides  Michx)  and  white  spruce 
(Picea  glauca  Voss)  are  the  main  tree  species.  Aspen  poplar,  a species 
capable  of  rapid  regeneration  in  disturbed  regions,  covers  the  largest 
areal  extent  of  the  study  area  (Thompson  et  al.  1978).  White  spruce  is 
the  climax  species  in  well  drained  areas.  Jack  pine  (Pinus  banksiana 
Lamb.)  occupies  the  eastern  sand  hills.  Larch  (Larix  laricina  Koch)  and 
black  spruce  (Picea  mariana  BSP)  grow  in  open  muskeg. 

The  Athabasca  South  Forest  region  occurs  in  the  northeast  part 
of  the  study  area  (Richardson  and  Firebag  uplands).  Jack  pine  inhabits 
dry,  sandy  areas  and  black  spruce,  wetlands. 

The  Northwest  Transitional  Forest  region  is  found  north  of 
Lake  Athabasca  within  the  Precambrian  Shield.  Black  spruce  is  the 
dominant  tree  species. 

The  Upper  Mackenzie  Forest  region  occurs  in  the  north  of  the 
study  area.  Alder  (Alnus  tenuifolia  Nutt.),  balsam  poplar  (Populus 
balsamifera  L.),  and  white  spruce  are  found  in  lowlands  bordering 
rivers.  Jack  pine  and  aspen  poplar  with  black  spruce  and  larch  occur  in 
the  uplands. 

1.2.4  Physiography 

The  two  major  physiographic  provinces  within  the  AOSERP  study 
area  are  the  Interior  Plains  and  the  Precambrian  Shield  (Atlas  of 
Alberta  1969).  Pettapiece  (1981)  has  subdivided  the  provinces  into 
regions,  sections,  and  districts.  Detailed  descriptions  of  these  units 
are  presented  by  Thompson  et  al.  (1978)  and  Turchenek  and  Lindsay 
(1982). 

The  physical  geography  of  a region  is  based  on  surficial 
deposits,  geology,  and  vegetation.  Chemical  characteristics  of  surface 
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water  differ  significantly  between  physiographic  regions  (Dillon  and 
Kirchner  1975;  Kirchner  1975)  and  geological  type  (Naumann  1929). 
Accordingly,  water  quality  constituents  of  surface  waters  within  the 
same  physiographic  region  should  be  similar. 

1.2.5  Geology 

1.2.5. 1 Bedrock  geology.  A bedrock  geology  map  of  the  AOSERP  area, 
presented  in  Figure  2,  was  derived  by  Turchenek  and  Lindsay  (1982). 

The  geology  of  the  AOSERP  study  area  is  characterized  by 
material  from  the  Pleistocene  and  Recent  eras  (Carrigy  1959;  McPherson 
and  Kathol  1977).  Since  surface  mining  for  oil  sands  does  not  extend 
beneath  the  McMurray  formation,  most  regional  studies  have  focussed  on 
rocks  of  the  Cretaceous  age  or  younger  (Mellon  and  Wall  1956).  Saline 
water  may  occur  in  the  basal  McMurray  formation,  presenting  ecological 
and  engineering  concerns  in  oil  sands  development. 

Cretaceous  rocks  (sandstone  and  shale)  dominate  most  of  the 
study  area.  Sandstones  including  the  oil  sands  are  exposed  in  the 
Athabasca  River  valley.  Shales,  with  some  sandstone,  occupy  the 
uplands.  Devonian  rocks  (limestone,  dolomite  and  gypsum)  occur  in  the 
north  and  Precambrian  granite,  the  extreme  northeast. 

1.2. 5. 2 Surficial  geology.  Thompson  et  al.  (1978)  developed  surficial 
geological  maps  (1:50  000)  for  the  AOSERP  study  area  that  provided 
detailed  information  on  deposits  and  land  forms.  Surficial  geological 
maps  also  have  been  prepared  for  specific  mining  areas  (McPherson  and 
Kathol  1977).  Bayrock  (1971,  1972a,  1972b)  and  Bayrock  and  Riemchen 
(1974)  mapped  (1:250  000)  surficial  geology  on  several  National 
Topographic  Sheets  (74C,  74L,  74M,  74P,  841).  Turchenek  and  Lindsay 
(1982)  prepared  a soil  parent  materials  map  (1:500  000)  and  detailed 
soil  maps  (1:126  720)  for  the  AOSERP  area. 
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Figure  2 . Bedrock  geology  of  the  AOSERP  study 

QreQ  . ( Adopted  from  Turchenek  and  Lindeoy  , 1982  ) 
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Surficial  geology  is  determined  by  glacial  and  post  glacial 
deposition  and  erosion.  Soil  classification  is  based  on  genetic  types 
of  material  on  which  soil  develops  (Turchenek  and  Lindsay  1982).  Three 
main  types  of  surficial  geology  in  the  AOSERP  study  area  include 
(Thompson  et  al . 1 978) : 

1.  Ground  and  hummocky  morain  (till)  - Birch  Mountains  and 
Methy  Portage  Plain; 

2.  Silt  and  clay  lake  deposits  - Birch  Mountains;  and 

3.  Silt  and  gravel  deposits  from  fluvial,  lacustrine,  and 
aeolian  sources  - lower  valley  areas,  eastern  study  area, 
and  delta. 

Surficial  geology  affects  the  rate  of  water  percolation, 
groundwater  storage,  surface  runoff,  land  erosion,  and  vegetation  cover. 
All  these  factors  affect  water  quality. 

1.2.6  Drainage  Systems 

A map  of  watersheds  associated  with  the  AOSERP  study  area 
presented  by  Neill  and  Evans  (1979)  is  reproduced  in  Figure  3.  Neill 
and  Evans  (1979)  consider  the  following  drainages  to  be  major  components 
of  the  surface  water  drainage  within  the  area: 

1.  The  Athabasca  River  enters  the  study  area  from  the 
southwest  and  extends  to  the  Athabasca  Delta. 

2.  The  Clearwater  River,  which  originates  in  Saskatchewan, 
enters  the  study  area  from  the  southeast  and  joins  the 
Athabasca  River  at  Fort  McMurray. 

3.  Major  tributary  streams  south  of  Fort  McMurray  include 
the  Horse,  Hangingstone,  and  Christina  rivers.  The  Horse 
River  enters  the  Athabasca  River.  The  Hangingstone  and 
Christina  rivers  are  tributaries  of  the  Clearwater. 

4.  Major  tributaries  flowing  from  the  east  to  the  Athabasca 
River  north  of  Fort  McMurray  include  the  Steepbank, 
Muskeg,  and  Fi rebag  rivers.  For  the  purpose  of  this 
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Figure  3.  Surface  water  drainage  of 
the  AOSERP  study  area  . 

( Adopted  from  Ntlll  and  Evans  , 1979  } 
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report,  the  Richardson  River  will  be  included  in  the 
eastern  region.  The  Richardson  River  originates  in 
Saskatchewan  and  enters  the  Athabasca  mainstem  in  the 
delta  area. 

5.  Major  tributaries  flowing  from  the  west  to  the  Athabasca 
River  north  of  Fort  McMurray  include  the  Beaver/Poplar, 
MacKay,  and  Ells.  Several  smaller  tributaries  drain  the 
east  slopes  of  the  Birch  Mountains. 

6.  Water  bodies  associated  with  the  Athabasca  Delta, 
including  lake  distributary  channels,  will  be  discussed 
in  the  section  of  this  report  concerning  the  northern 
region. 

1.2.7  Surface  Water  Hydrology 

Stream  gauging  data  for  tributaries  in  the  AOSERP  study  area 
were  complied  by  Loeppky  and  Spitzer  (1977),  Warner  and  Spitzer  (1979), 
and  Warner  (1979).  Neill  and  Evans  (1979)  provided  an  overview  of 
surface  water  hydrology  of  the  study  area.  Sekerak  and  Walder  (1980) 
included  some  hydrological  data  in  their  biophysical  inventory  of  major 
tributaries  in  the  AOSERP  area.  Schwartz  (1979)  summarized  the 
hydrology  of  the  Muskeg  River.  Locations  of  stream  gauging  stations 
within  the  study  area  are  indicated  in  Figure  4. 

Gauging  stations  on  the  Athabasca  River  above  Fort  McMurray 
and  the  Clearwater  above  Fort  McMurray  have  been  in  operation  si  ne'e 
1957.  The  downstream  gauge  at  Embarras  was  installed  in  1971. 
Short-term  records  (since  1975)  are  available  for  most  of  the 
tributaries  entering  the  mainstem. 

1.2.7. 1 Water  balance:  precipitation,  runoff,  and  evapotranspiration. 

Water  is  supplied  by  precipitation  and  lost  through  runoff  and 
evapotranspi ration.  Water  also  may  remain  as  surface  or  subsurface 
storage,  but  over  the  long  term  this  component  is  assumed  to  be 
negligible  (Neill  and  Evans  1979).  Thus,  a water  balance  equation 
may  be  expressed  as  follows: 
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P = R + ET 

where:  P = total  precipitation; 

R = runoff;  and 
ET  = evapotranspiration. 

When  values  for  all  three  components  are  known,  there  are  often 
Inconsistencies  in  the  water  balance  because  of  inadequate 
estimates  of  precipitation  and  evaporation  (Neill  and  Evans  1979; 
Hydrological  Atlas  of  Canada  1978).  Errors  in  runoff  estimates  (stream 
discharge  data)  are  minor. 

Neill  and  Evans  (1979)  presented  the  mean  annual  water  balance 
for  Fort  McMurray  and  the  Athabasca  Delta  based  on  data  provided  in 
the  Hydrological  Atlas  of  Canada  (1978).  In  the  AOSERP  study  area, 
evaporation  accounted  for  78  to  80%  of  precipitation;  accordingly, 
runoff  accounts  for  20  to  25%  of  precipitation.  This  high  evapotranspi- 
ration: runoff  ratio  results  from  high  summer  temperatures , long  hours  of 
sunshine,  low  humidity,  abundant  vegetation,  open  water,  and  muskeg 
(Neill  and  Evans  1979). 

The  annual  water  balance  determined  for  the  region  using  the 
limited  local  data  set  (1976-77)  generally  supported  the  findings 
presented  in  the  Hydrological  Atlas  of  Canada  (1978).  Runoff  and 
derived  evapotranspiration  (precipitation-runoff)  estimates  from 
the  limited  data  set  were  lower  than  the  long-term  means  (Neill  and 
Evans  1979). 

1.2. 7. 2 Interaction  of  surface  flows  with  groundwater.  The  Alberta 
Research  Council  conducted  a regional  hydrogeological  evaluation  of  the 
Athabasca  Oil  Sands  area  (Hackbarth  and  Nastasa  1979).  A total  of  75 
observation  wells  were  installed  at  15  locations  within  the  study  area 
during  1974-75  and  1975-76.  Hackbarth  and  Nastasa  (1979)  estimate  that 
groundwater  contributions  entering  the  Athabasca  River  between  Fort 
McMurray  and  Embarras  are  negligible. 
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1.2. 7. 3.  Streamflow  balance.  The  Athabasca  (17.2  km3/a)  and  Clearwater 
(4.3  km3/a)  rivers  provide  the  major  discharge  inputs  into  the  AOSERP 
study  area  (Neill  and  Evans  1979).  The  contribution  of  flows  from 
tributary  streams  are  relatively  minor  (Figure  5).  Flow  from  the 
southern  tributaries  (Horse  and  Hangingstone  rivers)  contribute  about 
0.4  km3/a.  Flow  contribution  from  the  west  bank  (0.8  km3/a)  is  less 
than  the  east  bank  (1.2  km3/a)  tributaries.  The  mean  annual  inflow  from 
the  Athabasca  River  system  to  Lake  Athabasca  is  23.9  km3/a. 

Several  northern  rivers  drain  directly  into  Lake  Athabasca, 
the  largest  of  which  is  the  Fond  du  Lac  River  (9.4  km3/a).  The  average 
outflow  from  Lake  Athabasca  and  the  delta  to  the  Slave  River  is  about 
45.3  km3/a. 

1.2. 7. 4 Response  of  streamflow  to  snowmelt  and  rainfall.  Neill  and 
Evans  (1979)  developed  runoff  snowmelt  and  rainfall  coefficients  for 
various  basins  within  the  AOSERP  study  area. 

Hydrographs  for  streams  in  the  area  depict  a snowmelt  response 
beginning  at  the  end  of  March  and  peaking  in  mid-April.  Estimated 
snowmelt  runoff  coefficients  (volume  of  runoff  divided  by  volume  of 
snowwater  equivalent  on  the  ground)  ranged  from  18%  (Muskeg  River)  to 
56%  (Joslyn  and  MacKay  rivers).  These  values  are  underestimates 
because  some  snowmelt  is  stored  in  the  soil  and  subsequently  used  by 
vegetation  or  is  lost  to  evaporation  (Neill  and  Evans  1979). 

Estimated  rainfall  runoff  coefficients  were  calculated  using 
the  following  formula  (Neill  and  Evans  1979): 

RR  = SF  - SR 

where:  RR  = rainfall  runoff  coefficient; 

SF  = stream  flow  (01  April  to  31  December  1976);  and 

SR  = snowmelt  runoff. 

Values  for  rainfall  runoff  coefficients  were  low  (4  to  19%) 
for  most  basins  except  for  those  in  the  southern  region  (Horse  River, 
28%;  Hangingstone,  32%).  Rainfall  runoff  coefficients  were  lower  than 
snowfall  runoff  coefficients  for  basins  on  the  west  bank  of  the 
Athabasca  River. 
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HORSE  a HANGtNGSTONE  RIVER 


Athabosco  River  inflow  to  lake  system  » 17.2+0.4+4.3+0.8  + 1.2=23.9 

Lake  outflow  to  Slave  River  : 23.9  +1.2  +2.6+9.4+8.2  = 45.3 

Units  of  km3  ( = I09  m3  ) 


Figure  5.  Mean  annuGl  streamflow  balance  . 

( Adapted  from  Neill  and  Evans  , 1979  ) 
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1 .2.7.5  Erosion  and  sedimentation.  Suspended  materials  follow  a 
seasonally  cyclical  pattern,  with  high  levels  occurring  during  high  flow 
periods.  Suspended  sediment  data  have  been  recorded  at  WSC  gauging 
stations  along  the  Athabasca  River  within  the  AOSERP  study  area  at  sites 
below  Fort  McMurray  and  at  Embarras.  Other  sediment  samples  have  been 
obtained  at  infrequent  intervals  throughout  the  study  area. 

Doyle  (1977)  described  the  suspended  sediment  regime  of  the 
Athabasca  River  using  suspended  sediment  rating  curves.  During  high 
flow  conditions,  the  upstream  section  of  the  Athabasca  River  within  the 
study  area  has  a higher  capacity  to  transport  sediments;  whereas  during 
lower  flows,  the  carrying  capacity  of  the  downstream  section  becomes  more 
important. 

The  average  annual  sediment  load  upstream  of  Fort  McMurray 
is  13.6  x 106  tonnes/a  (Thurber  et  al . 1975),  of  which  5%  is  bedload. 

The  average  annual  sediment  load  from  the  Clearwater  River  is 
0.9  x 106  tonnes/a.  Griffiths  and  Walton  (1978)  estimated  that  the 
average  annual  sediment  load  entering  the  Athabasca  River  between  Fort 
McMurray  and  Embarras  from  tributary  streams  represents  5%  of  the  load 
at  the  downstream  site. 

Most  tributaries  entering  the  Athabasca  maintream  have 
characteristically  stable  channels.  Little  erosion  of  the  bank  materials 
occurs  within  the  basin  because  of  dense  vegetation  cover  and  low  runoff. 
Also,  downstream  reaches  are  entrenched  in  bedrock.  Erosion,  however, 
is  prevalent  in  streams  draining  the  east  slopes  of  the  Birch  Mountains. 
Here,  there  is  extensive  slumping  in  the  headwater  reaches,  resulting 
in  eroded  material  being  transported  into  alluvial  fans  at  the  confluence 
with  the  Athabasca  mainstem  (Griffiths  and  Walton  1978;  Corkum  et  al . 1982). 

1.2. 7. 6 Impacts  of  oil  sands  development  on  hydrology.  Hackbarth  and 
Nastasa  (1979)  outline  hydrogeological  difficulties  associated  with  oil 
sands  development,  including  groundwater  depressurization , natural  gas 
occurrence,  saline  water,  and  induced  infiltration. 
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Oil  sands  development  places  increasing  demands  on  water 
supply.  Surface  waters  and  groundwaters  are  used  in  processing  and  are 
later  discharged  into  waterways.  Groundwater  associated  with  oil  sands 
must  be  disposed  of  before  the  sands  can  be  mined.  This  water  is  highly 
saline  and  frequently  toxic  to  aquatic  organisms  (Tsui  et  al.  1980). 

Depressurization  (reduction  of  pore  pressure)  and  dewatering 
(removal  of  pore  water)  are  required  to  ensure  pit  slope  stability  and 
dry  pit  conditions.  Chemical  characteristics  of  mine  depressurization 
water  (high  saline  levels,  range  of  heavy  metals,  and  organic 
constituents)  have  been  described  by  McMahon  et  al . (1977),  Lake  and 
Rogers  (1979),  and  Giles  et  al.  (1979).  Natural  gas  has  been  produced 
from  the  basal  McMurray  aquifer  and  may  present  problems  for 
depressurization  (Hackbarth  and  Nastasa  1979).  High  saline  levels  will 
be  most  severe  during  low  flow  (winter)  conditions  and  where  groundwater 
constitutes  base  flows. 

Infiltration  may  occur  from  the  Athabasca  River  to  oil  sands 
mines  which  are  located  close  to  the  mainstem  (Hackbarth  and  Nastasa 
1979).  Infiltation  may  reduce  saline  concentrations  of  groundwater 
during  mining  operations,  but  may  be  detrimental  once  operations  cease 
and  water  moves  from  the  basal  McMurray  aquifer  to  the  Athabasca  River. 
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2.  APPROACH 

This  report  is  a study  of  water  quality  parameters  within  and 
among  four  sub-basins  of  the  AOSERP  study  area  and  the  Athabasca 
mainstem.  The  sub-basins  include  south,  west,  and  east  drainages 
entering  the  Athabasca  River  between  Fort  McMurray  and  Embarras  Airport, 
and  the  Athabasca  Delta  drainage. 

The  objective  of  this  report  is  to  establish  relationships 
among  water  quality  parameters  which  can  be  used  to  reflect  changes  in 
water  quality  owing  to  developmental  impact  or  hydrological  regime. 
Within  each  region  or  sub-basin  such  relationships  will  be  established 
between  water  quality  and  surface  water  discharge.  Other  relationships 
will  be  established  (depending  on  available  data)  among  habitat 
characteristics  of  rivers,  land  formations  (ecodistrict,  geology) 
through  which  rivers  flow,  and  benthic  invertebrates  (environmental 
indicators)  inhabiting  the  rivers. 

This  report  is  organized  by  region,  concluding  with  an  overall 
summary  analysis.  The  aim  is  to  establish  relationships  among 
parameters  in  regions  without  perturbations  and  then  to  determine  the 
effect  of  point  source  effluents  on  the  Athabasca  mainstem. 

When  analysing  data  by  regions,  there  is  a difficulty  in  that 
the  study  area  may  be  so  homogenous  that  no  significant  relationships 
exist.  Also,  the  various  data  sets  have  not  always  been  retrieved  from 
the  same  location,  with  the  same  frequency,  or  within  the  same  time 
frame. 

2.1  METHODS 

The  following  information  will  be  presented  for  each  region: 

1 . base! ine  data; 

2.  analysis;  and 

3.  impact  of  development. 
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2.1.1  Baseline  Data 


2. 1.1.1  Land  formations.  Maps  of  watersheds,  ecoregions,  ecodistrict, 
physiology,  and  bedrock  geology  that  were  prepared  for  the  AOSERP  study 
area  by  other  workers  are  presented  and  enlarged  for  the  regions.  The 
delineation  of  watersheds  were  obtained  from  Neill  and  Evans  (1979). 
Ecoregion  and  ecodistrict  maps  as  well  as  physiographic  regions  were 
obtained  from  Turchenek  and  Lindsay  (1982)  based  on  work  by  Pettapiece 
(1981).  The  bedrock  geology  maps  were  derived  by  Turchenek  and  Lindsay 
(1982). 

2. 1.1.2  Discharge.  Discharge  data  (m3/s)  have  been  obtained  from  the 
Water  Survey  of  Canada  (Environment  Canada  1974)  for  the  network  of 
stations  within  the  AOSERP  study  area  (Figure  4).  The  gauging  stations, 
period  of  record,  and  plots  of  discharge  with  time  are  presented  for 
each  region  to  correspond  with  the  water  quality  sampling  record. 

2. 1.1.3  Water  quality.  The  water  quality  records  used  in  this  study 
were  restricted  to  those  monitoring  programs  conducted  by  Alberta 
Environment  throughout  the  AOSERP  region.  Procedures  for  the 
collection,  preservation,  and  transport  of  water  samples  are  outlined  in 
Akena  (1980).  Most  samples  were  analysed  by  Chemex  Labs  (Alberta)  Ltd. 
under  contracts  with  AOSERP  and  Alberta  Environment.  The  analytical 
methods  are  registered  with  the  National  Water  Quality  Data  Bank' 
(NAQUADAT)  (Appendix  10.1).  Bacteriological  samples  (fecal  and  total 
col i forms)  were  analysed  by  the  Public  Health  Laboratory,  Edmonton. 

Field  data  sheets  for  each  sampling  event  included  on-site 
water  quality  measurements  and  details  on  location  and  depth  at  which 
samples  were  retrieved.  This  information  as  well  as  tabulated  data  may 
be  obtained  from  the  Water  Quality  Control  Branch,  Alberta  Environment. 

2. 1.1.4  Benthic  invertebrates.  Aquatic  invertebrates  are  frequently 
used  as  environmental  indicators  of  water  quality  (McCart  and  Mayhood 
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1980).  Site-specific  studies  frequently  were  instituted  to  determine 
the  effects  of  effluents  on  the  abundance  and  species  composition  of 
invertebrates  by  analysing  samples  above  and  below  point  sources  as  well 
as  before  and  during  development.  Other  invertebrate  studies  were 
conducted  to  provide  baseline  data  in  watersheds  where  government  or 
industry  anticipated  development. 

Most  studies  were  conducted  by  private  consulting  firms  or 
government  agencies  under  contract  to  Alberta  Environment  and  industrial 
proponents.  With  permission,  data  for  analysis  have  been  retrieved  from 
AOSERP  reports. 

2.1.2  Data  Presentation  and  Analysis 

2. 1.2.1  Within  regions.  The  relationships  examined  within  each  region 
depended  on  which  parameters  were  measured  at  the  same  location  and 
time.  For  example,  there  were  more  water  quality  sampling  locations 
than  gauging  stations  in  the  study  area  and,  for  site-specific 
investigations,  there  were  more  biological  sampling  sites  than  water 
quality  sites. 

Tabular  summaries  were  prepared  representing  number  of 
samples,  and  median  and  maximum  values  for  water  quality  parameters. 
Values  representing  the  Alberta  Surface  Water  Objectives  (ASWO)  and 
Canadian  Drinking  Water  Guidelines  (CDWG)  also  are  provided. 

A correlation  analysis  (Sokal  and  Rohlf  1969)  was  performed 
using  median  values  for  water  quality  parameters  (routine,  nutrients, 
and  metals)  for  each  site  with  more  than  10  samples.  Results  of  the 
correlation  analysis  of  water  quality  parameters  were  used  to  select 
those  parameters  to  be  used  in  subsequent  regression  analysis  (Draper 
and  Smith  1966)  relating  water  quality  and  discharge. 

For  those  WSC  sites  with  continuous  discharge  data,  a data  set 
was  prepared  using  all  available  water  quality  data  obtained  at  or  near 
that  locale  and  corresponding  mean  daily  discharge  values  for  the  period 
of  record  during  which  water  quality  values  were  obtained. 
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Water  quality  parameters  vary  directly  (e.g.,  suspended 
solids)  or  indirectly  (e.g.,  specific  conductance)  with  discharge. 
Since  values  for  specific  conductance  and  total  dissolved  solids  were 
frequently  determined  by  calculation  of  ions,  it  seemed  inappropriate  to 
present  relationships  for  ions  which  are  autocorrelated  with 
conductivity. 

When  benthic  invertebrate  data  were  available,  an  analysis  of 
variance  (Sokal  and  Rohlf  1969)  was  used  to  test  if  significant 
differences  occurred  in  mean  abundance  of  invertebrates  among  sampling 
sites.  If  differences  occurred  in  mean  abundance,  the  a_  posteriori 
test,  Student  Newman's  Kuels  procedure  (Sokal  and  Rohlf  1969),  was  used 
to  determine  which  sites  were  similar  to  one  another.  Multiple 
regression  analysis  was  used  to  determine  which  few  (if  any)  of  the 
physical  and  land  formation  parameters  would  be  useful  in  predicting  the 
abundance  or  percentage  of  clean  water  organisms.  Codes  were  assigned 
to  the  various  land  formation  parameters. 

A monitoring  program  of  water  quality  and  biological 
parameters  was  conducted  on  the  Athabasca  mainstem  in  1981  and  1983. 
Results  from  these  studies  are  being  prepared  as  separate  reports  for 
Alberta  Environment. 

Differences  in  study  objectives,  habitat  type,  and  sampling 
bias  resulted  in  the  use  of  a variety  of  invertebrate  sampling 
techniques  in  the  AOSERP  region.  Sampling  differences  which  hinder  data 
comparisons  could  be  alleviated  if  guidelines  were  established  for 
monitoring  studies  in  erosional  and  depositional  habitats. 

2. 1.2.2  Among  regions.  Cluster  analysis  (Wishart  1975)  was  used  to 
determine  grouping  of  similar  sites  for  three  sets  of  water  quality 
parameters  (routine,  nutrients,  and  metals)  at  sites  throughout  the 
AOSERP  study  area.  Dendrograms  were  used  to  display  groupings  among 
sites. 

A principal  component  analysis  (PCA)  was  used  to  reduce  the 
number  of  water  quality  parameters  for  use  in  subsequent  analysis.  The 
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details  of  PCA  are  outlined  by  several  authors  including  Davis  (1973), 
Orloci  (1975),  and  Morrison  (1967).  Green  (1979)  discusses  the 
application  of  the  analyses  to  environmental  problems. 

A PCA  summarizes  the  information  from  a correlation  or 
covariance  (abundance  data)  matrix  in  terms  of  new  components.  The 
first  principal  component  accounts  for  the  greatest  proportion  of  total 
variation  among  sites,  with  successive  components  contributing 
progressi vely  smaller  portions  of  the  total  variation.  The  first  two  or 
three  principal  components  generally  account  for  most  of  the  variability 
in  the  original  data.  The  principal  component  scores  for  each  river 
site  can  be  used  in  subsequent  analysis  (in  place  of  the  original  data) 
with  little  loss  of  information. 

2. 1.2.3  Statistical  analysis.  In  this  study,  data  were  entered  into 
files  on  the  Michigan  Terminal  System  (MTS)  at  the  University  of 
Alberta.  Cluster  analysis  was  performed  using  Ward's  method  of  the 
Clustan  C-l  computer  package  (Wishart  1975).  Analyses  of  variance  and 
covariance  were  performed  using  the  APL  (A  Programing  Language)  library 
of  statistical  packages.  Correlation  and  multiple  regression  analyses 
were  performed  using  the  BMDP  (Biomedical  Computer  Programs-P  Series) 
package  of  computer  programs  (Dixon  and  Brown  1979). 

2.1.3  Impact  of  Development 

Increased  demands  on  surface  waters  associated  with  oil  sands 
development  include  recreational , domestic,  and  commercial  demands  on 
water  quality.  This  study  will  consider  diffuse  and  point  source 
effluents  into  the  Athabasca  mainstem.  Diffuse  inputs  occur  as  a result 
of  runoff  from  various  land  use  activities  or  weathering  of  surficial 
material.  The  effluent  from  the  Fort  McMurray  sewage  treatment  plant, 
saline  discharge  by  Syncrude  Canada  Ltd.  via  Poplar  Creek,  and  waste 
discharge  (industrial  and  lagoons)  from  the  Suncor  oil  sands  facility, 
all  of  which  enter  the  Athabasca  River  from  the  west  bank,  are  the  major 
cultural  point  source  effluents  within  the  AOSERP  study  area.  The 
effects  of  the  Clearwater  River  (the  major  tributary  of  the  Athabasca 
River)  on  the  mainstem  also  will  be  examined. 
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3.  ATHABASCA  MAINSTEM 

The  Athabasca  River  originates  in  the  Columbia  Glacier  of 
Jasper  National  Park  and  flows  unregulated  1450  km  to  Lake  Athabasca  via 
the  delta  distributaries. 

This  section  of  the  report  will  focus  on  the  lower  Athabasca 
River  from  the  Fort  McMurray  area  (upstream  of  the  confluence  with  the 
Horse  River)  to  the  Embarras  Airport,  a distance  of  about  190  km. 

3.1  BASIN  FEATURES 

3.1.1  Land  Formation 

The  Athabasca  River  flows  easterly  to  Fort  McMurray  where  the 
bedrock  geology  changes  from  Cretaceous  rocks  (shale  and  sandstone)  to 
Devonian  limestone  and  the  river  is  directed  northward.  With  the 
exception  of  a 35  km  stretch  of  Cretaceous  rocks  between  the  Ells  and 
Eymundson  rivers,  the  Athabasca  River  flows  over  Devonian  rocks  to  the 
delta.  The  Clearwater  River,  a major  tributary  of  the  Athabasca  River, 
enters  the  mainstem  at  Fort  McMurray. 

The  Athabasca  River  flows  through  the  Methy  Portage  and  Great 
Slave  (Muskeg-Embarras-Athabasca  Delta  ecodistricts)  plains  of  the 
Boreal  Mixedwood  Forest.  The  transition  between  the  southern  Methy 
Portage  and  northern  Great  Slave  plains  occurs  between  the  confluences 
of  the  Steepbank  and  Muskeg  rivers  with  the  mainstem. 

3.1.2  Watershed 

The  Athabasca  River  below  Fort  McMurray  (WSC  gauging  station 

07DA001)  drains  1.33  x 105  km2.  At  Embarras  (WSC  station  07DD001),  the 

5 2 

river  drains  1.52  x 10  km  . The  river  distance  between  the  two 
locations  is  180  km.  The  southern  lowlands  and  northern  plains  of  the 
river  valley  are  moderately  forested. 

3.1.3  Physical  Description 

The  river  reach  upstream  of  Fort  McMurray  is  gorge-like  with 
many  meanders  confined  within  steep  banks  and  rapids  (Grand,  Brule, 
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Bortes,  Cascade,  and  Mountain)  (Bond  and  Berry  1980).  The  Mountain 
Rapids  are  located  2 km  upstream  of  Fort  McMurray.  The  river  bed  is 
composed  of  coarse  gravel,  cobble,  and  boulders. 

There  are  high  clay  bluffs  with  exposed  bitumen  at  Fort 
McMurray  where  the  river  turns  northward.  The  river  gradient  is  gradual 
between  Fort  McMurray  and  Fort  MacKay  (0.144  m/km)  and  further  reduced 
between  Fort  MacKay  and  the  Embarras  Airport  (0.117  m/km)  (Kellerhals 
et  al.  1972). 

The  Athabasca  River  below  Fort  McMurray  has  a mean  channel 
width  of  448  m and  mean  depth  of  1.37  m.  The  channel  is  straight  with 
occasional  islands  and  mid-channel  bars.  The  riverbed  is  sand  with  some 
local  gravel  over  limestone.  Bank  materials  consist  of  silt  and  clay. 

At  the  Embarras  Airport,  the  river  channel  has  a mean  width  of 
384  m and  mean  depth  of  2.65  m.  The  channel  has  irregular  meanders  with 
occasional  islands,  point  bars,  and  mid-channel  bars.  The  riverbed 
consists  of  deep  sand;  the  bank  is  of  sand,  silt,  and  clay. 

3.2  BASELINE  DATA 

3.2.1  Discharge 

The  two  Water  Survey  of  Canada  gauging  stations  located  on  the 
Athabasca  River  within  the  study  area  are  below  Fort  McMurray  (07DA001) 
and  the  Embarras  Airport  (07DD001).  Stage  records  only  are  available  at 
Fort  McMurray  and  at  Shott  Island  upstream  of  the  Fi rebag  River. 

Mean  monthly  discharge  is  presented  for  the  McMurray 
(1976  to  1983)  and  Embarras  (1977  to  1983)  gauging  stations  in  Figure  6. 
This  period  of  record  corresponds  to  the  water  quality  sampling  regime. 

During  the  eight  years  of  record,  the  mean  annual  discharge 
for  the  Athabasca  River  below  Fort  McMurray  ranged  from  481  m3/s  (1981) 
to  773  m3/s  (1978).  With  one  exception  (August  1976),  maximum  mean 
monthly  flows  occurred  in  June  and  July  and  ranged  from  1010  to 
1990  m3/s.  Minimum  mean  monthly  "winter"  flows  ranged  from  107  m3/s 
(1982)  to  208  m3/s  (1980). 


EMBARRAS  AIRPORT  C OTDDOOl  ) 
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Figure  6.  Mean  monthly  discharge  at  Fort  McMurray  and  Embarras  Airport  on 
the  Athabasca  River 
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Mean  annual  long-term  discharge  for  the  Athabasca  River  at 
Embarras  Airport  is  estimated  to  be  766  m3/s  (Kellerhals  et  al . 1972). 
Actual  values  are  unavailable  because  of  interrupted  flow  measurements 
between  November  and  May.  Maximum  mean  monthly  flow  typically  occurred 
in  June  and  July,  with  flows  ranging  from  1170  to  1960  m3/s.  In  1978, 
maximum  flow  occurred  in  September  when  mean  monthly  discharge  was 
1755  m3/s. 

The  Hydrology  Branch,  Technical  Services  Division  of  Alberta 
Environment  (Alberta  Environment  1984)  conducted  a low  flow  analysis  of 
the  Athabasca  River  Basin  based  on  21  years  of  data  using  a water  year 
from  July  1 to  June  30.  The  1:10  year,  7-year  draught  for  the  Athabasca 
River  below  Fort  McMurray  (WSC  Station  07DA001)  is  119  cms.  The  7-day 
low  flows  (119  cms)  were  recorded  in  1961-62,  1967-68,  1970-71  and 
1981-82.  The  1981-82  low  flow  event  of  95.1  cms  for  December  16  to  22 
represented  the  lowest  sequence  for  the  Athabasca  River  below  Fort 
McMurray.  This  draught  represented  a 1 in  100  year  low  flow  event. 

3.2.2  Water  Quality 

Seven  water  quality  locations  were  regularly  monitored  along 
the  Athabasca  River  from  1976  to  1983.  Generally,  samples  were  taken 
from  the  centre  of  the  river  at  the  site  upstream  of  Fort  McMurray 
(CC0012)  and  at  Embarras  (DD0010),  the  furthest  downstream  site.  At  the 
other  five  stations  (DA0203,  DA0205,  DA0206,  DA0207,  and  DA0208), 
samples  were  taken  from  the  left  and  right  banks  of  the  river.  A total 
of  12  sites  was  examined. 

A brief  description  of  the  historical  water  quality  sampling 
stations  is  presented  in  Table  2.  A seasonal  summary  of  water  quality 
constituents  for  the  12  sites  is  presented  in  Appendix  10.2.  Seasons 
were  designated  on  the  basis  of  discharge  regimes:  winter  (Nov.  1 to 

April  15),  spring  (April  16  to  June  30),  summer  (July  1 to  August  31) 
and  Fall  (Sept.  1 to  Oct.  31). 

Since  such  a large  data  set  is  available  for  the  Athabasca 
mainstem,  the  discussion  of  routine,  nutrient,  and  metal  parameters  will 
be  as  follows: 


Table  2.  Description,  period  of  record,  and  number  of  samples  obtained  for 
water  quality  sampling  stations  along  the  Athabasca  River. 


Code 

(OOAT07-) 

Latitude  N 
Longitude  W 
(D  M S) 

Description 

Period  of 
Record 

Number  of 
Samples 

CC0012 

56 

43 

06 

Above  Fort  McMurray 

1976-83 

74 

111 

24 

11 

100  m above  Horse  River 

DA0200 

56 

44 

33 

Off  McDonald  Island 

1976 

15 

111 

23 

25 

DA0201 

56 

46 

48 

Km  10.5 

1976 

7 

111 

24 

18 

DA0202 

56 

58 

01 

Above  Suncor  plant 

1976 

8 

111 

27 

14 

DA0203 

56 

56 

21 

Above  Poplar  Creek 

1976-77 

69 

111 

26 

34 

km  30.6 

1981-83 

DA0204 

57 

01 

15 

At  old  A0SERP  dock 

1976-77 

11 

111 

29 

49 

km  42.3 

DA0205 

57 

04 

33 

Below  Suncor 

1976-77 

51 

111 

31 

59 

km  47.9 

1982-83 

DA0206 

57 

07 

49 

Above  Muskeg  River 

1976-77 

107 

111 

36 

18 

km  55.5 

1980-83 

DA0207 

57 

11 

16 

At  Fort  Mac Kay 

1976-83  ' 

111 

111 

37 

26 

DA0208 

57 

21 

56 

Above  Tar  River 

1976-77 

49 

111 

39 

43 

km  84.3 

1980-83 

continued  . . 
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Table  2.  Concluded. 


Code 

Latitude  N 
Longitude  W 
(D  M S) 

Description 

Period  of 
Record 

Number  of 
Samples 

DA0209 

57 

44 

23 

Above  Fi rebag  River 

1976-77, 

12 

111 

21 

58 

1980 

DA0210 

57 

55 

31 

Below  Fi rebag  River 

1976-78 

8 

111 

26 

41 

DD0010 

58 

12 

18 

At  Embarras  Airport, 

1977-83 

47 

111 

23 

24 

WSC  gauge 
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1.  Those  not  found  above  the  detection  limit; 

2.  Those  found  not  to  vary  longitudinally  or  seasonally 
( i . e . , in  space  or  time); 

3.  Those  exhibiting  seasonal  but  not  longitudinal  changes; 

4.  Those  exhibiting  longitudinal  but  not  seasonal  changes; 
and 

5.  Those  that  change  longitudinally  and  seasonally. 

3. 2. 2.1  Routine  parameters.  Routine  parameters  include  the  major 
inorganic  ions,  parameters  related  to  suspended  particulates,  organic 
constituents  such  as  total  oraanic  carbon,  phenolic  substances,  and  oil 
and  grease.  The  major  ions  reflect  the  influence  of  seasonal  discharge 
changes  on  the  river.  The  ions  and  suspended  substances  occur  in  water 
as  a result  of  natural  erosion  of  bedrock  and  surficial  materials.  Ions 
are  considered  to  be  conservative  parameters  ( i . e . , they  may  be  used  as 
tracers  in  the  river  system  because  they  are  not  affected  by. chemical  or 
biological  processes).  Chloride,  however,  may  vary  in  surface  waters 
because  of  the  use  of  chlorine  in  domestic  or  industrial  processes. 

The  surface  waters  of  the  Athabasca  River  are  characterized  as 
well  buffered  (elevated  pH  and  alkalinity  levels)  and  hard  to  very  hard, 
with  median  hardness  concentrations  ranging  from  84  to  177  mg/L 
(Appendix  10.2,  Tables  27  to  38). 

The  maximum  and  minimum  values  (expressed  as  me/L)  of  major 
ions  are  displayed  in  Figure  7.  Specific  conductance  or  conductivity, 
which  is  a measure  of  the  total  ion  composition,  expresses  the  degree  to 
which  water  conducts  an  electrical  current.  Changes  in  conductivity  are 
a function  of  mineral  concentration  in  the  river  water.  Thus,  increases 
in  conductivity  in  the  downstream  direction  reflect  increases  in  the 
mineral  content  and  sediment  load. 

The  dissolved  salts,  sodium  chloride  and  calcium  bicarbonate, 
vary  indirectly  with  discharge  ( i . e . , seasonally)  and  among  sites. 
Accordingly,  conductivity,  filterable  residue  (dissolved  solids), 
hardness  (a  measure  of  calcium  and  magnesium),  and  alkalinity  (a  measure 
of  carbonate,  bicarbonate,  and  hydroxides)  change  in  the  same  manner. 
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Figure  7. 


Maximum  and  minimum  values  ( me/L)  of 
major  ions  at  sites  along  the  Athabasca 
River  . 
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Although  calcium  bicarbonate  is  the  dominant  dissolved  salt, 
sodium  chloride  is  prevalent  at  several  sites.  Sources  of  sodium 
chloride  include  the  inputs  from  the  Clearwater  watershed,  saline 
groundwater  discharge  from  Syncrude  Canada  Ltd.'s  mining  operation  via 
Poplar  Creek  (upstream  of  site  DA0205),  and  saline  seepage  into  the 
mainstem  (e.g.,  Embarras  site  DD0010).  Elevated  chloride  levels  present 
at  site  DA0203  (right  bank)  reflect  the  chlorine  used  in  treatment  of 
domestic  sewage. 

Since  ion  constituents  vary  inversely  with  flow,  the  highest 
concentrations  are  present  during  low  winter  flows.  Magnesium  and 
sulphate  ions  did  vary  inversely  with  discharge  but  did  not  vary  among 
sites.  Potassium  did  ot  vary  with  discharge  or  among  sites. 

Total  organic  (TOC)  consists  of  dissolved  and  particulate 
organic  fractions.  Most  of  the  TOC  in  water  samples  is  composed  of 
humic  substances  resulting  from  the  degradation  of  plant  and  animal 
matter.  Decomposition  of  carbon  reduces  oxygen  concentrations  in 
surface  waters.  No  surface  water  guidelines  have  been  established  for 
TOC. 

Although  median  and  maximum  values  of  TOC  were  generally  lower 
in  winter  than  in  the  other  seasons,  no  real  differences  in  TOC  could  be 
discerned  at  different  sampling  times  within  the  open  water  season. 
Differences  in  TOC,  however,  occurred  among  sites.  Maximum  winter  TOC 
values  were  encountered  in  the  Fort  MacKay  area  (site  8,  35  mg/L;  site 
9,  63  mg/L).  These  differences  among  sites  were  not  pronounced  in  other 
seasons. 

Phenolic  compounds  at  very  high  concentrations  may  be  toxic  to 
fish.  At  lower  levels,  phenols  may  cause  fish  tainting  and  taste 
problems.  The  Alberta  Surface  Water  Objective  (ASWO)  for  phenols  is 
0.005  mg/L  and  the  Canadian  Drinking  Water  Guideline  (CDWG)  is 
0.002  mg/L. 

With  few  exceptions,  median  values  of  phenol ics  were  below  the 
surface  water  objectives.  Median  values  of  phenolic  compounds  exceeded 
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the  ASWO  at  the  Ells  River  (DA0208,  left  bank)  in  winter  (0.012  mg/L) 
and  at  the  Embarras  Airport  (DD0012)  in  spring  (0.008  mg/L). 

Maximum  values  of  phenols  exceeded  the  objectives  on  several 
occasions.  In  winter,  non-compliance  occurred  at  all  sites  except 
CC0012  (upstream  of  Fort  McMurray).  Maximum  values  were  deduced  to 
0.006  mg/L  at  the  Embarras  Airport.  In  spring  the  maximum  value  of 
phenols  was  recorded  at  the  upstream  control  site  CC0012  (0.041  mg/L). 
Somewhat  elevated  levels  also  were  noted  in  spring  at  site  DA0207  near 
Fort  MacKay  (left  bank,  0.016  mg/L;  right  bank,  0.012)  and  at  the 
Embarras  Airport  (0.011  mg/L).  In  summer,  maximum  phenolic  values 
exceeding  the  objectives  were  noted  at  site  DA0203  below  the  sewage 
treatment  plant  on  the  left  bank  (0.009  mg/L)  and  at  the  barge  tie-up  on 
the  right  bank  (0.015  mg/L).  Levels  were  exceeded  below  Suncor  (DA0205 
left  bank;  0.009  mg/L)  and  at  other  locations  downstream,  but  not  at  the 
Embarras  site.  During  the  fall,  non-compl iant  values  occurred  at 
several  locations  along  the  mainstem,  including  the  site  upstream  of 
Fort  McMurray. 

No  guidelines  have  been  established  for  oil  and  grease. 
Because  these  products  occur  on  the  surface  film,  the  concentrations  do 
not  vary  with  flow.  No  distinctive  pattern  was  noted  among  sites. 

3. 2. 2. 2 Nutrients.  Phosphorus  and  nitrogen  are  essential  nutrients 
which  may  limit  plant  growth. 

Data  representing  two  forms  of  phosphorus  (total  and 
orthophosphate)  are  presented  in  Appendix  10.2  (Tables  39  to  44). 
Because  total  phosphorus  includes  both  dissolved  and  particulate  forms, 
concentrations  vary  directly  with  flow.  The  particulate  or  apatite 
phosphorus  leached  from  rocks  is  generally  not  biologically  available  in 
rivers.  Orthophosphate  phosphorus  (PO^-P)  which  is  present  in 
fertilizers,  storm  runoff,  and  sewage  treatment  wastes,  is  biologically 
available. 

The  ASWO  for  total  phosphorus  is  50  ug/L.  Median  values  of 
total  phosphorus  were  below  this  level  in  water  and  fall  (low  flow 
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periods),  but  exceeded  the  objectives  during  spring  and  summer. 
Orthophosphates  did  not  vary  appreciably  with  season.  In  the  spring, 
maximum  values  were  higher  at  both  banks  of  sites  DA0203,  DA0205, 
DA0206,  and  DA0207  than  at  other  locations  (Appendix  10.2,  Tables  39 
to  44). 

Results  representing  three  forms  of  nitrogen  (nitrate/nitrite, 
ammonia,  and  total  kjeldahl  nitrogen)  are  presented  in  Appendix  10.2 
(Tables  39  to  44).  Inorganic  forms  of  nitrogen  include  ammonia, 
nitrite,  and  nitrate.  Total  kjeldahl  nitrogen  (TKN)  measures  both 
ammonia  and  organic  forms.  Total  nitrogen  is  calculated  by  the  addition 
of  TKN  and  nitrate/nitrite  nitrogen.  The  ASWO  for  total  nitrogen  is 
1.0  mg/L  (1000  ug/L). 

Nitrite  nitrogen  is  readily  converted  to  nitrate  nitrogen  in 
the  presence  of  oxygen.  Of  the  inorganic  forms,  ammonia  is  the  major 
form  discharged  from  sewage  treatment  plants  and  industrial  effluents. 
Sources  of  organic  nitrogen  include  watershed  runoff  and  sewage 
effluents. 

Nitrate/nitrites  vary  seasonally,  with  the  highest  levels 
occurring  in  winter  particularly  below  the  Fort  McMurray  sewage 
treatment  plant  (max.  560  ug/L)  and  in  the  Fort  MacKay  area  (max. 

500  ug/L). 

Although  ammonia  levels  are  reduced  in  the  fall  (median 
values,  6 to  43  ug/L;  maximum  values,  9 to  130  ug/L),  there  are  no 
substantial  differences  in  concentrations  during  other  seasons  ' 
(Appendix  10.2,  Tables  39  to  44).  Differences  in  ammonia  levels  do 
occur  among  sites.  In  winter,  elevated  levels  occur  below  Fort  McMurray 
(max.  level  440  ug/L).  In  spring,  concentrations  are  lower  at  the  Ells 
River  (DA0208)  and  Embarras  Airport  (DD0010)  sites  than  at  other 
locales.  A similar  pattern  occurs  in  summer  along  with  reduced  ammonia 
concentrations  downstream  of  Suncor  (DA0205).  There  was  insufficient 
information  on  TKN  to  determine  if  differences  occurred  among  seasons  or 
sites. 

Bacteriological  indicators  (total  and  fecal  coliforms)  reflect 
biological  activity.  Bacteria  are  ubiquitous  and  may  always  be  observed 
at  low  levels.  Sources  of  bacterial  activity  include  sewage  and 
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industrial  effluents  as  well  as  agricultural  runoff  contaminated  with 
animal  wastes.  Coliform  bacteria  is  an  indicator  of  the  possible 
presence  of  disease-producing  organisms.  Fecal  coliforms  are  less 
likely  to  grow  outside  the  intestinal  tract  of  warm-blooded  animals. 
Objectives  for  outdoor  recreation  that  excludes  direct  contact  with 
water  are  90%  of  the  samples  with  a total  coliform  density  less  than 
5 000  per  100  mL  and  fecal  coliform  density  less  than  1 000  per  100  mL. 
More  stringent  objectives  are  set  for  direct  contact  with  a limited 
number  of  samples  allowed  to  exceed  1 000  per  100  mL  total  coliforms  and 
200  per  100  mL  for  fecal  coliforms. 

Data  for  total  coliforms  and  fecal  coliforms  are  presented  in 
Appendix  10.2  (Tables  45  and  46).  Lower  levels  of  fecal  coliforms  were 
recorded  in  the  fall  than  in  other  seasons.  Relatively  higher  levels 
of  fecal  coliforms  occur  near  settlements  (Fort  McMurray  and  Fort 
MacKay)  in  winter,  spring,  and  summer.  Maximum  recorded  levels  of 
total  coliforms  were  elevated  at  all  sites  during  the  winter 
(540  t.c./lOO  mL),  summer  (540  to  2400  t.c./lOO  mL)  and  fall 
(690  t.c./lOO  mL).  In  spring  elevated  values  ( 2400  t.c./lOO  mL)  were 
noted  at  several  sites  including  the  site  upstream  of  Fort  McMurray. 

3. 2. 2. 3 Metals.  Of  the  suite  of  metals  which  are  typically  analysed, 
nine  exhibited  values  above  the  detection  limit  (aluminum,  arsenic, 
boron,  copper,  iron,  manganese,  nickel,  vanadium,  and  zinc)  Appendix 
10.2  (Tables  47  to  58). 

Aluminum  concentrations  varied  directly  with  discharge,  with 
elevated  levels  being  present  in  spring  and  summer.  Differences  in 
concentrations  occurred  among  sites,  with  the  highest  values  being 
recorded  at  sites  CC0012,  DA0206,  and  DA0207.  Aluminum  is  not 
considered  to  be  a public  hazard  and  guidelines  have  not  been 
established  for  the  metal. 

Median  values  of  arsenic  concentrations  did  not  differ  among 
sites  or  seasons.  Maximum  values  were  higher  (exceeding  the  surface 
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water  objective  of  10  mg/L)  at  sites  CC0012,  DA0206,  and  DA0207  in  the 
spring  than  elsewhere. 

Boron  did  not  change  seasonally  and,  although  concentrations 
were  elevated  at  sites  DA0206  and  DA0207,  all  values  were  below  the 
surface  water  objectives  of  500  ug/L. 

Copper  varied  directly  with  discharge,  with  maximum 
concentrations  exceeding  the  surface  water  objectives  of  20  ug/L  in 
spring  and  summer.  Median  values  of  copper  were  below  the  objective  in 
all  cases. 

Iron  levels  increased  with  increasing  flows.  With  the 
exception  of  a few  median  values  in  winter,  both  maximum  and  median 
values  of  iron  exceeded  the  surface  water  objectives  of  300  ug/L.  The 
guidelines,  however,  were  established  for  aesthetic  rather  than  health 
reasons  as  iron  discolors  clothing  and  stains  plumbing  fixtures. 

Manganese  levels  varied  directly  with  discharge,  exhibiting 
median  and  maximum  values  which  exceeded  the  surface  water  objective 
(50  ug/L)  in  the  spring  and  summer.  Maximum  values  were  non-compl iant 
in  fall  and  winter.  Guidelines  for  manganese  were  established  because 
of  unfavourable  taste  and  staining  properties  at  high  concentrations. 

Nickel  values  were  higher  during  spring  and  summer  than  in 
fall  or  winter.  Although  there  was  no  consistent  longitudinal  trend  in 
nickel  concentrations,  maximum  values  were  consistently  high  at  site 
CC0012,  upstream  of  Fort  McMurray.  Nickel  occurs  in  substantial  amounts 
in  oil  sands  and  its  presence  in  surface  waters  may  be  attributed  to'the 
weathering  and  leaching  of  oil  sands.  In  spring  and  summer,  bitumen 
frequently  can  be  observed  oozing  out  of  the  valley  walls  along  the 
mainstem,  particularly  in  the  high  cliffs  upstream  of  Fort  McMurray. 

Maximum  levels  of  vanadium  recorded  at  each  site  were  higher 
in  spring  and  summer  than  in  fall  or  winter.  Median  values  were 
consistently  low  (1  to  2 ug/L)  in  all  seasons.  Although  vanadium  occurs 
in  oil  sands,  concentrations  in  surface  waters  below  the  Suncor  plant 
were  within  the  range  of  values  recorded  at  other  sites. 

Median  levels  of  zinc  exhibited  higher  values  during  spring 
and  summer  than  in  other  seasons.  Although  median  values  were  below  the 


37 


surface  water  objective  of  50  ug/L,  maximum  values  at  some  sites 
exceeded  the  objective  in  all  four  seasons. 

In  summary,  the  metals  (with  the  exception  of  arsenic  and 
boron)  varied  directly  with  discharge,  indicating  that  concentrations 
are  associated  with  weathering  and  transport  of  bedrock  and  surficial 
material  rather  than  with  point  source  inputs.  No  consistent 
longitudinal  trend  was  noted  in  metal  concentrations.  High  values, 
however,  typically  occurred  upstream  of  Fort  McMurray  (DD0012)  as  well 
as  upstream  of  the  Muskeg  (DA0206)  and  Mac Kay  (DA0207)  rivers. 

3.3  DELINEATION  OF  IMPACTS 

The  numbers  and  locations  of  all  samples  collected  from  the 
Athabasca  River  are  presented  in  Table  3.  Principal  component  analysis 
(PCA),  however,  was  used  for  the  simultaneous  examination  of  water 
quality  parameters  at  those  12  specified  sites  having  sufficient  data 
(Table  4).  The  analysis  was  performed  on  median  and  maximum  values  of 
selected  water  quality  parameters  for  the  winter,  spring,  summer,  and 
fall  data  sets  (1976  to  1983).  Median  values  were  selected  to  represent 
typical  conditions;  maximum  values  were  selected  to  represent  the  worst 
situation. 

The  following  parameters  were  selected  for  use  in  the  PCA 
because  of  their  association  with  diffuse  and  point  source  effluents: 

1.  Specific  conductivity; 

2.  Bicarbonate  ion; 

3.  Chloride  ion; 

4.  Total  organic  carbon  (TOC); 

5.  Phenolic  compounds; 

6.  Oils  and  grease; 

7.  Orthophosphate  phosphorus  (PO^-P); 

8.  Ammonia; 

9.  Fecal  col i forms; 

10.  Nickel  ; and 

11.  Vanadium. 
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Table  3.  Number  and  location  of  samples  collected  at  sites  along  the 
Athabasca  River. 


Code 

(00AT07-) 

Location3 

Total 

Left 

Centre 

Right 

Not  Specified 

CC0012 

3 

34 

15 

22 

74 

DA0203 

30 

1 

31 

7 

69 

DA0205 

23 

0 

21 

7 

51 

DA0206 

46 

13 

41 

7 

107 

DA0207 

50 

8 

35 

18 

111 

DA0208 

20 

0 

21 

8 

49 

DD0010 

4 

20 

6 

17 

47 

dLooking  downstream. 
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Table  4.  Sample  sites  used  in  data  analysis. 


Site 

Code  (00AT07-) 

Location8 

1 

CC0012 

Centre 

2 

DA0203 

Left 

3 

DA0203 

Right 

4 

DA0205 

Left 

5 

DA0205 

Right 

6 

DA0206 

Left 

7 

DA0206 

Right 

8 

DA0207 

Left 

9 

DA0207 

Right 

10 

DA0208 

Left 

11 

DA0208 

Right 

12 

DD0010 

Centre 

a 

Looking  downstream. 


40 


The  ionic  parameters  (numbers  1,  2,  and  3)  which  vary 
inversely  with  discharge  reflect  changes  in  ion  concentration  between 
tributaries  and  the  mainstem  as  well  as  industrial  saline  discharge  from 
Poplar  Creek  during  the  open  water  season.  Parameter  4 (TOC)  is  a 
measure  of  biological  degradation.  Parameters  5 (phenols)  and  6 (oils 
and  grease)  can  occur  naturally  within  the  system,  but  elevated  levels 
are  often  associated  with  domestic  and  industrial  wastes.  Parameters  7 
to  9,  the  nutrients  and  fecal  coliforms,  are  associated  with  sewage 
effluents.  Although  total  kjeldahl  nitrogen  would  typically  be  used  as 
a measure  of  sewage  effluents,  there  was  not  sufficient  information 
available  to  use  this  parameter  in  the  analysis.  Instead,  ammonia,  the 
major  inorganic  form  of  nitrogen  discharged  from  sewage  plants,  was 
selected.  The  metals  (no.  10,  nickel  and  no.  11,  vanadium)  are 
associated  with  oil  sands. 

The  factor  scores  for  all  sites  are  presented  on  bivariate 
plots  of  the  first  and  second  principal  components.  Designated  site 
clusters  (i.e.,  those  sites  exhibiting  similar  water  quality  parameters) 
are  delineated  on  schematic  maps  of  the  river. 

3.3.1  Winter  (Median  Values) 

The  first  two  principal  components  (PC)  of  median  values  of 
selected  water  quality  parameters  in  the  winter  accounted  for  51.5%  of 
the  total  variance.  The  first  component  alone  accounted  for  30.6%  of 
the  variance.  To  interpret  the  results,  the  coefficients  (loadings)  of 
the  11  parameters  were  examined.  The  first  PC  is  interpreted  as  a 
comparison  of  bicarbonate  ions  and  fecal  coliforms  (with  high  positive 
coefficients)  to  PO^-P  (lowest  coefficient).  For  PC2,  chlorides  and 
ammonia  represent  the  highest  positive  loadings,  whereas  TOC  represents 
the  lowest  negative  loadings. 

The  cluster  of  sites  are  represented  as  five  groups  on  the 
bivariate  plot  and  corresponding  river  map  (Figures  8 and  9).  Sites  2 
(DA0203  left)  and  9 (DA0207  right),  which  are  grouped  together. 
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Principal  component  analysis  of  sample  sites  on 
the  Athabasca  River. 
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exhibited  highest  values  of  bicarbonate  ions  and  fecal  coliforms.  The 
sites  are  associated  with  the  Fort  McMurray  sewage  effluent  and  the  Fort 
MacKay  settlement. 

Cluster  2 (sites  4,  10,  and  11)  represents  an  intermediate 
recovery  zone  downstream  from  the  areas  of  enrichment.  Cluster  3 (sites 
6,  7,  8,  and  12)  is  interpreted  to  represent  recovered,  well -mixed  areas 
of  the  river. 

Sites  3 and  5,  along  the  true  right  bank  of  the  Athabasca 
River  downstream  from  the  Clearwater  River,  form  a distinct  group 
characterized  by  elevated  levels  of  PO^-P  and  low  levels  of  bicarbonate 
ions  and  fecal  coliforms. 

Site  1 (upstream  of  Fort  McMurray)  with  elevated  levels  of  TOC 
is  distinct  from  the  other  11  sites. 

3.3.2  Winter  (Maximum  Values) 

The  first  two  principal  components  resulting  from  the  analysis 
of  maximum  values  in  winter  accounted  for  58.7%  of  the  total  variance. 
PCI  contributed  35.9%  of  the  total  variability.  For  PCI,  conductivity 
(along  with  chloride)  and  nickel  exhibited  the  highest  and  lowest 
coefficients,  respectively.  For  PC2,  TOC  and  oil/grease  displayed  the 
highest  and  lowest  coefficients,  respecti vely . 

The  two  major  clusters,  numbers  1 (sites  2,  4,  6,  and  11)  and 
2 (sites  3,  5,  7,  and  10),  were  grouped  at  the  intersection  of  PCI  and 
PC2  (Figure  10).  The  other  four  sites  exhibited  distinct  features. 
Interestingly,  nickel  concentrations  were  highest  at  site  1 (upstream  of 
Fort  McMurray)  (Figure  10). 

3.3.3  Spring  (Median  Values) 

Together,  PCI  (35.7%)  and  PC2  (21.7%)  explained  57.4%  of  the 
total  variability  of  the  PCA  using  median  water  quality  parameters. 
Nickel  and  chloride  exhibited  the  highest  and  lowest  coefficients  for 
PCI.  Conductivity  (along  with  PO^-P)  and  chloride  had  the  highest  and 
lowest  loading  coefficients  for  PC2. 
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The  sampling  sites  differentiated  into  three  major  clusters 
(Figure  11).  Cluster  1 sites  (1,  6,  8,  9,  and  12)  were  characterized  by 
elevated  concentrations  of  nickel  and  low  levels  of  chloride. 

Cluster  2 consisted  of  two  subgroups.  One  subgroup  (sites  3, 
5,  and  7),  located  along  the  true  right  bank  downstream  of  the 
Clearwater  River,  exhibited  high  chloride  values.  The  other  subgroup 
(sites  10  and  11)  had  lower  levels  of  chloride. 

Cluster  3 included  site  2 (below  the  Fort  McMurray  sewage 
treatment  plant)  and  site  4 (downstream  from  Suncor).  Water  samples 
from  these  locales  were  high  in  conductivity  and  PO^-P. 

3.3.4  Spring  (Maximum  Values) 

PCI  (43.7%)  and  PC2  (24.82%)  accounted  for  68.5%  of  the  total 
variability  in  the  data  set,  representing  maximum  parameter  values  in 
spring. 

Three  major  clusters  of  sample  sites  were  evident  (Figure  11). 
Cluster  1 (sites  2,  3,  and  4)  exhibited  maximum  levels  of  conductivity 
and  PO^-P,  reflecting  the  inputs  from  the  Clearwater  River  and  Fort 
McMurray  sewage  plant.  The  highest  PO^-P  levels  throughout  the  year 
were  recorded  in  the  spring.  Cluster  2 represents  sites  with  high 
concentrations  of  chloride  and  phenolic  material  (sites  5,  10,  11,  and 
12).  Cluster  3 sites  (6,  7,  8 and  9)  exhibited  high  levels  of  TOC  and 
nickel.  Nickel  concentrations  were  higher  in  the  spring  than  in  any 
other  season.  Site  1 was  characterized  by  the  high  concentrations  of 
phenols,  TOC,  and  nickel. 

3.3.5  Summer  (Median  Values) 

PCI  (32.8%)  and  PC2  (16.4%)  accounted  for  about  one  half  of 
the  total  variability  in  the  data  set. 

Site  3 was  distinct  from  other  site  groupings  because  of  the 
high  concentrations  of  chloride  ions  (Figure  12).  The  remaining  11 
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Figure  II.  Principal  component  analysis  of  sample  sites  on 
the  Athabasca  River. 
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Figure  12.  Principal  component  analysis  of  sample  sites  on 
the  Athabasca  River. 
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sites  formed  three  clusters.  Cluster  1 (sites  1,  4,  and  8)  exhibited 
elevated  concentrations  of  vanadium.  Within  the  cluster,  site  8 had 
characteristically  high  conductivity  levels  and  site  4 had  high  TOC 
concentrations.  Cluster  2 (sites  2,  5,  6,  and  12)  was  characterized  by 
elevated  levels  of  TOC.  Cluster  3 sites  (7,  9,  10,  and  11)  had 

relatively  high  levels  of  conductivity. 

3.3.6  Summer  (Maximum  Values) 

PCI  (27.6%)  and  PC2  (23.1%)  accounted  for  about  51%  of  the 
total  variability  in  the  data  set.  Chloride  and  metals  (nickel  and 
vanadium)  values  differentiated  site  groupings  along  PCI  in  a similar 
manner  to  the  spring  (median  value)  data  set.  In  contrast,  site  12  (the 
Embarras  Airport  locale)  rather  than  site  3 exhibited  the  highest 
chloride  values.  Site  3 (elevated  fecal  coliforms)  and  site  8 (high 
bicarbonate  and  metal  concentrations)  also  were  distinct  from  other 
groups  (Figure  12). 

The  largest  cluster  of  five  sites  (1,  2,  4,  6,  and  9)  had 
characteristically  high  concentrations  of  nickel  and  vanadium.  Metal 
concentrations,  however,  were  lower  in  the  summer  than  spring  as 
discharge  levels  begin  to  decline. 

3.3.7  Fall  (Median  Values) 

Because  all  sites  exhibited  the  same  low  value  of  vanadium 
(1  ug/L),  the  parameter  was  deleted  from  the  data  set  for  the  analysis. 

The  first  two  PCs  accounted  for  62.6%  of  the  variability  in 
the  data.  Bicarbonate  ions  and  PO^-P  differentiated  site  clusters  along 
the  first  principal  component  (accounting  for  32.4%  of  the  variability). 
Conductivity  and  nickel  concentrations  distinguished  sites  along  the 
second  component. 

Site  5 was  distinguished  from  all  other  sites  because  of 
relatively  high  levels  of  PO^-P  (Figure  13).  The  three  clusters  of  the 
remaining  sites  formed  a patchwork  of  interlocking  sites  along  the 
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Figure  13.  Principal  component  analysis  of  sample  sites  on 
the  Athabasca  River. 
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length  of  the  mainstem  (Figure  9).  For  example,  sites  1,  8,  and  12  were 
grouped  together  because  concentrations  of  nickel  at  these  sites  were 
greater  than  1 mg/L.  Left  and  right  bank  pairs  (sites  2/3  and  10/11) 
were  grouped  because  of  similarities  in  conductivity  levels. 

3.3.8  Fall  (Maximum  Values) 

All  11  parameters  were  used  in  the  PCA  of  the  fall  (maximum 
values)  data  set.  PCI  (33.2%)  and  PC2  (27.6%)  accounted  for  about  61% 
of  the  variability  in  the  data  set. 

Conductivity,  bicarbonate  ions,  and  PO^-P  were  useful  to  a 
limited  extent  in  isolating  clusters  of  sample  sets  (Figure  13).  Eight 
sites  were  incorporated  into  one  loose  cluster  on  the  basis  of 
conductivity.  Sites  2,  4,  and  5 (located  downstream  from  the  Fort 
McMurray  sewage  treatment  plant)  formed  one  group  characterized  by 
elevated  levels  of  PO^-P.  Site  1,  exhibiting  high  concentrations  of  TOC 
and  PO^-P,  was  isolated  from  other  river  sites. 

3.3.9  Summary 

The  following  items  summarize  the  results  of  the  principal 
component  analyses  of  impact  related  water  quality  parameters  at  12 
sites  along  the  Athabasca  River. 

1.  There  was  a consistent  seasonal  pattern  in  the  parameters 
which  distinguished  sites  along  the  first  component. 
During  low  flow  periods  (winter  and  fall)  bicarbonate 
ions  and  PO^-P  were  the  parameters  that  differentiated 
sample  sites.  In  addition,  concentrations  of  nickel  were 
useful  in  distinguishing  sites  in  winter  (maximum 
values).  During  high  flow  periods  (spring  and  summer), 
TOC,  chloride  ions,  nickel,  and  vanadium  differentiated 
groups  of  sample  sites. 

2.  Change  in  ion  concentration  was  attributed  to  inputs  from  - 
either  the  Clearwater  River  or  natural  saline  seepage 
rather  than  the  industrial  saline  discharge  via  Poplar 
Creek. 
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3.  Parameters  which  serve  as  markers  for  sewage  effluents 
are  PO^-P,  TOC  and  fecal  conforms.  Municipal  wastes  at 
Fort  McMurray  were  continuously  discharged.  The 
downstream  effects  of  treated  sewage  were  limited  to  the 
mainstem  upstream  of  the  confluence  with  the  Muskeg 
River. 

4.  The  metals  nickel  and  vanadium  are  associated  with  oil 
sands  in  the  AOSERP  region.  The  concentrations  of  nickel 
and  vanadium  were  higher  at  periods  of  increased  flow 
(spring  and  summer),  indicating  the  significant  effect  of 
weathering  of  rocks.  Sites  exhibiting  elevated 
concentrations  of  nickel  and  vanadium  in  surface  waters 
were  associated  with  the  natural  occurrence  of  bitumen  in 
the  river  bank  and  valley  wall  (e.g.,  CC0012)  rather  than 
with  industrial  effluents  (site  DA0205  left  bank).  Some 
metals,  however,  which  exhibit  very  low  concentrations  in 
water  may  be  associated  with  sediments. 
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4.  THE  NORTH 

The  northern  region  includes  the  drainage  basins  within  the 
Athabasca  Delta.  The  delta  includes  the  following  three  subregions: 

(1)  the  distributaries  (Embarras  River,  Fletcher  Channel,  Goose  Island 
Channel,  and  Big  Point  Channel)  which  enter  Lake  Athabasca  from  the 
Athabasca  River;  (2)  outflow  from  Lakes  Claire  and  Mamawi  which 
originate  from  the  Birch  River  system;  and  (3)  Lake  Athabasca  outflows 
of  Chenal  des  Quatre  Fourches  and  Riviere  des  Rochers  which  join  the 
Peace  River  to  form  the  Slave  River. 

4.1  BASIN  FEATURES 

4.1.1  Land  Formations 

Land  formations  of  the  Athabasca  Delta  region  have  been 
discussed  by  Bayrock  (1971)  and  Turchenek  and  Lindsay  (1982). 

The  delta  was  formed  about  10  000  years  ago  during  the 
recession  of  the  Pleistocene  ice  age.  Although  the  delta  now  separates 
Lakes  Claire  and  Athabasca,  the  two  lakes  were  formerly  one  basin. 

All  Athabasca  River  sediments  are  deposited  in  the  delta. 
Sedimentation  from  streams  entering  Lake  Claire  constitutes  less  than 
15%  of  the  delta  sediments.  Coarse  sands  representing  the  bottom  load 
are  deposited  first.  The  sands  are  subsequently  covered  by  fine 
sediments  which  settle  out  of  the  water  column. 

Most  of  the  region,  with  the  exception  of  the  Athabasca 
distributaries,  occurs  within  the  Northern  Boreal  Forest.  The 
distributaries  lie  within  the  Boreal  Mixedwood  Forest.  Riparian 
forests,  deciduous  shrubs,  and  fen  communities  are  prevalent  throughout 
the  delta. 

The  bedrock  geology  of  the  north  is  characterized  by  granite 
rocks  of  the  Archean  in  the  east  and  Devonian  rocks  (limestone  and 
shale)  in  the  west  (Figure  14).  The  Archean  constitutes  the  earlier  or 
older  rocks  of  the  Precambrian  era.  The  line  separating  the  two 
geological  types  passes  through  Mamawi  Lake. 
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Figure  14.  Bedrock  geology  of  the  delta  area 

( Bosed  on  Bayrock  , 1971  ) 
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4.1.2  Watershed 

The  Athabasca  River  enters  Lake  Athabasca  through  a system  of 
channels  including  the  Embarras  River,  Fletcher  Channel,  Goose  Island 
Channel,  and  Big  Point  Channel  (Figure  15).  Flow  and  sediment 
deposition  increase  in  channels  from  west  to  east,  resulting  in  a 
migration  of  river  outflow  to  the  east. 

Bayrock  (1971)  forecasted  the  breakthrough  of  the  Athabasca 
River  into  the  Embarras  River  with  flow  ultimately  being  redirected  into 
Lakes  Mamawi  and  Claire.  To  prevent  this  flow  transfer,  a cutoff  was 
constructed  near  Embarras  Portage  in  1972. 

Flow  from  the  Birch  River  system  through  Lakes  Claire  (area, 
1160  km2;  volume,  1.2  km3)  and  Mamawi  (area,  160  km2;  volume,  0.16  km3) 
enters  the  Peace  River  (via  Chenal  des  Quatre  Fourches)  and  Lake 
Athabasca  (area,  7850  km2;  estimated  volume,  200  km3).  Reverse  flows 
can  occur  in  the  Chenal  des  Quatre  Fourches  bringing  Peace  River  water 
into  Lake  Mamawi. 

Most  of  the  Lake  Athabasca  outflow  occurs  through  Riviere  des 
Rochers  and  Chenal  des  Quatre  Fourches  (Neill  et  al . 1981).  Revillon 
Coupe  is  a sub-distributary  of  the  Riviere  des  Rochers. 

4.1.3  Physical  Description  of  Water  Bodies 

Neill  et  al . (1981)  summarized  hydraulic  conditions  of  the 
channels  and  lakes  of  the  delta  using  satellite  imagery  and  colour 
aerial  photography. 

Operation  of  the  Bennett  Dam,  located  on  the  Peace  River  in 
British  Columbia,  has  altered  the  physical  features  of  the  Peace 
Athabasca  Delta  by  controlling  water  levels  and  flow  regimes.  The 
Bennett  Dam  was  constructed  in  1967  and  the  filling  of  the  reservoir, 
Willi ston  Lake,  was  completed  in  1972,  witholding  a volume  of  60  km3 
from  the  Peace  River  in  Alberta.  Prior  to  dam  closure,  water  levels  in 
the  Peace  River  at  the  confluence  of  the  Athabasca  outflow  were  high 
enough  to  restrict  or  reverse  flow,  thereby  flooding  low  lying  areas, 
and  thus  enhancing  nutrients  and  organic-  levels.  The  control  structure 
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Figure  15.  Delta  area  surface  water  drainage 

( Adapted  from  Neill  et  al.  , 1961  ) 
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reduced  water  levels  about  3 metres  at  the  confluence  of  the  Peace  River 
and  Riviere  des  Rochers.  Annual  maximum  water  levels  in  Lake  Athabasca 
were  reduced  about  1.2  metres. 

Mitigation  measures  to  maintain  water  levels  included 
construction  of  weirs  on  delta  channels.  A temporary  dam  was 
constructed  in  1971  on  the  soutwest  arm  of  Chenal  des  Quatre  Fourches, 
the  control  outlet  of  Lake  Mamawi . The  weir,  which  was  designed  to 
maintain  water  levels  in  Lakes  Claire  and  Mamawi,  was  frequently  washed 
out  and  rebuilt  until  the  project  was  abandoned  in  December  1975.  Water 
levels,  however,  have  remained  the  same  in  Lakes  Claire,  Mamawi,  and 
Athabasca  because  slopes  of  the  connecting  channels  are  flat.  In  June 
1975  and  March  1976,  permanent  weirs  were  established  on  Riviere  des 
Rochers  ana  Revillon  Coupe  to  control  the  outflow  from  Lake  Athabasca 
when  water  levels  on  the  Peace  River  are  low. 

4.2  BASELINE  DATA 

4.2.1  Discharge 

There  are  no  continuous  gauging  stations  on  the  Athabasca 
River  distributaries,  although  infrequent  flow  measurements  have  been 
taken  by  the  Water  Survey  of  Canada.  Most  (75%)  of  the  flow  from  the 
Athabasca  mainstem  enters  the  delta  via  the  Goose  Island  and  Big  Point 
channels  (Figure  16)  (Neill  et  al . 1981). 

The  average  annual  outflow  from  Lake  Athabasca  is  about 
45  km3.  The  river  systems  contributing  to  this  outflow  are  the 
Athabasca  (52%),  Fond  du  Lac  (21%),  and  small  tributaries  (21%) 
(Figure  5)  (Neill  and  Evans  1979).  Outflows  from  Lakes  Claire  and 
Mamawi,  which  for  the  most  part  bypass  Lake  Athabasca,  contribute  less 
than  3%  of  the  average  annual  Lake  Athabasca  outflow  (Neill  et  al. 
1981). 

4.2.2  Water  Quality 

Water  quality  of  lakes  and  channels  in  the  northern  region  is 
likely  influenced  by  the  control  structure  on  the  Peace  River  and  by 
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LAKE  ATHABASCA 


Figure  16.  Approximate  average  distribution  of  Athabasca 

River  flow  among  four  main  delta  distributaries  . 

( Adapted  from  Neill  et  al.  , 1981  ) 


B'9  Point  channel 
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differing  geology.  Water  quality  values  also  may  differ  between  lake- 
draining and  ri ver-draining  channels.  A description  of  the  water 
quality  sampling  sites  is  presented  in  Table  5. 

4.2.2. 1 Routine  parameters.  A summary  of  routine  parameters  for  the 
region  is  presented  in  Appendix  10.3  (Table  59). 

Surface  waters  of  the  northern  region  are  neutral  or  alkaline. 
Lowest  pH  values  were  recorded  in  waters  overlying  granite  rocks, 
including  Lake  Athabasca  (6.7),  Riviere  des  Rochers  (6.8),  and  Jackfish 
Creek  (7.0).  These  water  bodies  also  exhibited  the  lowest  buffering 
capacity  in  the  region  with  median  values  ranging  from  33  to  70  mg/L. 
Median  values  of  total  alkalinity  for  other  waterbodies  ranged  from 
99  to  118  mg/L.  Water  samples  from  the  devonian  (limestone)  area 
exhibited  higher  median  values  of  hardness  (147  to  162  mg/L)  than  waters 
overlying  granite  (35  to  113  mg/L).  Similar  trends  in  the  relationship 
between  water  quality  constituents  and  geology  (i.e.,  lower  values  in 
waters  overlying  granite  than  limestone)  were  noted  for  turbidity, 
colour,  filterable  residue,  and  conductivity  (Appendix  10.3,  Table  59). 

The  relative  contribution  of  major  ions  is  presented  in 
Figure  17.  The  total  ion  concentration  exhibited  by  sites  located  on 
granite  is  lower  than  sites  located  in  limestone  areas.  The  relative 
contribution  of  ions  also  differs  between  geological  types.  In  granite 
areas,  calcium  bicarbonate  is  the  dominant  dissolved  salt  in  drainage 
systems.  In  limestone  areas,  both  sodium  chloride  and  calcium' 
bicarbonate  salts  are  prevalent.  Lake  Claire  also  exhibits  elevated 
levels  of  sulphate  ions. 

Because  soil  layers  overlying  limestone  are  thicker  than 
layers  overlying  granite,  values  of  total  organic  content  are 
correspondingly  higher  at  sample  sites  on  the  devonian  than  the  Archean 
geological  type  (Appendix  10.3,  Table  59). 

Median  values  of  phenolic  materials  were  below  the  Alberta 
Surface  Water  Objective  (ASW0)  of  0.005  mg/L.  Maximum  values  exceeded 
the  objectives  at  Big  Point  Channel  (0.032),  Chenal  des  Quatre  Fourches 
(0.019  mg/L),  and  Riviere  des  Rochers  (0.022  mg/L). 


Table  5.  Description,  period  of  record,  and  number  of  samples  obtained  for  water  quality  sampling 
stations  within  the  northern  region.  See  Figure  4 for  location  of  sites. 
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Figure  17 


Maximum  and  minimum  values  (me/L)  of  major 
ions  in  the  northern  region  . 


62 


No  guidelines  have  been  established  for  oil  and  grease. 
Median  values  ranged  from  0.5  to  1.2  mg/L;  maximum  values  ranged  from 
1.3  to  3.8  mg/L. 

4. 2. 2. 2 Nutrients.  Phosphorus  and  nitrogen  are  the  major  nutrients 
controlling  plant  growth. 

Total  phosphorus  (particulate  and  dissolved  fractions)  and 
orthophosphate  phosphorus  (PO^-P)  are  the  two  forms  which  are  presented 
in  Appendix  10.3  (Table  60).  Since  total  phosphorus  varies  directly 
with  discharge,  values  are  elevated  in  spring  and  summer  during  high 
runoff  periods.  Median  values  of  total  phosphorus  were  below  the  ASWO 
of  50  ug/L  in  Lake  Athabasca,  Jackfish  Creek,  and  Riviere  des  Rochers. 
Median  values  at  the  other  sites  ranged  from  51  to  80  ug/L.  Maximum 
values  ranged  from  40  to  930  ug/L. 

Orthophosphates  are  a measure  of  natural  or  cultural 
enrichment.  Median  values  ranged  from  4 to  14  ug/L.  Maximum  levels  of 
PO^-P  in  Lake  Athabasca  (9  ug/L)  and  Richardson  Lake  (19  ug/L)  in 
granite  areas  were  substantially  lower  than  in  Mamawi  Lake  Channel 
(280  ug/L)  and  Lake  Claire  (180  ug/L)  in  limestone  areas.  Maximum  P0^-P 
levels  in  Big  Point  Channel  (235  ug/L)  and  Chenal  des  Quatre  Fourches 
(360  ug/L)  were  higher  than  maximum  levels  in  other  channels  (18  to 
70  ug/L). 

The  ASWO  for  total  nitrogen  (total  kjeldahl  nitrogen  and 
nitrate/nitrite)  is  1.0  mg/L  (1000  ug/L).  Values  for  total  kjeldahl 
nitrogen  (TKN  is  a measure  of  organic  nitrogen  and  ammonia),  ammonia, 
and  nitrate/nitrite  nitrogen  are  presented  in  Appendix  10.3  (Table  60). 

Maximum  values  of  TKN  for  all  water  bodies  in  the  region 
ranged  from  1.20  to  9.15  mg/L,  exceeding  the  ASWO  for  total  nitrogen. 
Elevated  median  values  were  recorded  in  water  bodies  of  the  limestone 
region  (1.06  to  1.62  mg/L)  as  well  as  in  Richardson  Lake  (1.40  mg/L). 

Ammonia  and  nitrate/nitrite  nitrogen  are  the  two  inorganic 
forms  of  nitrogen.  Median  values  of  ammonia  ranged  from  0.03  to 
0.06  mg/L.  The  highest  ammonia  value  was  recorded  in  Big  Point  Channel 
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(1.10  mg/L).  Median  values  of  nitrate/nitrite  nitrogen  ranged  from  less 
than  0.01  to  0.07  mg/L.  Maximum  values  ranged  from  0.08  to  0.80  mg/L, 
with  the  highest  value  recorded  in  Riviere  des  Rochers. 

4. 2. 2. 3 Metals.  Summary  data  for  16  metals  obtained  from  the 
monitoring  program  in  the  northern  region  are  presented  in  Appendix  10.3 
(Table  61).  Surface  water  objectives  have  been  established  for  12  of 
the  16  metals.  Objectives  have  not  been  established  for  aluminum, 
cobalt,  nickel,  or  vanadium  because  they  are  usually  of  no  public  health 
significance  (McNeely  et  al . 1979). 

Median  and  maximum  values  were  below  the  ASW0  in  all 
waterbodies  for  arsenic  (10  ug/L),  boron  (500  ug/L),  cadmium  (10  ug/L), 
chromium  (50  ug/L),  lead  (50  ug/L),  selenium  (10  ug/L),  and  silver 
(50  ug/L). 

Median  values  of  aluminum  at  all  sites  ranged  from  120  to 
645  ug/L  and  so  were  within  the  typical  surface  water  levels  of  less 
than  1.0  mg/L  (1000  ug/L)  (McNeely  et  al.  1979).  The  highest  value  of 
aluminum  (30.4  mg/L)  was  recorded  at  Prairie  Creek.  Other  maximum 
values  ranged  from  440  to  6900  ug/L). 

The  presence  of  cobalt  in  surface  waters  results  from 
weathering  of  surficial  material.  Median  values  of  cobalt  for  all 
sample  sites  were  less  than  the  detection  limit  of  2 ug/L.  Maximum 
values  ranged  from  less  than  2 to  4 ug/L. 

Although  copper  is  an  essential  element  for  plants  and 
animals,  elevated  levels  of  copper  in  water  are  distasteful.  The 
threshold  levels  for  taste  perception  range  from  1 to  5 mg/L  (McNeely  et 
al . 1979).  The  highest  value  of  copper  was  recorded  in  samples  from  Big 
Point  Channel  (150  ug/L)  and  Rivere  des  Rochers  (260  ug/L). 

The  ASWOs  for  iron  (300  ug/L)  and  manganese  (50  ug/L)  were 
established  for  aesthetic  reasons  since  elevated  levels  of  these 
constituents  can  alter  the  taste  of  water  as  well  as  stain  clothes  and 
plumbing  fixtures.  Median  and  maximum  values  of  iron  exceeded  the 
objective  at  all  sample  sites.  Although  median  values  of  manganese  for 
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all  sites  were  within  the  ASWO,  maximum  values  exceeded  the  objective  at 
all  sites. 

Non-compliance  with  surface  water  objectives  was  also  noted 
for  mercury.  All  sites  (except  Richardson  Lake)  exceeded  the  objective 
of  0.1  ug/L. 

Nickel  and  vanadium  may  occur  in  surface  waters  as  a result  of 
weathering  of  surficial  material  or  as  wastes  from  industrial  processing 
such  as  fossil  fuels.  Median  values  of  nickel  ranged  from  less  than 
2 to  4 ug/L.  A sample  from  Big  Point  Channel  exhibited  the  highest 
nickel  levels  (30  ug/L).  Median  values  of  vanadium  were  less  than  or 
equal  to  1 ug/L  at  all  sites.  The  highest  value  of  vanadium  was 
recorded  at  Big  Point  Channel  (20  ug/L)  and  Chenal  des  Quatres  Fourches 
(32  ug/L). 

The  median  values  of  zinc  at  all  sites  were  below  the  ASWO  of 
50  ug/L.  With  the  exception  of  Prairie  River,  maximum  values  of  sample 
sites  exceeded  the  objectives.  Although  zinc  is  not  toxic  to  man, 
elevated  levels  (5  mg/L)  will  alter  the  taste  of  water  (McNeely  et  al . 
1979). 

4.2.3  Benthic  Invertebrates 

Gallup  et  al.  (1971)  conducted  a survey  of  planktonic  and 
benthic  invertebrates  of  the  Peace  Athabasca  Delta  during  the  summer  of 
1971.  Earlier,  Rawson  (1947)  studied  the  invertebrates  of  Lake 
Athabasca. 

Gallup  et  al.  (1971)  collected  benthic  invertebrates  with  an 
Ekman  grab,  the  contents  of  which  were  washed  through  a 600  micron 
sieve.  Highest  standing  crops  (more  than  500  mg  wet  weight/sample)  of 
invertebrates  were  recorded  in  the  lower  reaches  of  rivers  entering 
lakes  (Birch  River  and  Mamawi  Creek)  and  channels  between  Lakes  Claire 
and  Mamawi.  A mean  standing  crop  of  100  mg  wet  weight/sample  was 
recorded  for  the  other  sites.  Production  was  greater  near  shorelines 
than  in  open  water  areas  of  lakes. 

Pontoporeia  hoyi , an  amphipod,  dominated  the  bottom  fauna  of 
the  western  end  of  Lake  Athabasca  in  numbers  (61%)  and  wet  weight  (70%) 
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(Gallup  et  al.  1971).  Clams  (Sphaeridae) , snails  (Lymnaeidae)  and 
midges  dominated  the  fauna  in  other  water  bodies. 

4.3  RELATIONSHIPS  AMONG  PARAMETERS 

Because  only  stage  (water  level)  data  are  available  at  sites 
in  the  delta,  relationships  among  parameters  focussed  on  water  quality 
constituents  alone.  A correlation  analysis  was  performed  between 
conductivity  (a  known  correlate  of  discharge)  and  other  water  quality 
constituents  using  median  values  from  each  sampling  site  (Appendix  10.3, 
Tables  59  to  61).  Results  showed  that  six  constituents  (total  organic 
carbon,  total  phosphorus,  total  kjeldahl  nitrogen,  aluminum,  copper,  and 
iron)  were  significantly  correlated  with  conductivity  at  the  0.01  level. 

Relationships  between  each  of  the  six  constituents  and 
conductivity  were  estimated  using  all  available  data  from  sample  sites 
representing  water  bodies  entering  Lake  Athabasca  from  southern 
distributaries  (Jackfish  Creek,  Big  Point  Channel),  channel  inputs  from 
lentic  areas  in  the  west  (Prairie  Creek,  Mamawi  Lake  Channel),  and 
channels  leaving  the  northwest  end  of  Lake  Athabasca  (Chenal  des  Quatre 
Fourches,  Riviere  des  Rochers).  Of  the  six  sample  sites,  one  (Big  Point 
Channel)  receives  river  water  and  the  other  five  receive  drainage  from 
lakes. 

Regression  analyses  were  performed  using  linear,  log/linear, 
and  log/log  transformations  of  the  data.  Equations  from  those 
transformations  yielding  the  best  fit  (i.e.,  least  variability  and 
highest  R2)  were  selected  and  presented  in  Table  6. 

There  were  significant  relationships  between  TOC  and 
conductivity  for  Jackfish  Creek,  Prairie  Creek,  and  Big  Point  Channel 
(Table  6).  The  TOC  varies  inversely  with  conductivity  in  channels  which 
drain  lakes  (Jackfish  and  Prairie  creeks),  but  directly  with  Big  Point 
Channel,  a distributary  of  the  Athabasca  River.  These  relationships 
suggest  that  there  are  seasonal  differences  in  biological  activity 
(indicated  by  TOC)  between  rivers  and  lakes. 

There  was  a significant  log/log  positive  relationship  between 
total  phosphorus  and  conductivity  in  Big  Point  Channel.  No  significant 
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Table  6.  Summary  of  regression  equations  of  total  organic  carbon  (TOC),  total 
phosphorus  (TP),  total  kjeldahl  nitrogen  (TKN),  aluminum  (Al), 
copper  (Cu),  ana  iron  (Fe)  with  conductivity  for  water  bodies  of  the 
northern  region.  Equations  with  the  best  fit  of  the  data  were  of 
the  following  format: 


1 inear 

Y = A + Bx 

(1) 

log/1  inear 

X 

CO 

cu 

<c 

II 

>- 

(2) 

log/log 

Y = AxB 

(3) 

where:  A = intercept,  B = slope,  Y = dependent  water  quality 

variable;  and 

X = independent  variable  (conductivity). 


Equation 

Format 

Water  Body  (1,  2,  3) 

df 

(n-2) 

Y 

A 

B 

p 

R2  (HckB^O) 

Jackfish 


Creek 

1 

10 

TOC 

-0.43 

0.046 

.69 

L0.001 

(DD0090) 

3 

10 

TP 

2.16 

0.583 

.17 

ns 

1 

10 

TKN 

362.60  - 

2.861 

.11 

ns 

3 

8 

Al 

1.00x10"° 

2.480 

.54 

L0.025 

1 

10 

Cu 

4.15 

-0.009 

.06 

ns 

3 

10 

Fe 

450.79 

0.168 

.01 

ns 

Big  Point 

Channel 

3 

41 

TOC 

527.95/- 

-0.709 

.23 

L0.01 

(DD0212) 

3 

44 

TP 

1.12x10° 

-1.742 

.35 

L0.001 

1 

30 

TKN 

589.48^ 

0.163 

.00 

ns 

2 

41 

Al 

1.97x10° 

-0.008 

.25 

L0.001 

1 

44 

Cu 

-2.4427 

0.034 

.02 

ns 

3 

44 

Fe 

6.50x10' 

-1.999 

.23 

L0.001 

Input  from  W: 

Mamawi  Lake 

1 

10 

TOC 

7.60 

0.012 

.03 

' ns 

Channel 

3 

11 

TP 

11.88 

0.268 

.01 

ns 

(KF0100) 

1 

11 

TKN 

838.22 

0.644 

.01 

ns 

1 

8 

Al 

2433.72 

-3.363 

.10 

ns 

3 

11 

Cu 

2.73 

0.197 

.00 

ns 

2 

11 

Fe 

383.37 

0.002 

.08 

ns 

Input  from  W:  ~ 

Prairie  Creek 

3 

10 

TOC 

8.66x10  " 

1.268 

.56 

L0.01 

(KF0140) 

1 

10 

TP 

181.08. 

-0.164 

.15 

ns 

3 

10 

TKN 

1.45x10* 

2.75x10^ 

-0.368 

.02 

ns 

3 

8 

Al 

-3.485 

.29 

ns 

2 

10 

Cu 

5.49, 

0.007 

.09 

ns 

2 

10 

Fe 

1.03x10* 

-0.003 

.27 

ns 

continued  . 


Table  6.  Concluded 
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Equation 

Format 

Water  Body  (1,  2,  3) 

df 

(n-2) 

Y 

A 

B 

p 

R2  ( Ho : B 1=0 ) 

Output: 

Riv  des 

3 

43 

TOC 

1.20 

0.339 

.06 

ns 

Rochers 

1 

46 

TP 

”3.01 

0.395 

.08 

ns 

(NA0030) 

3 

45 

TKN 

136.05 

0.280 

.01 

ns 

2 

45 

A1 

524.79 

-0.005 

.02 

ns 

2 

48 

Cu 

2.28 

0.004 

.04 

ns 

2 

48 

Fe 

511.83 

0.002 

.00 

ns 

Output: 

Chenel  des 

2 

25 

TOC 

8.25 

0.0006 

.01 

ns 

Quatres 

1 

31 

TP 

207.05 

-0.294 

.04 

ns 

Fourches 

1 

29 

TKN 

-1118.26 

6.650 

.14 

L0.05 

(KFU101) 

1 

29 

A1 

1357.63 

-2.034 

.03 

ns 

3 

31 

Cu 

250.14 

-0.681 

.06 

ns 

1 

31 

Fe 

3841.05 

-6.761 

.06 

ns 

L = less  than 
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relationships,  however,  were  established  for  the  other  five  channels 
which  drain  lakes.  The  association  between  total  phosphorus  and 
conductivity  in  rivers  reflects  the  direct  relationship  between  apatite 
phosphorus  and  discharge. 

There  was  a significant  positive  relationship  between  TKN  and 
conductivity  in  Chenal  des  Quatres  Fourches,  the  output  channel  draining 
the  northwest  end  of  Lake  Athabasca.  There  were  no  significant 
relationships  between  the  two  parameters  in  the  other  five  channels. 

Significant  relationships  between  aluminum  and  conductivity 
were  noted  in  the  channels  which  enter  Lake  Athabasca  from  the  southern 
Athabasca  River  drainage,  but  not  from  the  lentic  inputs  from  the  west 
or  the  two  outlet  channels.  There  was  an  inverse  relationship  between 
aluminum  and  conductivity  in  Big  Point  Channel,  but  a direct 
relationship  between  parameters  in  Jackfish  Creek,  the  channel  which 
drains  Richardson  Lake. 

No  significant  relationships  were  established  between  copper 
and  conductivity  in  any  of  the  drainages. 

There  was  a significant  negative  relationship  between  iron 
concentrations  arid  conductivity  in  Big  Point  Channel  but  not  for  any  of 
the  lake-draining  channels.  The  negative  relationship  between  iron  and 
conductivity  implies  a direct  association  between  iron  and  discharge. 
Thus,  the  increased  iron  concentrations  in  surface  waters  during  high 
flows  resulted  from  weathering  of  surficial  materials. 
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5.  THE  EAST 

This  region  includes  watersheds  east  of  the  Athabasca  mainstem 
between  the  Clearwater  River  in  the  south  and  the  delta  in  the  north,  a 
range  in  latitude  of  about  1°45‘.  The  four  major  watersheds  of  the 
eastern  region  are  the  Firebag,  Muskeg,  Steepbank,  and  Richardson.  The 
Richardson  drains  north  into  the  delta;  the  others  drain  west  into  the 
Athabasca  River. 

5.1  BASIN  FEATURES 

5.1.1  Land  Formations 

Land  formations  of  the  eastern  region,  which  have  been 
delineated  by  Turchenek  and  Lindsay  (1982),  are  presented  in  Figure  18. 

The  eastern  region  lies  entirely  within  the  Boreal  Mixedwood 
Forest.  This  forest  type  is  characterized  by  aspen  and  white  spruce, 
with  jack  pine  present  on  well-drained  sandy  soils  and  black  spruce  on 
wetland  communities. 

The  bedrock  geology  of  the  northern  portion  of  the  region  is 
characterized  by  Devonian  rocks  (limestone  and  shale)  of  the  Paleozoic 
era  and  in  the  extreme  northeast  by  granite  rocks  of  the  Archean. 
Younger  Cretaceous  rocks  of  sandstone  and  shale  occur  in  the  southern 
half  of  the  region. 

The  Steepbank  River  originates  in  the  Muskeg  River  uplands  and 
flows  through  level  and  undulating  plains  before  entering  the  Athabasca 
River  in  the  McMurray  lowland.  Although  the  Muskeg  and  Richardson 
rivers  originate  in  upland  areas,  the  major  portion  of  each  watershed 
flows  through  undulating  plains.  The  Firebag  watershed  lies  within  the 
Great  Slave  Plain. 

5.1.2  Watershed 

The  following  is  a brief  description  of  the  four  major 
drainages  within  the  eastern  region  (Figure  19). 

The  Steepbank  River  watershed  (1425  km2)  encompasses  the  land 
drained  by  the  mainstem  and  the  North  Steepbank  River.  The  Steepbank 
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Figure  18.  Lond  formotions  of  the  eastern  region. 

( Adopted  from  Turehenek  ond  Lindeoy  , I9B2  ) 
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Figure  19. 


East  surface  water  drainage  . 

( Adopted  from  Nolll  end  Evan*,  1979  ) 
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River  flows  120  km  northwest  to  its  confluence  with  the  Athabasca  River. 
The  North  Steepbank  River  drains  an  area  of  525  km2  before  joining  the 
Steepbank  River  about  32  km  upstream  from  the  Athabasca  mainstem. 

The  upper  reaches  of  the  watershed  are  meandering,  with 
extensive  beaver  activity  impounding  90  to  95%  of  the  channel  (Griffiths 
1973).  River  gradient  increases  downstream  from  the  junction  of  the 
North  Steepbank  river  and  mainstem,  forming  alternating  sections  of 
riffles  and  pools.  Outcroppings  of  the  bitumen-impregnated  McMurray 
Formation  and  the  Devonian  Waterways  Formation  occur  in  the  downstream 
15-km  stretch  of  the  Steepbank  River. 

The  Muskeg  River  watershed  (1480  km2)  drains  an  area  of  forest 
and  muskeg.  The  Muskeg  River  originates  in  the  Muskeg  Mountains  uplands 
and  flows  110  km  to  its  confluence  with  the  Athabasca  River.  Hartley 
Creek  (drainage  area,  325  km2)  joins  the  Muskeg  River  33  km  upstream  of 
its  confluence.  Stanley  Creek  (66  km2),  a minor  tributary,  flows 
southeasterly  to  the  Muskeg  River.  Kearl  Lake  (surface  area,  5.4  km2; 
maximum  depth,  2 m)  occurs  within  the  watershed.  The  Kearl  Lake 
tributary  enters  the  Muskeg  River  from  the  southeast  between  Stanley  and 
Hartley  Creeks. 

The  upper  reaches  of  the  Muskeg  River  are  well  drained  as  the 
river  flows  over  steep  channel  slopes.  As  the  river  flows  to  the 
southwest,  gradient  is  reduced  and  muskeg  is  a characteristic  feature  of 
the  watershed.  The  river  gradient  increases  in  the  downstream  16  km 
section  where  oil  sands  of  the  McMurray  Formation  and  limestone  of  the 
Devonian  Waterways  Formation  are  exposed  (Northwest  Hydraulics 
Consulting  Ltd.  1975).  Bitumen-impregnated  sands  also  are  exposed  in 
the  lower  reaches  of  Hartley  Creek. 

The  Fi rebag  River  watershed  (6138  km2)  drains  into  the 
Athabasca  River,  113  km  north  of  Fort  McMurray.  The  Marguerite  River 
(1746  km2)  and  the  Firebag  River  above  the  confluence  of  the  Marguerite 
are  the  two  major  sub-basins  of  the  watershed.  The  confluence  of  the 
sub-basins  is  about  45  km  upstream  from  the  mouth  of  the  Firebag. 

Headwater  regions  for  the  sub-basins  occur  outside  of  the 
AOSERP  study  area.  Upstream,  watercourses  meander  through  the  Firebag 
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Plain;  beaver  dams  are  frequent.  Stream  gradient  is  moderate  and 
bitumen  deposits  are  exposed  along  the  Fi rebag  and  Marguerite  Rivers. 
Downstream  of  the  confluence  of  the  Marguerite,  river  gradient  is  low  as 
the  meandering  river  enters  the  Athabasca  Plain. 

Groundwater  contribution  to  the  Fi rebag  watershed  is 
substantial,  maintaining  base  flows  and  moderating  water  temperatures  in 
summer  (Northwest  Hydraulics  Consulting  Ltd.  1975). 

The  Richardson  River  watershed  is  located  on  granite  rocks  of 
the  Precambrian  Shield.  The  river  drains  a watershed  area  of  2925  km2 
of  which  56%  is  located  in  the  province  of  Saskatchewan  (Northwest 
Hydraulics  Consulting  Ltd.  1975).  The  Richardson  River  forms  a 
sand/silt  delta  where  it  enters  the  Athabasca  River,  about  35  km 
upstream  of  Lake  Athabasca.  Sand  riverbed  and  bank  material  in  the 
downstream  Richardson  River  reach  enhances  percolation  of  precipitation 
to  groundwater  sources.  Numerous  lakes  within  the  watershed  also  reduce 
surface  runoff  and  maintain  high  base  flows. 

5.1.3  Physical  Description  of  Study  Streams 

Erosional  and  depositional  habitats  can  be  discerned  from 
elevation  profiles  of  rivers  (Figure  20).  Sekerak  and  Walder  (1980) 
present  elevation  profiles  and  associated  channel  characteristics  for 
the  Fi rebag.  Muskeg,  and  Steepbank  rivers.  Although  no  information  is 
available  on  channel  characteristics  for  the  Richardson  River,  general 
features  may  be  inferred  from  the  elevation  profile  of  the  river  (NTS 
maps  74E,  F,  and  L). 

The  Richardson  River  was  arbitrarily  subdivided  into  reaches 
based  on  the  elevation  profile  presented  in  Figure  20.  The  major 
portion  of  the  river  has  characteristical ly  steep  gradients  (Reach  4, 

5 m/km;  Reach  3,  9.1  m/km;  Reach  2,  4.6  m/km).  Steep  gradients  are 
typical  of  erosional  habitats  which  are  characterized  by  fast  flowing 
waters  and  coarse  substrates.  The  downstream  reach  with  a gradient  of 
0.5  m/km  is  likely  characterized  by  placid  flows  and  fine  substrates. 

The  Firebag  River  is  subdivided  into  four  upstream,  steep 
gradient  reaches  and  two  downstream,  low  gradient  reaches.  Riffle  areas 
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DISTANCE  FROM  THE  ATHABASCA  RIVER  (km) 


Figure  20.  Elevation  profiles  of  rivers  in  the 
eastern  region  . 

( Adapttd  Stkerak  and  Waldar  , I960  ) 
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and  coarse  substrates  predominate  in  the  upstream  reaches;  pools  and 
fine  substrates  predominate  in  the  downstream  reaches.  High  banks  are 
prevalent  along  the  entire  length  of  the  river,  particularly  in  Reaches 
1 and  6. 

The  Muskeg  River  channel  is  characterized  by  steep  gradients 
in  the  upstream  and  downstream  reaches.  The  river,  however,  is  low  in 
gradient  for  most  of  its  length  (Reach  4).  Beaver  dams  and  associated 
pools  are  prevalent  within  Reach  4.  Mud  and  silt  characterize  the 
greater  portion  of  the  Muskeg  River,  shifting  to  coarse  gravel  and 
boulders  in  the  downstream  15  km. 

Sekerak  and  Walder  (1980)  restricted  their  investigation  of 
the  Steepbank  River  to  the  downstream  reaches.  Here,  erosional  habitats 
prevail.  Although  Reach  5 is  depositional  in  nature,  there  is  a shift 
from  fine  to  coarse  substrates  and  from  pools  to  riffles  in  the 
downstream  45  km  of  the  river.  Reach  2 is  characterized  by  canyon  walls 
60  km  high. 

5.2  BASELINE  DATA 

5.2.1  Discharge 

Five  Water  Survey  of  Canada  gauging  stations  within  the 
eastern  region  were  selected  for  study  (Table  7;  Figure  4).  The 
stations  were  located  near  the  mouth  of  the  major  rivers  (Steepbank, 
Muskeg,  Hartley,  Firebag,  and  Richardson)  and  the  recorded  flow 
represented  surface  water  drainage  from  each  of  the  watersheds. 

Mean  daily  discharge  (1  in  5 days)  is  presented  for  the 
Steepbank  (1976  to  1980),  Muskeg  (1976  to  1981),  Hartley  Creek  (1976  to 
1981),  Firebag  (1976  to  1979)  and  Richardson  (1977  to  1979)  rivers  in 
Figures  21  to  25.  The  discharge  years  were  selected  for  illustration 
because  they  corresponded  to  the  water  quality  sampling  regime. 

Of  the  five  rivers,  the  Firebag,  with  the  largest  drainage 
area  (6138  km2),  contributed  the  greatest  flow  to  the  Athabasca 
mainstem.  The  mean  annual  discharge  for  the  Firebag  River  was 
22.0  m3/s , 30.0  m3/s,  and  27.4  m3/s  in  1977,  1978,  and  1979, 


Table  7.  Water  Survey  of  Canada  gauging  stations  in  the  east. 
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Figure  21.  Daily  discharge  (1  in  5)  in  the  Steepbank  River,  1976-1980. 
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Figure  23.  Daily  discharge  (I  in  5)  in  Hartley  Creek,  1976-1981. 
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Figure  24.  Daily  discharge  (I  in  5)  in  the  Firebag  River,  1976-1979. 
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Figure  25.  Daily  discharge  (I  in  5)  in  the  Richardson 
River,  1977-1979. 
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respectively.  Mean  flow  data  were  unavailable  in  1976  because  of 
interrupted  records.  Maximum  flows  occurred  in  April  and  May,  with 
maximum  mean  monthly  flows  ranging  from  37.7  to  113.5  m3/s  between  1976 
and  1979.  Elevated  flows  occurred  occasionally  in  summer  or  fall  as  a 
result  of  storms.  Winter  base  flows  ranged  from  7.5  to  11.0  m3/s 

The  mean  annual  discharge  for  the  Richardson  River  was 

15.2  m3/s , 17.0  m3/s , and  16.1  m3/s  in  1977,  1978,  and  1979, 
respectively.  Maximum  mean  monthly  discharge  occurred  in  May,  with 
flows  ranging  from  18.4  to  28.9  m3/s  during  the  three  years.  Elevated 
flows  may  occur  in  the  fall  as  a result  of  flood  conditions  (e.g.,  1976, 
1978).  Winter  base  flows  ranged  from  9 to  15  m3/s. 

The  mean  annual  discharge  for  the  Muskeg  River  varied 

considerably , ranging  from  1.3  to  5.9  m3/s  during  a six-year  period 
(1976  to  1981).  Maximum  mean  monthly  flows  occurred  during  September  in 
1978  (21.9  m3/s)  and  1980  (15.5  m3/s).  Maximum  mean  monthly  flows 
during  spring  runoff  (April/May)  ranged  from  2.7  to  21.4  m3/s  during  the 
six  years.  Winter  base  flows  ranged  from  0.2  to  0.5  m3/s  (Figure  22). 

The  discharge  pattern  for  Hartley  Creek,  a major  tributary  of 
the  Muskeg  River,  followed  that  of  the  larger  river.  Discharge 

differences  occurred  only  in  amplitude  (Figure  23). 

Mean  annual  flow  in  the  Steepbank  River  ranged  from  3.0  to 
6.4  m3/s  during  the  five-year  period  (1976  to  1980).  Maximum  mean 
monthly  flow  occurred  in  September  (1978  and  1980)  as  a result  of 
rainstorms.  Maximum  mean  monthly  flows  during  spring  runoff  ranged  from 
3.6  to  30.1  m3/s  during  the  five  years.  Winter  base  flows  ranged  from 
0.2  to  0.5  m3/s. 

5.2.2  Water  Quality 

Water  quality  sampling  sites  within  the  eastern  region  are 
presented  in  Table  8.  Only  sites  with  sufficient  data  (more  than  10 
samples)  were  selected  for  analysis.  Data  used  in  the  analysis  are 
summarized  in  Appendix  10.4  (Tables  62  to  64).  Water  quality  data  were 
available  for  each  of  the  major  rivers  within  the  region  including 
Steepbank,  Muskeg,  Hartley,  Firebag,  and  Richardson  rivers. 
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Table  8.  Description,  period  of  record  and  number  of  samples  obtained  for 
water  quality  sampling  stations  within  the  eastern  region.  See 
Figure  4 for  location  of  sites. 


Latitude  N 

Code  Longitude  W Period  of  Number 

Watershed  (00AT07-)  (D  M S)  Description  Record  of  Samples 


Muskeg 


DA0079 

57 

08 

05 

Muskeg  River  near 

1976-78 

10 

111 

36 

08 

mouth 

DA0080 

57 

11 

30 

Muskeg  River  at 

1976-81 

77a 

111 

34 

05 

at  WSC  gauge 

DA0081 

57 

08 

03 

Muskeg  River  above 

1976-81 

17a 

111 

35 

53 

Shell  pumping  above 
Hartley  Creek 

DA0082 

57 

14 

18 

Hartley  Creek 

1976-81 

10 

111 

24 

55 

3.2  kmupstream  from  mouth 

DA0083 

57 

21 

08 

Stanley  Creek 

1976-77 

lla 

111 

22 

44 

DA0084 

57 

18 

14 

Kearl  Lake  tributary 

1976-78 

16a 

111 

22 

30 

1.6  km  upstream  from 
Muskeg  River 

DA0085 

57 

25 

00 

Muskeg  River  11.2  km 

1976-81 

26a 

111 

13 

16 

upstream  from  Stanley 
Creek 

DA0086 

57 

22 

48 

Muskeg  tributary 

1976-77 

10 

111 

16 

44 

5.6  km  upstream  from 
Stanley  Creek  and 
0.8  km  upstream  from 
Muskeg  River. 

- 

DA0087 

57 

11 

21 

Hartley  Creek  SW  fork 

1976-77 

9 

111 

23 

44 

0.4  km  from  junction 
with  SE  fork 

DA0088 

57 

09 

23 

Hartley  Creek,  SE 

1976-77 

9 

111 

23 

27 

fork  0.4  km  from 
junction  with  SW  fork 

DA0089 

57 

15 

08 

Tributary  to  Muskeg 

1977 

7 

111 

21 

54 

River,  4.8  km  upstream 

from  Hartley  Creek 


continued  . 


Table  8.  Concluded. 
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Latitude  N. 


Code 

Longitude  W. 

Period  of 

Number 

Watershed 

(00AT07-) 

(D 

M 

S) 

Description 

Record 

of  Samples 

DA0090 

57 

15 

34 

Hartley  Creek, 

1976-81 

53d 

111 

27 

53 

0.4  km  above 
confluence  with 
Muskeg  River,  WSC 
gauge 

DA0091 

57 

15 

42 

Tributary  leading 

1976-77 

10 

111 

19 

18 

to  Kearl  Lake 

DA0092 

57 

16 

15 

Kearl  Lake  outlet 

1976-78 

10 

111 

15 

01 

DA0093 

57 

16 

28 

Kearl  Lake  inlet 

1976-77 

7 

111 

13 

30 

DA0094 

57 

20 

41 

Muskeg  River, 

1976-77 

6 

111 

07 

50 

22.5  km  upstream 
from  Stanley  Creek 

DA0095 

57 

06 

06 

Hartley  Creek,  SW 

1976-77 

10 

111 

23 

07 

fork  16  km  from 
junction  with  SE 
fork 

DA0096 

57 

04 

22 

Hartley  Creek,  SE 

1976-77 

8 

111 

11 

18 

fork  21  km  from 
junction  with  SW 
fork 

Steepbank 

DA0060 

57 

00 

17 

Steepbank  River, 

1976-80 

39a 

111 

24 

53 

7.2  km  upstream  from 

mouth,  WSC  gauge 

Fi rebag 

DC0010 

57 

38 

30 

Firebag  River,  WSC 

1976-79 

31a 

111 

10 

30 

gauge 

DC0011 

57 

44 

35 

Firebag  River  near 

1977 

3 

111 

21 

03 

mouth 

DC0020 

57 

17 

20 

Lost  Creek,  0.8  km 

1976-78 

9 

110 

27 

50 

above  mouth,  WSC 
gauge 

Richardson 

DD0020 

58 

21 

48 

Richardson  River 

1977-79 

15a 

111 

14 

14 

WSC  gauge 

a 

Sites  with  sufficient  numbers  to  be  analysed. 
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5.2.2. 1 Routine  parameters.  A summary  of  routine  parameters  for  the 
region  is  presented  in  Appendix  10.4  (Table  62). 

Surface  waters  of  the  Richardson  River  are  relatively  low  in 
pH  levels  (6.7  to  7.8),  poor  in  buffering  capacity  (alkalinity,  maximum 
47  mg/L)  and  low  in  conductivity  (maximum  100  uS/cm).  These  water 
quality  values,  which  are  characteristic  of  surface  waters  flowing  over 
granite  rock  of  the  Precambrian  Shield,  clearly  differed  from  waters  of 
other  drainages  (Steepbank,  Firebag,  and  Muskeg)  which  lie  within 
sedimentary  rock  of  the  Interior  Plain.  Surface  waters  of  the 
Steepbank,  Muskeg,  and  Firebag  watersheds  were  relatively  high  in  pH 
values  (max.  range  8.1  to  8.7),  well  buffered  (alkalinity  maximum  184  to 
362  mg/L)  with  elevated  conductivity  levels  (max.  range  320  to  704 
uS/cm).  Values  of  filterable  residue  (total  dissolved  solids)  varied 
directly  with  conductivity. 

The  relative  contribution  of  major  ions  is  presented  in 
Figure  26.  The  extent  of  the  ion  concentration  reflected  in  the  Stiff 
diagrams  illustrates  the  low  conductivity  and  restricted  range  of  ion 
values  in  the  Richardson  River  compared  to  levels  in  the  Steepbank, 
Muskeg,  and  Firebag  watersheds.  Calcium  bicarbonate  was  the  dominant 
dissolved  salt  in  all  watersheds. 

Other  differences  between  granite  and  sedimentary  rock  include 
the  thickness  of  overlying  top  soil  and  weathering  of  rock  material. 
Soil  layers  in  the  Precambrian  Shield  are  thin  and  the  rock  material 
erodes  slowly.  Soil  layers  overlying  sedimentary  rock  are  thick  and 
high  in  organic  content;  slumping  of  river  banks  contributes  to  the 
sediment  load.  Thus,  median  values  of  total  organic  carbon  (3  mg/L)  and 
total  inorganic  carbon  (10  mg/L)  observed  in  the  Richardson  River  were 
noticeably  lower  than  corresponding  values  for  the  other  watersheds. 
The  muskeg  drainage  of  the  Firebag,  Muskeg,  and  Steepbank  contributes  to 
the  organic  carbon  in  the  three  watersheds. 

Non-fil terable  residue  (suspended  solids)  tends  to  vary 
directly  with  discharge  levels.  Thus,  maximum  residue  values  are 
recorded  during  spring  runoff  or  other  storm  events.  Again,  values  for 
the  Richardson  River  were  consistently  lower  than  the  other  watersheds 
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Figure  26.  Maximum  and  minimum  values  (me/L)  of 
major  ions  in  rivers  of  the  east. 
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(Appendix  10.4,  Table  62).  The  maximum  value  of  non-fil terable  residue 
was  recorded  in  Hartley  Creek  (459  mg/L). 

Colour  is  a measure  of  water  clarity.  Maximum  values  were 
lower  in  the  Richardson  (50  units)  and  Fi rebag  (90  units)  than  in  the 
Steepbank  (180  units),  Muskeg  (140  units),  or  Hartley  Creek  (140  units). 

Sources  of  phenolic  material  include  industrial  and  municipal 
wastes  as  well  as  decaying  plant  material.  Values  for  the  Richardson 
River  were  below  the  detection  limit.  Maximum  recorded  values  for  the 
other  sampling  sites  exceeded  the  Alberta  surface  water  objectives 
(0.0005  mg/L).  Similarly,  oil  and  grease  values  for  the  Richardson 
River  (max.  0.7  mg/L)  were  substantially  lower  than  for  the  other  rivers 
(max.  range  2.1  to  12.9  mg/L).  Highest  oil  and  grease  levels  were 
recorded  for  the  Muskeg  River. 

5. 2. 2.2  - Nu-trients.  Phosphorus  and  nitrogen  are  the  major  nutrients 
controlling  plant  growth.  Accordingly,  the  constituents  in  surface 
water  vary  seasonally. 

Total  phosphorus  and  orthophosphate  phosphorus  (PO^-P)  are  the 
two  forms  of  phosphorus  recorded  in  Appendix  10.4  (Table  63). 

Although  total  phosphorus  increases  with  increasing  discharge,  it  is  the 
particulate  not  the  dissolved  portion  which  increases  substantially  as  a 
result  of  runoff.  PO^-P  increases  as  a result  of  enrichment  (natural  or 
domestic  sewage).  Total  phosphorus  values  of  the  Richardson  River  (max. 
67  ug/L)  were  less  than  the  other  rivers  within  the  eastern  region  (max. 
range  138  to  330  ug/L),  reflecting  the  low  levels  of  erosion  from 
granite  material.  The  P0^-P  values  for  the  Richardson  River  (max.  42 
ug/L)  were  within  the  range  recorded  for  the  other  sampling  sites  (max. 
range  28  to  150  ug/L).  Highest  values  were  recorded  in  the  Steepbank 
Rivers. 

Nitrate/nitrite,  ammonia,  and  total  kjeldahl  nitrogen  are  the 
three  forms  of  nitrogen  regularly  recorded  in  surface  waters.  Total 
kjeldahl  nitrogen  is  a measure  of  total  organic  nitrogen  plus  ammonia. 
Inorganic  nitrogen  includes  ammonia  and  nitrate/nitrite  nitrogen. 
Nitrogen  forms  vary  seasonally  as  a function  of  dilution  at  high 
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discharge  levels,  temperature,  uptake  by  plants,  and  loss  to  the 
atmosphere. 

Richardson  River  ammonia  values  (max.  0.1  mg/L)  were  lower 
than  the  other  sampling  sites  (0.2  to  3.9  mg/L).  Values  for  other 
nitrogen  forms  in  the  Richardson  River  were  within  the  range  recorded  at 
other  sampling  sites  (Appendix  10.4,  Table  63). 

5. 2. 2. 3 Metals.  Concentrations  of  metals  in  surface  waters  are 
attributed  to  weathering  of  surficial  material  and  groundwater 
discharge,  as  well  as  municipal  and  industrial  wastes.  The  metals, 
which  were  predominant  in  surface  waters  of  the  eastern  region,  were 
iron,  manganese,  aluminum,  and  arsenic  (Appendix  10.4,  Table  64). 
Concentrations  of  boron,  cadmium,  chromium,  lead,  selenium,  and  silver 
in  the  rivers  were  all  less  than  the  Alberta  Surface  Water  Objectives. 

Median  values  of  iron  concentrations  exceeded  the  Alberta 
objectives  (300  ug/L  or  0.3  mg/L)  for  all  five  eastern  tributaries.  The 
range  of  iron  concentrations  for  the  rivers  are  as  follows: 


Code 

River 

Range  Fe  (mg/ 

DA0060 

Steepbank 

0.42  to  3.70 

DA0080 

Muskeg 

0.35  to  10.50 

DA0090 

Hartley 

0.40  to  10.70 

DC0010 

Firebag 

0.07  to  2.06 

DD0020 

Richardson 

0.36  to  0.82 

Iron  concentrations  typically  vary  directly  with  discharge, 
and  suspended  material  is  a result  of  weathering  processes.  The  maximum 
iron  concentration  of  the  Richardson  River  was  lower  than  the  other 
tributaries  and  was  related  presumably  to  the  lower  erosion  levels 
characteristic  of  granite  rock. 

Manganese  is  similar  to  iron  in  its  chemical  behaviour  and  is 
often  found  in  association  with  iron  (McNeely  et  al.  1979).  Maximum 
manganese  concentrations  in  surface  waters  of  the  Steepbank  (500  ug/L), 
Muskeg  (2300  ug/L),  Hartley  (2150  ug/L),  and  the  Firebag  (790  ug/L) 
rivers  exceeded  the  Alberta  Surface  Water  Objectives  of  50  ug/L. 


88 


Manganese  concentrations  in  the  Richardson  River  were  within  the 
provincial  guidelines. 

Aluminum  concentrations  in  surface  waters  are  limited  to 
soluble  salts.  The  maximum  values  of  the  eastern  tributaries  (except 
the  Steepbank  River)  were  less  than  1000  ug/L  (Appendix  10.4,  Table  64) 
which  is  within  the  range  of  aluminum  concentrations  in  most  surface 
waters  (McNeely  et  al . 1979).  The  maximum  aluminum  value  (5600  ug/L) 
was  recorded  in  the  Steepbank  River  during  runoff  (1978  May  03). 

Non-compliance  with  surface  water  objectives  was  noted  for 
arsenic,  copper,  and  mercury.  Arsenic  concentrations  exceeded  the  ASW0 
(0.01  mg/L)  in  the  Steepbank,  Muskeg,  Hartley,  and  Firebag  Rivers.  Most 
elevated  values  were  recorded  in  winter  months.  Mercury  concentrations 
exceeded  the  objectives  (0.1  ug/L)  in  all  five  major  tributaries. 

5.2.3  Benthic  Invertebrates 

A summary  of  benthic  invertebrate  studies  conducted  in  streams 
throughout  the  eastern  region  is  presented  in  Table  9.  Because  several 
projects  were  initiated  by  proponents  of  oil  sands  development  including 
Alsands  (Mayhood  et  al.  1981;  Corkum  and  McCart  1981),  Petro  Canada 
(Corkum  and  McCart  1981)  and  Sandalta  (O'Neill  et  al.  1982),  sampling 
sites  were  frequently  restricted  to  lease  areas.  Other  invertebrate 
studies  were  conducted  to  provide  baseline  data  prior  to  oil  sands 
development  or  as  part  of  larger  projects  on  aquatic  resources  (Barton 
and  Wallace  1980;  Walder  et  al.  1980). 

Clearly,  there  was  no  consistant  approach  to  the  study  of 
invertebrates  with  regard  to  type  of  sampler,  sampling  effort,  mesh  size 
through  which  samples  were  washed,  or  sample  frequency.  Also,  no  one 
sampler  can  be  used  with  equal  efficiency  in  erosional  and  depositional 
habitats.  For  example,  kick  nets,  cylinder  and  surber  samplers  were 
used  in  erosional  habitats;  Ekman  grab  and  airlift  pumps  were  used  in 
depositional  habitats.  Sampler  selection  is  discussed  in  detail  by 
Merritt  et  al . (1978)  and  Mayhood  et  al . (1981). 


Table  9.  Summary  of  benthic  invertebrate  studies  in  streams  of  the  east. 
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Mayhood  et  al.  1980  Mul tipi  ate  8 07,  08,  10e  E/D  Ephemeroptera , Hydra 


Sampler 

Water  Body  Reference  Sample  Year  Mesh  Size  Sites  Month  Habitat  Dominant  Taxa 
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et  al.  (1979)  250  urn  Ephemeroptera 


Sampler 

Water  Body  Reference  Sample  Year  Mesh  Size  Sites  Month  Habitat  Dominant  Taxa 
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Relative  abundance  not  indicated,  but  list  of  taxa  provided  in  original  report. 

Multiplate  (Hester-Dendy)  artificial  sampler. 

Multi  pi  ate  samplers  positioned  in  river  during  May/June  and  retrieved  in  July,  August,  and  October; 
many  samples  in  August  and  October  were  lost. 

= Mesh  size  not  reported. 
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Differences  in  faunal  composition  recorded  in  eastern 
tributaries  reflect  differences  between  erosional  and  depositional 
habitats  (Table  9).  Erosional  habitats  are  characterized  by  shallow, 
swift  flowing,  well  oxygenated  water  with  a streambed  of  coarse 
substrates.  Depositional  habitats  are  characterized  by  deep,  slow 
flowing  water  often  low  in  oxygen  and  with  a streambed  of  fine  (organic, 
sand,  silt)  substrates.  Clean  water  organisms  such  as  Ephemeroptera , 
Plecoptera,  Trichoptera,  and  Orthocladiinae  occur  in  erosional  habitats. 
Oligochaeta,  Chironominae , and  Mollusca  are  typically  found  in 
depositional  habitats. 

Because  of  differences  in  sampling  techniques  and  level  of 
taxonomic  identification,  and  lack  of  access  to  the  raw  data,  a summary 
of  analyses  of  invertebrates  or  factors  influencing  their  distribution 
will  not  be  presented. 

5.3  RELATIONSHIPS  AMONG  PARAMETERS 

The  objective  was  to  develop  relationships  between  water 
quality  constituents  and  discharge.  Since  specific  conductance  is  a 
function  of  discharge,  a correlation  analysis  was  performed  between 
conductivity  and  other  water  quality  parameters  to  select  those  few 
significant  parameters  from  the  suite  available  for  subsequent  analysis 
with  reference  to  discharge. 

Results  from  the  correlation  analysis  performed  using  median 
values  of  water  quality  parameters  (listed  in  Appendix  10.4,  Tables  62 
to  64)  indicated  that  ammonia,  iron,  and  manganese  were  significantly 
correlated  with  conductivity  (P<0.01).  These  four  parameters  were  used 
to  develop  regression  relationships  with  discharge,  using  data  from  each 
sampling  date.  Relationships  were  developed  for  the  Richardson 
(DD0020) , Firebag  (DC0010),  Muskeg  (DA0080),  and  Steepbank  (DA0060) 
rivers  using  linear,  log/linear,  and  log/log  transformations  of  the 
data. 

Equations  with  the  best  fit  (i.e.,  least  variability  and 
highest  R2)  were  selected  and  are  presented  in  Table  10.  The 
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Table  10.  Summary  of  regression  equations  of  conductivity  (SC),  ammonia  (NHL), 
iron  (Fe),  and  manganese  (Mn)  with  discharge  for  the  major  eastern 
rivers.  Equations  with  the  best  fit  of  the  data  were  of  the 
following  format: 


1 inear 

Y = A + Bx 

(1) 

log/linear 

Y = AeBx 

(2) 

log/log 

Y = AxB 

(3) 

where:  A = intercept,  B = slope,  Y = dependent  variable  (water 
quality);  and 

X = independent  variable  (discharge),  df  = degrees 
of  freedom. 


Equation 

Format  df  p 

Water  Body  (1,  2,  3)  (n-2)  Y A B R2  (HckB^O) 


Richardson 


(DD0020) 

3 

13 

SC 

1 

11 

NH, 

2 

12 

Fe' 

1 

12 

Mn 

Firebag 

3 

31 

SC 

(DC0010) 

1 

24 

NH 

1 

26 

Fe' 

1 

26 

Mn 

Muskeg 

3 

66 

SC 

(DA0080) 

3 

66 

NH 

3 

66 

Fe' 

3 

66 

Mn 

Steepbank 

(DA0060) 

3 

37 

SC 

2 

33 

NH 

1 

37 

Fe 

2 

37 

Mn 

216.80 

-0.334 

.47 

L0.01 

88.38 

-4.067 

.33 

L0.05 

481.06 

0.011 

.05 

ns 

53.85 

-1.293 

.27 

ns 

531.13 

-0.332 

.69 

L0.001 

73.70 

-0.855 

.11 

ns 

420.44 

10.560 

.35 

L0.001 

29.77 

0.398 

.16 

L0.05 

498.70 

-0.397 

.69 

L0.001 

205.00 

-0.835 

.29 

L0.001 

1822.56 

-0.445 

.36 

L0.001 

260.86 

-1.012 

.44 

L0.001 

637.78 

-0.649 

.89 

L0.001 

61.81 

-0.101 

. .21 

L0.01 

835.07 

54.381 

.21 

L0.01 

28.25 

0.091 

.40 

L0.001 

L = less  than. 
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significance  of  each  regression  equation  was  tested  by  determining  if 
the  slope  differed  significantly  from  zero. 

There  was  a significant  inverse  relationship  between 
conductivity  and  discharge  for  four  rivers.  Thus,  conductivity  levels 
were  high  during  low  flow  periods  and  vice  versa.  The  four  regression 
analyses  were  compared  using  analysis  of  covariance.  However,  since 
residual  variances  among  rivers  were  heterogeneous  (Bartlett's  test  for 
more  than  two  groups),  further  comparisons  of  slope  and  intercept  were 
inappropriate. 

Ammonia  concentrations  were  also  inversely  correlated  with 
discharge,  with  highest  concentrations  of  ammonia  being  present  in 
winter  months.  The  regression  equation  between  variables  was  not 
significant  for  the  Firebag  River,  but  significant  associations  were 
established  for  the  Richardson,  Muskeg,  and  Steepbank  rivers.  Because 
the  residual  variances  among  the  Richardson,  Muskeg,  and  Steepbank 
rivers  were  heterogeneous,  a statistical  comparison  among  the  regression 
lines  was  not  conducted.  There  was,  however,  a trend  for  ammonia  levels 
to  decrease  more  rapidly  with  discharge  in  the  Richardson  River  than  in 
the  Muskeg  or  Steepbank  rivers. 

There  were  no  significant  predictable  relationships  between 
either  iron  and  discharge  or  manganese  and  discharge  for  the  Richardson 
River.  Iron  and  manganese  typically  increase  in  surface  waters  as  a 
result  of  weathering  of  the  surficial  material.  Because  soil  layers 
overlying  granite  rocks  are  usually  thinner  than  soils  overlying  ' 
sedimentary  rocks,  no  significant  association  between  parameters  was 
anticipated  in  the  Richardson  watershed. 

There  was  a significant  direct  relationship  between  iron  and 
manganese  with  discharge  in  the  Firebag  and  Steepbank  rivers. 
Surprisingly,  there  was  an  inverse  relationship  between  these  parameters 
in  the  Muskeg  River,  indicating  that  groundwater  inputs  may  have 
elevated  levels  of  iron  and  manganese. 
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5.4  IMPACT  OF  DEVELOPMENT 

Several  short-term  monitoring  programs  were  established 
throughout  the  Muskeg  basin  in  response  to  anticipated  oil  sands 
development.  Although  development  has  ceased  because  of  unfavourable 
economic  conditions,  one  company  (A1 sands  Energy  Limited)  had  cleared 
and  drained  land  from  both  mine  and  plant  sites  in  1980.  A study 
conducted  by  Mayhood  et  al . (1981)  was  designed  to  determine  the  effect 
of  mine  and  plant  drainage  on  the  Muskeg  River. 

Two  events  related  to  the  development  occurred  during  the 
environmental  study: 

1.  Flood.  Overflow  of  plant  and  mine  site  drainage  ditches 
(April  1980)  resulted  in  the  inundation  of  surrounding 
lands  and  subsequent  elevated  levels  of  suspended 
sediments  and  turbidity  in  the  Muskeg  River. 

2.  Eutrophication.  Fertilizer  application  to  the  mine  site 
drainage  area  resulted  in  temporary  nutrient  enrichment 
in  the  Muskeg  River. 

After  these  events,  changes  were  noted  in  the  benthic  algae 
and  invertebrates  in  the  Muskeg  River.  Changes  in  the  benthic 
invertebrate  community  which  occurred  during  the  June  sampling  period 
may  have  resulted  from  the  April  flood.  In  June,  total  benthic 
invertebrate  biomass  and  the  abundance  of  detritivores  and  predators 
were  reduced  in  the  Muskeg  River  about  100  m downstream  from  the 
drainage  site.  Although  algal  biomass  (cell  volume)  was  enhanced  as  a 
result  of  enrichment,  no  change  in  algal  community  composition  was 
noted.  These  biological  responses  were  short  term,  with  levels 
returning  to  pre-flood  conditions  on  subsequent  sampling  trips.  A 
settling  pond  and  stepped  drainage  ditch  were  constructed  by  the 
industry,  improving  the  quality  of  the  drainage  water  during  flood 
events. 
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6.  THE  WEST 

Watersheds  which  are  west  of  the  Athabasca  mainstem  include 
the  Beaver-Poplar,  MacKay,  Ells,  and  several  tributaries  draining  the 
east  slopes  of  the  Birch  Mountains.  Sampling  programs  were  not 
established  in  the  Mclvor  or  Buckton  watersheds,  both  of  which  are 
located  in  the  northern  portion  of  the  region. 

6.1  BASIN  FEATURES 

6.1.1  Land  Formations 

Land  formations  of  the  western  region  delineated  by  Turchenek 
and  Lindsay  (1982)  are  presented  in  Figure  27. 

The  Boreal  Mixedwood  and  Northern  Boreal  ecoregions  occur 
within  the  western  region.  The  Northern  Boreal  Ecoregion  occupies  the 
northwest  portion  of  the  study  area  which  includes  the  Mclvor  basin  and 
Namur,  Gardinar,  and  Eaglenest  lakes,  as  well  as  the  headwater  reaches 
of  several  streams.  The  Boreal  Mixedwood  Forest  occupies  the  remaining 
area  of  the  region. 

The  east  slope  streams  originate  in  the  Birch  Mountain  Uplands 
of  the  Northern  Boreal  Forest  where  there  are  extensive  areas  of  organic 
soils.  The  downstream  portions  of  the  rivers  lie  within  the  Boreal 
Mixedwood  Forest,  flowing  through  poorly  drained  soils  of  the  Buckton 
Plain.  The  Pierre,  Eymundson,  and  Asphalt  rivers  flow  through  the 
Embarras  Plain  to  the  Athabasca  River.  The  downstream  portion  of  th'e 
Calumet  lies  within  the  Muskeg  River  Plain. 

Peat  deposits  occupy  60  to  80%  of  the  Muskeg  River  Plain. 
Here,  the  surficial  deposits  overlie  the  oil  sands  of  the  McMurray 
Formation  or  limestone  of  the  Devonian  Waterways  Formation.  In 
contrast,  the  northern  Embarras  Plain  is  characterized  by  sandy  outwash 
deposits  with  a few  (ca.  10%)  peatland  areas  and,  in  some  places,  large 
migrating  dunes.  Underlying  Devonian  rock  results  in  the  formation  of 
sinkholes. 

The  Ells  and  the  MacKay  rivers  flow  through  the  same 
physiographic  regions  from  headwater  (Wabasca-Athabasca  Lowland)  and 
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27.  Land  formations  of  the  western  region  . 

( Adopt.d  from  Turchtntk  and  Llndtay,  1982  ) 
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middle  (Methy  Portage  Plain)  to  downstream  (Great  Slave  Plain)  reaches. 
The  Poplar-Beaver  Creek  drainage  lies  within  the  Methy  Portage  Plain 
(upstream  and  downstream  reaches)  and  the  Great  Slave  Plain  (middle 
reaches) . 

Peat  deposits  and  associated  wetland  vegetation  occur  in  80  to 
90%  of  the  headwater  drainage  area  of  the  Ells  and  MacKay  rivers. 
Peatlands  occupy  only  20  to  40%  of  the  Methy  Portage  Plain  (specifically 
the  Dover  Plain).  Soils  of  the  Dover  Plain  are  characteristically 
saline.  The  Muskeg  River  Plain  of  the  Great  Slave  Plain  is 
characterized  by  glaciofluvial  sands  overlying  oil  sands  and  limestone. 

6.1.2  Watersheds 

Four  major  watersheds  occur  along  the  west  bank  of  the 
Athabasca  River,  north  of  Fort  McMurray,  including  (Figure  28): 

1.  Beaver-Poplar  Creek; 

2.  MacKay  River  (including  Dunkirk  and  Dover  rivers); 

3.  Ells  River; 

4.  Several  rivers  draining  the  east  slope  of  the  Birch 
Mountains  (Tar  River,  Pierre  River,  Asphalt  Creek,  and 
Eymundson  Creek). 

The  Beaver-Poplar  Creek  Diversion  was  constructed  in  1977  to 
divert  water  from  the  mining  operations  of  Syncrude  Canada  Ltd.  A dam 
constructed  on  the  slow  flowing  Beaver  River  formed  the  Beaver 
Reservoir.  Diversion  canals  drain  water  from  the  reservoir  to  Ruth  Lake 
and  from  the  lake  into  Poplar  Creek  Reservoir.  Water  draining  the 
reservoir  enters  (via  a spillway)  Poplar  Creek  which  subsequently  flows 
into  the  Athabasca  River.  Four  streams  which  originally  entered  Beaver 
Creek  now  enter  a channel  which  is  called  the  West  Interceptor  Ditch. 
Three  of  the  streams  flow  north  of  the  plant  site  to  the  original  Beaver 
Creek.  The  fourth  stream  flows  southward  to  the  Beaver  Creek  Reservoir. 

The  MacKay  River  watershed  (5517  km2)  is  the  largest  basin  in 
the  west  bank  region  (Sekerak  and  Walder  1980).  The  MacKay  River  flows 
northeasterly  for  200  km  with  an  average  gradient  of  0.2%  and  drains  an 
area  of  2350  km2.  The  Dunkirk  (2183  km2)  and  Dover  (984  km2)  sub-basins 
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Figure  28.  West  surface  water  drainage 

( Adapted  from  Neill  and  Evone , 1979  ) 
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drain  muskeg  areas.  The  Dunkirk  flows  150  km  southeasterly  and  enters 
the  MacKay  River  89  km  upstream  from  its  confluence  with  the  Athabasca 
River.  The  Dover  flows  northeasterly  and  enters  the  MacKay  River  16  km 
upstream  from  its  mouth. 

The  Ells  River  watershed  (2707  km2)  has  a mean  gradient  of 

0.2%  (Northwest  Hydraulics  Consulting  Ltd.  1975;  Hickman  et  al . 1980). 
Joslyn  (264  km2)  and  Chelsea  (202  km2)  creeks  are  the  two  major 
sub-basins  within  the  watershed.  Lakes  in  the  Birch  Mountain  headwaters 
also  drain  (1191  km2)  into  the  Ells  River.  The  high  storage  capacity  of 
the  basin,  owing  to  large  lake  and  muskeg  areas,  moderates  surface 
runoff  and  maintains  high  base  flows  in  winter. 

East  draining  streams  of  the  Birch  Mountain  uplands  include: 

1.  Tar  River  (324  km2); 

2.  Calumet  River  (179  km2); 

3.  Pierre  River  (137  km2); 

4.  Unnamed  River  (337  km2);  and 

5.  Eymundson  Creek  (303  km2). 

Permafrost  conditions  probably  contribute  to  valley  wall 
instability  in  the  Birch  Mountain  region.  Characteristic  stream  bank 
slumping  results  in  high  sediment  loads  for  the  rivers. 

6.1.3  Physical  Description  of  the  Study  Streams 

Beaver  Creek  Reservoir  (max.  depth  8 m)  and  Poplar  Reservoir 
(max.  depth  16  m)  are  the  two  lentic  areas  formed  by  the  diversion 
scheme  of  Syncrude  Canada  Ltd.  Groundwater  discharge  from  the  mine 
enters  the  Beaver  Creek  Reservoir  and  ultimately  enters  the  Athabasca 
River  via  diversion  canals  through  Ruth  Lake,  Poplar  Creek  Reservoir, 
and  Poplar  Creek. 

Prior  to  diversion,  the  lower  section  of  Poplar  Creek  was 
characterized  as  a slow  moving  meandering  stream  with  sand/silt 
substrates.  This  section  of  Poplar  Creek  was  subsequently  channelized 
to  receive  increased  flows  (Noton  and  Chymko  1978). 

The  MacKay  River  is  characterized  by  depositional  and 
erosional  habitats.  Sekerak  and  Walder  (1980)  subdivided  the  MacKay 
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River  into  seven  reaches  including  five  relatively  long  (26  to  46  km) 
upstream  and  two  short  (1.4  to  3 km)  downstream  reaches  (Figure  29). 
The  two  depositional  upstream  reaches  of  the  MacKay  River  are 
characterized  by  a gradual  slope  (0.5  to  0.7  m/km),  pools,  fine 
substrates,  and  stable  banks.  The  middle  three  reaches  are  steeper  (1.3 
to  2.4  m/km)  with  a series  of  riffles  and  pools,  variable  substrates 
depending  on  current  regime,  and  unstable  banks.  The  two  downstream 
reaches  of  the  river  flow  through  a canyon  with  high  (40  to  50  m)  valley 
walls.  The  current  regime  is  swirling  and  the  streambed  is  composed  of 
coarse  substrates.  Bitumen  is  present  in  valley  walls  and  stream 

substrates. 

The  Ells  River  is  variable  in  slope  throughout  its  length 
(Sekerak  and  Walder  1980).  The  gradient  of  the  headwater  reaches 
(ca.  20  km)  which  drain  the  Birch  Mountains  are  very  steep  (5.0  to 

7.2  m/km)  with  a series  of  riffles  and  pools  as  well  as  coarse 
substrates.  In  the  mid  portion  of  the  Ells  River  (ca.  130  km)  the 
gradient  is  reduced  (1.9  to  4.1  m/km),  but  erosional  conditions 
(riffles,  coarse  substrates)  prevail.  There  is  a shift  to  depositional 
habitats  in  the  downstream  section  of  the  river. 

Corkum  et  al.  (1982)  outlined  the  reach  characteristics  of 
streams  draining  the  east  slope  of  the  Birch  Mountains.  The  headwaters 
and  mid  regions  of  the  Pierre  and  Asphalt  rivers  are  characterized  by 
erosional  habitats  (steep  gradient,  coarse  substrates,  riffles).  In 
contrast,  upper  Eymundson  Creek  and  the  Calumet  River  are  depositional 
for  most  of  their  lengths.  The  channel  of  the  Calumet  is  ill -defined 
and  beaver  activity  has  impounded  large  areas  of  surrounding  forests. 
All  streams  are  depositional  at  their  confluence  with  the  Athabasca 
River. 

6.2  BASELINE  DATA 

6.2.1  Discharge 

Of  the  several  gauging  stations  established  within  the  western 
region,  three  were  selected  for  study  (Table  11;  Figure  4).  The 
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Figure  29.  Elevation  profiles  of  rivers  in  the  western  region  . 
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Table  11.  Water  Survey  of  Canada  river  gauging  stations  in  the  west. 


Station 

Description 

Drainage 
Area  (km2) 

Gauge  Location 
(Lat  N/Long  W) 

Period  of 
Record 

Gauge 

Type3 

07DA005 

Beaver  River  near  454 

57 

06 

00 

1961-66 

RS  Reg  75 

Fort  MacKay 

111 

38 

00 

1972-74 

MC 

1975 

M 

07DA007 

Poplar  Creek 

151 

56 

54 

50 

1972 

M Reg  76 

111 

27 

35 

1973-83 

RC 

07DA010 

Ells  River  below 

1380 

57 

22 

30 

1975-78 

RC 

Gardiner  Lakes 

112 

33 

40 

07DA011 

Unnamed  Creek 

274 

57 

39 

31 

1975-79 

RC 

near  Fort  McMurray 

111 

31 

11 

1980-83 

RS 

07DA012 

Asphalt  Creek 

148 

57 

32 

20 

1975-77 

RC 

111 

40 

36 

07DA013 

Pierre  River 

123 

57 

27 

55 

1975-77 

RC 

111 

39 

14 

07DA014 

Calumet  River 

183 

57 

24 

12 

1975-77 

RC 

111 

40 

57 

07DA015 

Tar  River 

301 

57 

21 

14 

1975-77 

RC 

111 

45 

29 

07DA016 

Joslyn  Creek 

257 

57 

16 

27 

1975-79 

RC 

111 

44 

30 

1980-83 

RS 

07DA017 

Ells  River  near 

2410 

57 

16 

04 

1973-83 

RC 

mouth 

111 

42 

51 

07DA018 

Beaver  River 

165 

56 

56 

29 

1975-77 

RC 

above  Syncrude 

111 

33 

54 

1979-83 

RC 

07DA019 

Tar  River 

103 

57 

29 

05 

1976-77 

RC 

(upper  station) 

112 

01 

10 

07DB001 

MacKay  River 

5550 

57 

12 

38 

1972-83 

RC 

near  Fort  MacKay 

111 

41 

36 

07DB002 

Dover  River  near 

963 

57 

10 

12 

1975-77 

RC 

mouth 

111 

47 

38 

dR  = recording  gauge;  S = seasonal  operation;  C = continuous  operation. 
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selected  stations  were  downstream  locations  on  the  Poplar,  MacKay,  and 
Ells  rivers  where  flow  was  measured  using  a continuous  recording  gauge. 

Mean  daily  discharge  (1  in  5 days)  is  presented  for  Poplar 
Creek  (1976  to  1982),  the  MacKay  River  (1976  to  1980)  and  the  Ells  River 
(1976  to  1978)  (Figures  30,  31,  and  32).  The  period  of  discharge  record 
was  selected  to  correspond  to  the  water  quality  sampling  regime. 

Mean  annual  discharge  for  Poplar  Creek  ranged  from  0.35  to 
1.79  m3/s  during  a seven-year  period  (1976  to  1982).  Maximum  mean 
monthly  flows  occurred  during  late  summer  in  1976  (September, 

5.94  m3/s),  1978  (September,  9.20  m3/s)  and  1980  (August,  5.54  m3/s). 
Maximum  mean  monthly  flows  during  spring  runoff  ranged  from  0.72  to 
9.93  m3/s  during  the  seven  years.  There  is  no  flow  during  winter. 

Since  the  Poplar/Beaver  Creek  diversion,  the  discharge  regime 
at  Poplar  Creek  is  typically  bimodal.  During  the  open  water  season, 
saline  groundwater  that  enters  the  Beaver  Creek  Reservoir  is  diluted  and 
eventually  flows  into  the  Athabasca  River  via  Poplar  Creek.  In  the 
fall,  water  levels  in  the  Beaver  Creek  Reservoir  are  lowered,  increasing 
flow  in  Poplar  Creek.  Subsequently,  no  reservoir  (saline)  release 
occurs  over  the  winter  months  despite  the  continous  groundwater 
discharge  to  the  reservoir. 

The  MacKay  River  watershed  has  the  largest  drainage  area 
(5550  km2)  within  the  western  region.  The  mean  annual  discharge  ranged 
from  5.89  to  18.17  m3/s  during  the  five-year  period  (1976  to  1980). 
Maximum  mean  monthly  flows  ranged  from  14.57  to  65.94  m3/s  over  the  five 
years.  Peak  flows  occurred  as  a result  of  spring  runoff  (1976,  1979)  or 
summer  rainstorms  (1977,  1978,  1980).  Minimum  winter  flows  ranged  from 
0.2  to  0.4  m3/s. 

Maximum  mean  monthly  flow  in  the  Ells  River  occurred  in  April 
(21.5  m3/s) , July  (7.81  m3/s),  and  May  (16.22  m3/s)  in  1976,  1977,  and 
1978,  respectively.  Minimum  daily  discharge  ranged  from  0.40  to 
0.91  m3/s  during  the  three  years. 

Flow  from  most  of  the  streams  originating  in  the  Birch 
Mountains  (Tar,  Calumet,  Pierre,  Asphal t/Eymundson)  enters  the  Athabasca 
mainstem  directly.  Joslyn  Creek  enters  the  Ells  River  before  draining 
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Figure  30.  Daily  discharge  ( I in  5 days)  in  Poplar  Creek  ( WSC  Station  No.  07DA007), 
1976-1982  . 
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into  the  mainstem.  Patterns  of  discharge  regime  for  the  Pierre  River 
and  Asphalt  Creek  show  that  flow  contribution  of  these  streams  is  minor 
compared  to  other  rivers  of  the  western  region  (Figures  33  and  34). 

6.2.2  Water  Quality 

The  description  of  all  water  quality  sampling  sites  and  period 
of  record  are  presented  in  Table  12.  A tabular  summary  of  water  quality 
constituents  is  presented  for  20  sites  for  which  there  were  sufficient 
data  (more  than  10  samples)  (Appendix  10.5,  Tables  63  to  65). 

6.2.2. 1 Routine  parameters.  The  surface  waters  of  the  western  region 
are  neutral  to  alkaline  (Appendix  10.5,  Table  65).  Most  surface  waters 
in  the  region  are  classified  as  hard  (121  to  180  ug/L)  to  very  hard 
(>180  mg/L)  (McNeely  et  al . 1979).  Softer  waters  were  characteristic  of 
sites  along  the  Ells  River  and  Thickwood  Creek  as  indicated  by  median 
hardness  values  ranging  from  59  to  90  mg/L. 

Total  salt  concentration  as  a measure  of  specific  conductivity 
was  exceptionally  high  in  Poplar  Creek,  Beaver  River,  Calumet  River,  and 
Eymundson  Creek  compared  to  other  rivers  in  the  region  (Figures  35,  36, 
and  37).  The  elevated  saline  concentrations  in  these  rivers  are  a 
function  of  groundwater  inputs. 

Dissolved  sodium  chloride  salts  were  dominant  in  the  Beaver 
River,  Calumet  River,  and  Poplar  Creek.  Sodium  and  calcium  bicarbonates 
were  prevalent  in  the  Dover  and  MacKay  rivers;  magnesium  and  bicarbonate 
were  dominant  in  Eymundson  Creek.  In  other  streams  of  the  region, 
calcium  bicarbonate  was  the  major  dissolved  salt. 

Suspended  particulate  material  in  the  water  was  measured  by 
turbidity  and  non-fi 1 terable  residue  (NFR)  tests.  Streams  draining  the 
Birch  Mountains  had  elevated  median  values  of  these  constituents 
(Appendix  10.5,  Table  65).  Bank  slumping  along  these  streams 
contributes  to  the  suspended  particulates.  Elevated  maximum  levels  of 
turbidity  and  NFR  were  noted  at  sites  on  the  Beaver  River,  Bridge  Creek, 
and  the  MacKay  River  (Hwy.  63),  but  likely  were  related  to  soil  erosion 
as  a result  of  industrial  activity  or  road  construction. 
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Figure  33.  Daily  discharge  ( I in  5 days  ) in  the  Asphalt 
(WSC  Station  No.  DA0I20  ) , 1975-1977  . 


Figure  34.  Daily  discharge  ( I in  5 days)  in  the  Pierre 
(WSC  Station  No.  DA0I30),  1975-1977. 
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Table  12.  Description,  period  of  record,  and  number  of  samples  obtained  for 
water  quality  sampling  stations  within  the  western  region.  See 
Figure  4 for  location  of  sites. 


Latitude  N 


Code 

Longitude  W 

Period  of 

Number  of 

Watershed 

(00AT07-) 

( D 

M 

S ) 

Description 

Record 

Samples 

Poplar  Creek 

DA0070 

56 

54 

50 

Poplar  Creek 

1976-83 

39a 

111 

27 

16 

at  Hwy.  63 

Beaver  River 

DA0179 

57 

0 

58 

Beaver  River 

1976 

2 

Bridge  Creek 

111 

37 

28 

inside 

Syncrude 

DA0180 

56 

56 

29 

Beaver  River 

1977-78 

18a 

111 

33 

54 

above  Syncrude 

DAO 181 

57 

6 

54 

Beaver  River 

1976-78 

21a 

111 

37 

22 

at  Hwy.  63 

DA0182 

57 

7 

17 

Bridge  Creek 

1976-80 

32a 

111 

37 

30 

Diversion  at 
Hwy.  63 

DA0183 

57 

7 

9 

Bridge  Creek 

1980-81 

17a 

111 

36 

20 

before  confl . 
Athabasca 

Mac Kay 

DB0011 

57 

12 

38 

MacKay  River 

1976-80 

40a 

111 

41 

36 

at  Hwy.  63 

DB0020 

57 

10 

12 

Dover  River 

1976-78 

17a 

111 

47 

38 

3.2  km  above 
confl . MacKay 

DB0030 

56 

51 

20 

Dunkirk  River 

1976-79 

29a 

112 

42 

40 

near 

Ft.  MacKay 

DB0040 

56 

53 

55 

Thickwood 

1976-78 

15a 

112 

10 

15 

Creek  above 
confl . MacKay 

Ells 

DA0098 

57 

18 

23 

Ells  River 

1976-79 

13a 

111 

40 

20 

near  the  mouth 

DA0100 

57 

22 

30 

Upper  Ells 

1976-79 

27a 

112 

33 

40 

River  below 

Gardiner  L. 
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Table  12.  Continued. 


Watershed 


Code 

(00AT07-1 ) 


Latitude  B 
Longitude  W 


Period  of  Number 


( D 

M 

S ) 

Description 

Record 

of  Sampl 

DA0160 

57 

16 

27 

Joslyn  Creek 

1976-79 

22* 

111 

44 

30 

3.2  km  above 
Ells  River 

DA0170 

57 

16 

4 

Lower  Ells 

1976-78 

18* 

111 

42 

51 

River  3.2  km 
above  conf 1 . 
Joslyn 

DA0110 

57 

39 

31 

Unnamed  Creek 

1976-79 

17* 

111 

31 

11 

north  of  Fort 
McMurray 

DAO 111 

57 

39 

31 

Small  trib. 

1976 

1 

111 

21 

12 

to  Unnamed 
Creek 

DA0115 

57 

39 

30 

Unnamed  Creek 

1976 

2 

111 

31 

14 

near  the  mouth 

DA0150 

57 

21 

14 

Lower  Tar 

1976-78 

15* 

111 

45 

29 

River 

DA0151 

57 

19 

20 

Tar  River 

1976 

3 

111 

40 

55 

near  mouth 

DA0190 

57 

29 

5 

Upper  Tar 

1976-78 

12* 

112 

1 

10 

River 

DA0130 

57 

27 

55 

Pierre  River 

1976-78 

- 12* 

111 

39 

14 

DA0131 

57 

27 

1 

Pierre  River 

1976 

4 

111 

37 

29 

near  mouth 

DA0120 

57 

32 

20 

Asphalt  Creek 

1976-77 

10* 

111 

40 

36 

near  Fort  MacKay 

DA0121 

57 

33 

10 

Eymundson  Ck 

1976-77 

11* 

111 

39 

20 

upstream  of 
Asphalt 

DA0122 

57 

29 

27 

Eymundson  at 

1976 

1 

111 

34 

7 

the  mouth 

DA0140 

57 

24 

12 

Calumet  River 

1976-78 

13* 

111 

40 

57 

near  Fort 

Rivers 
draining 
east  slopes 
of  Birch 
Mountains 


Mac Kay 
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Table  12.  Concluded. 


Watershed 

Code 

(00AT07-1 ) 

Latitude  B 
Longitude  W 
( D MS) 

Description 

Period 

Record 

of  Number 
of  Samples 

DA0141 

57  24  38 

111  39  57 

Calumet  River 
near  the  mouth 

1976 

6 

a 

Sites  with  sufficient  samples  to  be  used  in  analysis. 
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DA0098 

ELLS  RIVER  , NEAR  MOUTH 
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UPPER  ELLS  RIVER 


DA0I60 

JOSLYN  CREEK 


D&OI70 

LOWER  ELLS  RIVER 


Figure  36.  Maximum  and  minimum  values  ( me/L)  of  major 
ions  in  the  Ells  and  MacKay  drainages  . 
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DAOIIO 

UNNAMED  CREEK,  N.  OF  ASPHALT  CREEK 


DA0I20 

ASPHALT  CREEK 


DAOI2I 

EYMUNDSON  CREEK 


DA0I30 

PIERRE  RIVER 


DA0I40 

CALUMET  RIVER 


DA0I50 

LOWER  TAR  RIVER 


DA0I90 

UPPER  TAR  RIVER 


CATIONS  (me/L)  ANIONS  (me/L) 


Figure  37.  Maximum  and  minimum  values  (me/L)  of  major 
ions  in  streams  draining  the  east  slopes  of  the 
Birch  Mountains  . 
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Total  organic  carbon  is  a measure  of  biological  activity. 
Elevated  levels  were  prevalent  in  surface  waters  throughout  the  region, 
but  particularly  in  the  Calumet  and  Beaver/Bridge  Rivers  (Appendix  10.5, 
Table  65). 

The  presence  of  phenols  as  well  as  oil  and  grease  in  surface 
waters  is  the  result  of  industrial  and  municipal  activity.  Elevated 
levels  of  these  parameters  occurred  at  river  sites  associated  with  oil 
sands  development.  Median  values  of  phenolic  material  for  Poplar  Creek 
(DA0070)  and  Bridge  Creek  (DA0183)  were  0.003  mg/L  and  0.006  mg/L, 
respectively . Median  values  for  all  other  river  sites  were  less  than  or 
equal  to  0.001  mg/L.  The  median  value  of  oil  and  grease  for  Bridge 
Creek  (DA0183)  was  1.4  mg/L.  Median  values  for  all  other  river  sites 
were  less  than  1.0  mg/L. 

6. 2. 2. 2 Nutrients . Major  forms  of  phosphorus  and  nitrogen  are 
presented  in  Appendix  10.5  (Table  66). 

Elevated  concentrations  of  total  phosphorus  were  exhibited  in 
the  Asphalt,  Eymundson,  Calumet,  lower  Tar,  and  Josyln  Rivers.  The 
median  (113  to  150  ug/L)  and  maximum  (360  to  720  ug/L)  values  of  total 
phosphorus  for  these  streams  draining  the  Birch  Mountains  may  reflect 
erosional  processes  within  the  watershed.  Similarly  high  concentrations 
of  total  phosphorus  were  recorded  for  the  Beaver  and  Dunkirk  rivers 
(Appendix  10.5,  Table  66). 

Orthophosphate  phosphorus  (PO^-P),  an  indicator  of  enrichment, 
exhibited  the  highest  median  values  in  the  Calumet  (37  ug/L)  and  MacKay 
(24  ug/L)  rivers.  Median  values  at  all  other  sites  ranged  from  less 
than  10  to  20  ug/L.  The  highest  maximum  values  were  recorded  at  sites 
on  the  Beaver  River  above  Syncrude  Canada  Ltd.  (235  ug/L)  as  well  as  on 
the  Dunkirk  (180  ug/L),  MacKay  (100  ug/L),  and  Joslyn  (110  ug/L)  rivers. 

Total  kjeldahl  nitrogen  (TKN)  is  a measure  of  inorganic 
(ammonia)  and  organic  nitrogen.  Median  TKN  values  for  rivers  in  the 
western  region  ranged  from  0.78  to  2.45  mg/L.  Sites  with  the  highest 
median  TKN  concentrations  were  located  on  the  Calumet  (2.45  mg/L)  and 
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Eymundson  (1.86  mg/L)  rivers.  Maximum  TKN  concentrations  ranged  from 
1.68  to  7.85  mg/L.  High  concentrations  were  recorded  for  the  upper  Ells 
(7.85  mg/L),  Asphalt  (5.60  mg/L),  Eymundson  (5.10  mg/L),  lower  Tar 
(4.40  mg/L),  and  Calumet  (4.23  mg/L)  rivers. 

Median  ammonia  concentrations  ranged  from  0.03  to  0.35  mg/L. 
Sites  on  the  Calumet  River  (0.35  mg/L)  and  Beaver  River  at  Hwy.  63 
(0.19  mg/L)  exhibited  the  highest  median  concentrations  of  ammonia. 
Maximum  concentrations  of  ammonia  ranged  from  0.05  to  1.14  mg/L,  with 
sites  on  Poplar  Creek  (1.14  mg/L)  and  the  Calumet  River  (0.88  mg/L) 
exhibiting  the  highest  values. 

Nitrate/nitrite  nitrogen  is  another  form  of  inorganic 
nitrogen.  Median  concentrations  of  nitrate/nitrite  nitrogen  ranged  from 
less  than  0.01  to  0.09  mg/L.  Sites  on  the  lower  (0.09  mg/L)  and  upper 
(0.06  mg/L)  Ells  River  exhibited  highest  median  values.  The  maximum 
concentrations  of  nitrate/nitrite  nitrogen  ranged  from  0.06  to  1.0  mg/L. 
Highest  values  were  recorded  in  Thickwood  Creek  (1.0  mg/L),  a tributary 
of  the  Mackay  River,  and  the  MacKay  mainstem  (0.69  mg/L). 

The  east-draining  rivers  of  the  Birch  Mountains  were  more 
productive  in  terms  of  total  phosphorus  and  organic  nitrogen  than  other 
rivers  of  the  western  region.  In  contrast,  the  Beaver/Poplar  and  MacKay 
drainages  exhibited  higher  concentrations  of  P0^-P  and  inorganic 
nitrogen  than  the  other  rivers. 

6. 2. 2. 3 Metals.  Summary  data  of  metals  in  surface  waters  of  the 
western  region  are  presented  in  Appendix  10.5  (Table  67).  Of  these, 
seven  are  considered  to  be  potentially  toxic  when  present  in  excessive 
amounts  including  arsenic,  cadmium,  lead,  mercury,  selenium,  and  silver. 

The  median  range  of  arsenic  concentrations  in  the  western 
rivers  was  from  less  than  0.05  to  1.2  ug/L,  well  below  the  surface  water 
objective  of  10  ug/L.  Maximum  values  ranged  from  1.5  to  17.3  ug/L. 
Lower  Tar  River  (17.3  ug/L),  Poplar  Creek  (16  ug/L)  and  the  MacKay  River 
at  Hwy.  63  (13  ug/L)  exhibited  maximum  values  which  exceeded  the 
objective. 
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The  ASWO  for  cadmium  is  10  ug/L.  The  median  values  are  less 
than  or  equal  to  1 ug/L  for  all  rivers.  Maximum  values  for  most  sites 
ranged  from  1 to  6 ug/L.  One  site.  Poplar  Creek  (16  ug/L),  exceeded 
the  ASWO. 

The  median  values  of  chromium  ranged  from  less  than  3 to 
4 ug/L,  well  below  the  objective  of  50  ug/L.  However,  the  range  of 
maximum  values  of  chromium  was  5 to  205  ug/L.  The  three  sites  which 
exceeded  the  objectives  were  located  within  the  Beaver/Poplar  drainage 
(Appendix  10.5,  Table  67). 

Median  values  for  lead  concentrations  were  less  than  or  equal 
to  2 ug/L  for  all  sites.  Maximum  values  for  19  of  20  sites  ranged  from 
3 to  48  ug/L.  The  maximum  lead  concentration  at  Poplar  Creek  (175  ug/L) 
exceeded  the  ASWO  of  50  ug/L. 

Median  values  of  mercury  concentration  at  all  sites  were  less 
than  0.1  ug/L,  below  the  ASWO  of  1 ug/L.  Maximum  values  of  mercury 
concentration  which  exceeded  the  objective  were  recorded  at  five  sites 
including  the  Bridge  (6.4  ug/L),  Calumet  (4.6  ug/L),  Eymundson 
(2.7  ug/L),  MacKay  (2.1  ug/L),  and  Poplar  (1.6  ug/L)  rivers. 

Dissolved  selenium  concentrations  for  median  and  maximum 
values  at  all  20  river  sites  were  below  the  ASWO  of  10  ug/L. 

Median  values  of  silver  for  river  sites  within  the  region  were 
less  than  1 to  2 ug/L.  The  range  of  maximum  concentrations  for  18  of 
20  sites  was  1 to  9 ug/L.  Poplar  Creek  (28  ug/L)  and  the  lower  Tar 
(70  ug/L)  exhibited  high  silver  concentrations  compared  to  the  other 
18  sites.  Only  the  lower  Tar  River  exceeded  the  ASWO  of  50  ug/L. 

In  summary,  the  median  concentrations  of  potentially  toxic 
substrates  were  below  the  ASWO  for  all  rivers  examined  within  the 
western  region.  Of  the  seven  metals,  dissolved  selenium  was  the  only 
constituent  with  maximum  values  below  the  ASWO.  In  several  instances 
the  objectives  were  exceeded.  Poplar  Creek  is  of  interest  because 
maximum  concentrations  exceeded  the  ASWO  for  five  potentially  toxic 
metals  including  arsenic,  cadmium,  chromium,  lead,  and  mercury. 
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Other  non-toxic  metals  detected  in  surface  waters  are 
frequently  related  to  weathering  of  rock  material.  This  group  consists 
of  aluminum,  boron,  cobalt,  copper,  iron,  manganese,  nickel,  vanadium, 
and  zinc. 

Median  values  of  aluminum  concentrations  in  western  rivers 
ranged  from  40  to  328  ug/L  with  the  exception  of  Asphalt  Creek  (1095 
ug/L).  Maximum  concentrations  ranged  from  210  to  14  700  ug/L.  Since 
aluminum  is  not  considered  to  be  a health  problem  in  surface  waters, 
objectives  have  not  been  established. 

The  median  values  of  boron  concentrations  ranged  from  60  to 
340  ug/L,  below  the  ASW0  of  500  ug/L.  The  maximum  values  ranged  from 
100  to  1080  ug/L.  Sites  with  boron  concentrations  exceeding  the  ASW0 
were  located  within  the  Beaver/Poplar,  MacKay,  and  Birch  Mountain  east 
slope  drainages  (Appendix  10.5,  Table  67). 

No  surface  water  objectives  exist  for  cobalt.  Median  values 
ranged  from  less  than  1 to  2 ug/L.  Maximum  values  at  most  sites  within 
the  region  ranged  from  1 to  14  ug/L.  Asphalt  (28  ug/L)  and  Poplar  (31 
ug/L)  creeks  exhibited  higher  values. 

Copper  concentrations  (median  values)  ranged  from  1 to  8 ug/L, 
below  the  ASW0  of  20  ug/L.  Maximum  values  at  most  sites  ranged  from 
4 to  33  ug/L  with  two  sites,  lower  Tar  River  (114  ug/L)  and  the  MacKay 
River  (180  ug/L),  exhibiting  notably  higher  concentrations. 

The  ASW0  for  iron  (300  ug/L)  and  manganese  (50  ug/L)  were 
established  for  aesthetic  reasons  because  of  taste  and  staining 
properties  of  the  metals  at  high  concentrations.  Median  (range  300  to 
4325  ug/L)  and  maximum  (range  1650  to  52  500  ug/L)  values  of  iron 
exceeded  the  objectives  at  all  river  sites.  Median  manganese  values 
ranged  from  19  to  370  ug/L.  Here,  14  of  20  sites  representing  streams 
draining  the  Birch  Mountains  and  the  Beaver/Poplar  watershed  exceeded 
the  objectives.  Maximum  values  (range  69  to  6550  ug/L)  exceeded  the 
objectives  at  all  sites. 

Nickel,  vanadium,  and  zinc  are  frequently  associated  with  oil 
sands.  No  ASW0  exist  for  nickel  or  vanadium.  Median  values  of  nickel 
concentrations  ranged  from  <2  to  12  ug/L.  Maximum  nickel 
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concentrations  ranged  from  <2  to  46  ug/L,  with  Asphalt  and  Eymundson 
Creeks  exhibiting  the  highest  values.  Median  values  of  vanadium  were 
low  (range <1  to  1 ug/L).  Maximum  concentrations  ranged  from  1 to 
24  ug/L,  with  the  Beaver  River  at  Hwy.  63  exhibiting  the  highest  value. 

The  ASWO  for  zinc  is  50  ug/L.  Median  values  for  all  sites 
(range  3 to  30  ug/L)  were  within  the  objectives.  Maximum  values  (24  to 
155  ug/L)  exceeded  the  objectives  at  sites  within  all  watersheds 
throughout  the  region. 

6.2.3  Benthic  Invertebrates 

A summary  of  benthic  invertebrate  studies,  which  were 
conducted  in  the  Beaver/Poplar  and  MacKay  rivers  and  streams  of  the 
Birch  Mountain  drainages,  is  presented  in  Table  13.  Sampler  type,  mesh 
size,  and  sample  frequency  differed  among  studies  conducted  in  the  three 
drainages. 

Surveys  conducted  on  the  MacKay  River  (McCart  et  al.  1978)  and 
streams  of  the  Birch  Mountains  (Corkum  et  al . 1982)  were  baseline 
studies.  These  studies  provide  a base  against  which  future  changes  can 
be  compared.  The  benthic  surveys  conducted  on  the  Beaver/Poplar 
drainage  to  determine  the  effects  of  stream  discharge  on  the  biota 
(Noton  and  Chymko  1978;  Boerger  1983,  1984)  were  impact  studies. 

Taxonomic  differences  of  invertebrates  among  sample  sites  in 
the  baseline  studies  reflected  differences  in  habitat  type  (erosional 
versus  depositional ) . Corkum  et  al.  (1982)  categorized  organisms  into 
taxonomic  and  functional  feeding  groups  and  related  the  groups  to 
habitat  types.  The  relative  abundance  of  invertebrates  was  more  similar 
among  sites  within  the  same  habitat  than  among  sites  within  the  same 
drainage. 

Noton  and  Chymko  (1978)  summarized  changes  in  aquatic 
resources  of  water  bodies  associated  with  the  Beaver/Poplar  Creek 
diversion  scheme.  The  initial  effect  of  the  spillway  discharge  was  an 
increase  in  the  relative  abundance  of  benthic  invertebrates  in  Poplar 
Creek  downstream  from  the  confluence  with  the  spillway. 
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Habitat;  E,  Erosional;  D,  depositional . 

Artificial  sampler:  cylindrical  wire  barbecue  baskets  (length:  27  cm,  diameter:  18  cm), 

each  filled  with  20  to  30  rocks. 

Artificial  sampler:  rock-filled  wire  baskets  (McCart  et  al . 1977). 
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Boerger  (1983,  1984)  showed  that  following  the  initial 
enhanced  response  by  the  biota  to  the  spillway  discharge,  the  rate  of 
increase  invertebrate  density  and  taxa  declined.  Boerger  (1983) 
suggested  that  increases  in  biota  were  a function  of  increased  flow, 
temperature,  and  detritus  from  the  epi limnetic  discharge  over  the 
spillway  as  well  as  substrate  changes  resulting  from  channelization. 
Saline  mine  depressurization  water  released  into  Beaver  Creek  Reservoir 
had  no  significant  effect  on  the  benthos  in  Poplar  Creek. 

6.3  RELATIONSHIPS  AMONG  PARAMETERS 

To  select  possible  water  quality  parameters  which  may  have 
significant  associations  with  discharge,  a correlation  analysis  was 
performed  between  conductivity  (a  known  correlate  of  discharge)  and 
other  water  quality  constituents.  Variables  that  would  be 
autocorrelated  with  conductivity  (ions)  were  not  included  in  the 
analysis. 

Results  from  the  correlation  analysis  using  median  values  of 
water  quality  constituents  showed  that  only  total  inorganic  carbon  (TIC) 
was  significantly  correlated  with  conductivity  at  the  0.01  level. 

Relationships  among  conductivity,  TIC,  and  discharge  were 
estimated  using  data  from  all  sampling  dates  for  the  Poplar,  Ells,  and 
MacKay  rivers.  Regression  analyses  were  performed  using  linear, 
log/linear,  and  log/log  transformations  of  the  data.  Equations  from 
those  transformations  yielding  the  best  fit  (i.e.,  least  variable  and 
highest  R2)  were  selected  and  are  presented  in  Table  14. 

There  was  an  inverse  significant  relationship  between  both  TIC 
and  conductivity  with  discharge  for  all  three  drainages.  Thus,  levels 
of  conductivity  and  TIC  were  high  when  discharge  levels  were  low. 

Because  the  residual  variances  among  the  rivers  were 
heterogeneous  (Bartlett's  test,  p<0.05),  comparisons  among  intercepts 
and  slopes  of  the  three  regressions  were  considered  inappropriate. 

There  was,  however,  a trend  for  conductivity  levels  to 
decrease  more  rapidly  with  discharge  in  Poplar  Creek  and  the  MacKay 


123 


Table  14.  Summary  of  regression  equations  of  conductivity  (SC)  and  total 

inorganic  carbon  (TIC)  with  discharge  for  Poplar  Creek,  and  Ells  and 
MacKay  rivers.  Equations  with  the  best  fit  of  the  data  were  of  the 
following  log/log  format: 

Y = AXB 

where:  A = intercept,  B = slope,  Y = dependent  variable  (water 

quality),  X = independent  variable  (discharge). 


River 

df 

(n-2) 

Y 

A 

B 

R2 

p 

(Ho:  B:  = 0) 

Poplar 

57 

SC 

723.4 

-0.588 

.31 

L0.001 

(DA0070) 

27 

TIC 

38.6 

-0.227 

.32 

L0.001 

Ells 

16 

SC 

311.1 

-0.302 

.73 

L0.001 

(DA0170) 

16 

TIC 

27.5 

-0.289 

.59 

L0.001 

MacKay 

38 

SC 

801.9 

-0.465 

.83 

L0.001 

(DA0140) 

32 

TIC 

76.4 

-0.465 

.86 

L0.001 

L = less  than. 
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River  than  in  the  Ells  River.  The  TIC  appeared  to  decrease  more  rapidly 
with  discharge  in  the  MacKay  River  than  in  the  other  two  rivers. 
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7.  THE  SOUTH 

The  southern  region  encompasses  an  area  of  4700  km2  extending 
south  of  Fort  McMurray  from  the  Athabasca  and  Clearwater  rivers  to  the 
AOSERP  boundary.  There  are  ten  rivers,  several  tributaries,  and  lakes 
within  the  region.  A detailed  report  on  Gregoire  Lake  is  in  progress 
(Akena  in  prep.)  and  will  not  be  considered  in  this  report.  This  study 
focuses  on  six  major  rivers  (Horse,  Hangingstone,  Saline,  Surmount, 
Christina,  and  Clearwater)  for  which  there  is  an  adequate  water  quality 
data  set.  References  will  be  made  to  other  streams  in  the  region 
depending  on  available  data. 

7.1  BASIN  FEATURES 

7.1.1  Land  Formations 

Land  formations  of  the  southern  region,  adapted  from  Turchenek 
and  Lindsay  (1982),  are  presented  in  Figure  38.  The  region  occurs 
within  the  Boreal  Mixedwood  Forest,  much  of  which  is  poorly  drained  low 
lying  muskeg. 

The  bedrock  geology  of  the  Athabasca  and  Clearwater  river 
valleys  is  characterized  by  limestone  and  shale.  This  narrow  band  at 
the  northern  boundary  of  the  region  constitutes  the  Waterways  Formation 
of  the  Devonian  period  (Paleozoic  era).  Most  of  the  southern  region, 
however,  is  characterized  by  younger  Cretaceous  rocks  (sandstone  and 
shale)  of  the  Mesozoic  era. 

Most  of  the  tributaries  in  the  south  are  low  gradient  rivers 
that  meander  through  flood  plains.  The  smaller  tributaries,  which  drain 
muskeg  areas,  consist  of  a series  of  beaver  ponds  (Tripp  and  Tsui  1980). 
The  downstream  reaches  of  the  six  study  rivers  lie  within  the  Methy 
Portage  Plain,  an  undulating  plain  with  low  to  moderate  relief.  White 
spruce  and  aspen  forests  are  most  common  in  the  upland  regions  of  the 
plain.  Wetlands  of  low  lying  areas  are  characterized  by  black  spruce 
and  organic  soils. 

Where  the  Horse  River  enters  the  Athabasca  mainstem  and  the 
Hangingstone  and  Christina  join  the  Clearwater,  the  channels  cut  through 


Figure  38. 


Land  formations  of  the  southern  region  . 

( odoptod  from  Tarehonok  and  Linduoy  !••£) 
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steep  valleys.  Shrubs  occupy  the  bottom  lands;  spruce  and  aspen,  the 
valley  walls.  No  vegetation  occurs  on  the  steep  valley  walls. 

7.1.2  Watershed 

Major  drainage  systems  in  the  region  include  the  Clearwater, 
Christina,  Hangingstone,  and  Horse  rivers  (Figure  39).  The  following 
description  of  these  rivers  is  summarized  from  Griffiths  1973;  Northwest 
Hydraulics  Consulting  Ltd.  1975;  Tripp  and  McCart  1979;  and  Tripp  and 
Tsui  1980. 

The  Clearwater  River  drainage  basin  is  31  598  km2  of  which 
16  572  km2  occurs  within  Alberta  (Northwest  Hydraulics  Consulting  Ltd. 
1975).  The  Clearwater  flows  a distance  of  about  65  km  within  the  AOSERP 
study  area  before  its  confluence  with  the  Athabasca  River  at  Fort 
McMurray . 

The  large  number  of  lakes  within  the  Clearwater  basin  reduces 
runoff  and  moderates  flood  peaks  (Northwest  Hydraulics  Consulting  Ltd. 
1975).  However,  during  spring  breakup,  ice  jams  may  develop  on  the 
Athabasca  and  Clearwater  rivers,  causing  flood  conditions  in  Fort 
McMurray . 

The  Christina  River,  a major  tributary  of  the  Clearwater 
River,  drains  an  area  of  13  380  km2  and  flows  a distance  of  326  km 
(Griffiths  1973;  Northwest  Hydraulics  Consulting  Ltd.  1975).  Only  the 
downstream  27  km  of  the  Christina  River  occurs  within  the  AOSERP  study 
area. 

The  Gregoire  River  is  the  major  tributary  of  the  Christina 
River  within  the  AOSERP  study  area.  The  Gregoire  flows  83  km  from 
Gregoire  Lake  to  the  confluence  of  the  Christina  River,  draining  978  km2 
of  the  Christina  basin.  The  Gregoire  River  meanders  over  fine  substrate 
until  near  its  mouth,  where  it  cuts  through  a steep  valley  to  the 
Christina  River. 

The  Hangingstone  River  basin  drains  an  area  of  1074  km2 
(Northwest  Hydraulics  Consulting  Ltd.  1975).  The  river  originates  near 
the  southern  boundary  of  the  AOSERP  study  area  and  flows  98  km  before 
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Figure  39.  South  surface  water  drainage. 

(adapted  frem  Neill  end  Evans  ) 
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joining  the  Clearwater  River  near  Fort  McMurray.  In  the  upper  reach, 
stream  gradient  is  reduced  and  the  river  meanders  over  fine  sandy 
substrates.  Stream  gradient  increases  again  near  the  mouth  as  the  river 
cuts  through  a steep  valley  to  the  floodplain  of  the  Clearwater  River. 

The  Horse  River  basin  drains  an  area  of  2154  km2,  most  of 
which  occurs  outside  of  the  AOSERP  study  area.  The  Horse  River 
originates  in  muskeg,  flowing  200  km  before  joining  the  Athabasca  River 
at  Fort  McMurray.  About  1 km  upstream  of  the  confluence,  outcrops  of 
the  McMurray  Formation  occur  along  the  valley,  forming  steep  bituminous 
banks. 


7.1.3  Physical  Description  of  Study  Streams 

One  important  aspect  of  the  study  streams  is  the  shift  from 
depositional  to  erosional  habitats  near  their  confluence  with  the 
Athabasca  and  Clearwater  mainstems.  The  increase  in  slope  near  the 
mouths  of  these  streams  is  observed  in  the  elevation  profiles  presented 
in  Figure  40. 

The  tributary  streams  (the  Surmount  River  being  the  exception) 
originate  in  the  plains  and  meander  over  mud  substrates  and  through 
beaver  ponds,  with  occasional  stretches  of  fast  flowing  waters  occurring 
over  coarse  substrates.  The  Hangingstone  has  a steeper  gradient  along 
its  length  than  the  other  tributaries.  Running  water  over  coarse 
substrates  in  confined  channels  near  the  mouth  of  the  tributary  streams 
is  more  typical  of  headwater  reaches  in  eastern  deciduous  forests  or  the 
foothills  of  the  Rocky  Mountains.  The  difference,  however,  is  in  the 
increased  stream  width,  higher  flows,  less  overhanging  vegetation,  and 
occasional  slumping  of  valley  walls  in  the  downstream  areas  of  the 
southern  AOSERP  tributaries. 

7.2  BASELINE  DATA 

7.2.1  Discharge 

The  Water  Survey  of  Canada  has  operated  gauging  or  stage 
recorders  at  six  locations  within  the  south  region,  including  four  sites 
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Figure  40.  Elevation  profiles  of  southern  tributaries, 

{ Adapted  from  Tripp  apd  Tsui  , I960  ) 
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on  the  Clearwater  River  and  one  on  each  of  the  Horse  and  Hangingstone 
rivers  (Table  15).  The  Hangingstone  near  Fort  McMurray  (07CD004)  and 
the  Clearwater  River  near  Draper  (07DC001)  were  selected  for  the  study 
because  continuous  recording  gauges  were  operated  over  the  longest 
period,  including  the  dates  when  water  samples  were  collected.  The 
location  of  the  Hangingstone  and  Clearwater  stations  are  indicated  on 
Figure  4. 

Mean  daily  discharge  (1  in  5 days)  is  presented  for  the 
Hangingstone  (1976  to  1982)  and  Clearwater  (1976  to  1981)  rivers  in 
Figures  41  and  42,  respecti vely. 

Discharge  levels  in  the  Hangingstone  River  typically  peaked  at 
spring  breakup,  decreased  over  the  summer  months,  and  increased  during 
autumn  freshets  before  declining  to  winter  base  flows.  Peak  discharge 
levels  were  typically  22  to  24  m3/s;  higher  levels  were  recorded  in  1976 
and  1980.  Winter  base  flows  ranged  from  0.1  to  0.4  m3/s. 

The  same  discharge  pattern  of  fluctuations  was  evident  in  the 
Clearwater  River.  Here,  maximum  discharge  levels  were  about  300  m3/s 
with  levels  of  up  to  600  m3/s.  In  1976  and  1977  summer  flows  also  were 
elevated.  Winter  base  flows  ranged  from  40  to  100  m3/s. 

7.2.2  Water  Quality 

The  description  and  period  of  record  for  the  water  quality 
sampling  stations  are  presented  in  Table  16.  The  station  locations  also 
are  indicated  in  Figure  4. 

Water  quality  data  are  summarized  for  the  seven  sampling 
stations  in  the  region,  1976  to  1983  (Appendix  10.6,  Tables  68  to  70). 
The  data  are  presented  and  discussed  by  parameter  group  (routine, 
nutrients,  and  metals). 

7.2.2. 1 Routine  parameters.  This  parameter  grouping  includes 
inorganic  ions,  suspended  paticulates,  organic  constituents,  phenols, 
and  oil  and  grease.  Summary  data  are  presented  in  Appendix  10.6 
(Table  66). 


Table  15.  Water  survey  of  Canada  gauging  stations  in  the  southern  region. 
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Figure  41.  Daily  discharge  (I  in  5)  in  the  Hangingstone  River,  1976-1982. 


WATER  QUALITY  SAMPLING  DATE 
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Figure  42.  Doily  discharge  ( I in  5 days)  in  the  Clearwater  River  at 
Draper,  1976  -1981. 
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Table  16.  Description  and  period  of  record  for  water  quality  sampling 
stations  within  the  southern  region. 


Latitude  N 


Code 

Longitude  W 

Period 

Number 

(00AT07- 

) (D 

M 

S) 

Description 

of  Record 

of  Sampli 

CC0030 

56 

42 

20 

Horse  R.  at  Abasands  Park 

1976-77 

5 

111 

23 

40 

3.2  km  above  confluence 
with  Athabasca  River 

CC0300 

56 

43 

06 

Horse  R.  near  Fort  McMurray 

1977-79 

15 

111 

23 

40 

100  m above  confluence  with 
Athabasca  River 

CD0040 

56 

42 

18 

Hangingstone  R.  at  Hwy.  63 

1976-83 

67 

111 

21 

20 

CD0100 

56 

42 

10 

Saline  Creek 

1976-77 

8 

111 

20 

38 

CD0100 

56 

41 

53 

Clearwater  R. 

1976-81 

50 

111 

19 

02 

2 km  above  Waterways 

CE2000 

56 

27 

01 

Surmount  Creek  about 

1978-83 

47 

111 

03 

45 

3.2  km  above  Gregoire 
Lake 

CE2191 

56 

39 

30 

Christine  R.  above  confluence 

1978-79 

10 

111 

02 

30 

with  Clearwater  R. 
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Alkalinity  (max.  range  191  to  407  mg/L)  and  hardness  (max. 
range  177  to  355  mg/L)  concentrations  were  elevated  in  the  Horse, 
Hangingstone,  Saline,  and  Christina  rivers.  Lower  levels  of  these 
parameters  occurred  in  the  Clearwater  (alkalinity  max.  95  mg/L;  hardness 
max.  109  mg/L),  and  Surmount  rivers  (alkalinity  max.  132  mg/L;  hardness 
max.  132  mg/L).  Maximum  levels  of  pH  ranged  from  8.0  to  9.6  for  all 
rivers. 

The  relative  dominance  of  major  ions  for  each  water  quality 
site  is  presented  in  Figure  43.  The  figure  represents  the  relative 
contribution  of  each  major  ion  and  an  approximation  of  the  total  ion 
content. 

The  ranking  of  ion  dominance  (measured  in  me/L)  in  five  of 
seven  rivers  of  the  south  region  is  as  follows: 

cations:  calcium  > magnesium  > sodium 

anions:  bicarbonate  >sulphate  >chloride 

Exceptions  include  the  Clearwater  and  Christina  rivers  where 
sodium  and  chloride  ions  predominated.  Sodium  ions  also  were  prevalent 
in  the  Hangingstone  River.  Total  ion  concentration  appears  to  be 
related  to  bedrock  geology  ( i . e . , specific  conductance  is  highest  in  the 
Hangingstone  River  and  Saline  Creek).  The  three  water  quality  sites 
with  elevated  sodium  and  chloride  ions  are  located  within  Devonian  rock 
formations. 

Ion  concentration  is  inversely  related  to  discharge.  The 
maximum  and  minimum  ion  values  presented  in  Figure  43  represent  two 
discharge  regimes.  Maximum  ion  concentrations  are  evident  during  winter 
base  flows.  Minimum  concentrations  occur  at  spring  runoff  or  during 
flood  conditions. 

Colour  is  a measure  of  water  clarity  and  has  no  direct 
chemical  significance  (Hem  1970).  The  brown  colour  of  rivers  in  the 
AOSERP  region  is  due  to  the  inflow  of  humic  and  fulvic  acids  from 
surrounding  muskeg.  The  maximum  colour  values  in  the  Hangingstone 
(160  mg/L),  Clearwater  (120  mg/L),  Surmont  (200  mg/L),  and  Christina 
(130  mg/L)  rivers  reflect  muskeg  characteristics  of  the  watershed. 
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Figure  43.  Maximum  and  minimum  values  (me/L) 
of  major  ions  in  rivers  of  the  south  . 
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Total  organic  carbon  is  related  to  the  biological  activity  in 
the  river.  No  differences  were  detected  among  these  productive  rivers 
(Appendix  10.6,  Table  68). 

Phenolic  compounds  at  very  high  levels  may  be  toxic  to  fish 
and  at  lower  levels  may  cause  fish  tainting  and  taste  problems.  Sources 
of  phenols  include  industrial  and  municipal  wastes  as  well  as  decaying 
plant  material.  Maximum  values  of  recorded  phenol  levels  in  the 
Hangingstone,  Saline,  Clearwater,  and  Christina  rivers  exceeded  the 
Alberta  surface  water  objective  of  0.005  mg/L. 

Sources  of  oil  and  grease  include  industrial  and  municipal 
wastes.  Maximum  values  (6.5  mg/L)  were  detected  in  the  Hangingstone 
River.  Surface  water  objectives  have  not  been  established  for  oil  and 
grease. 

7. 2. 2. 2 Nutrients.  Phosphorus  is  frequently  considered  to  be  the 
nutrient  that  controls  plant  life.  Total  phosphorus  and  orthophosphate 
phosphorus  (PO^-P)are  the  two  forms  of  phosphorus  presented  in 
Appendix  10.6  (Table  69). 

Total  phosphorus  values  were  high,  with  median  values 
exceeding  the  objective  of  50  ug/L  at  all  sites.  A maximum  value  of 
1300  ug/L  was  recorded  at  the  Hangingstone  River,  1980  August  04.  Most 
values  at  this  site  ranged  from  30  to  230  ug/L,  with  higher  values 
occurring  in  spring  during  runoff. 

Orthophosphates  are  the  potentially  available  form  of 
phosphorus.  Typically,  concentrations  are  lowest  in  summer  when  uptake 
by  plants  occurs.  Maximum  concentrations  occur  in  fall  and  winter  when 
plants  die.  Highest  maximum  values  were  recorded  in  the  Horse 
(100  mg/L),  Hangingstone  (123  ug/L),  and  Saline  (220  ug/L)  rivers. 

Ammonia,  nitrate/nitrite  nitrogen,  and  total  kjeldahl  nitrogen 
(TKN)  are  the  three  forms  of  nitrogen  presented  in  Appendix  10.6 
(Table  69).  Forms  of  inorganic  nitrogen  are  ammonia  and  nitrate/nitrite 
nitrogen.  Total  kjeldahl  nitrogen  is  total  organic  nitrogen  plus 
ammonia.  In  the  presence  of  oxygen,  ammonia  may  be  converted  to  nitrate 
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by  a process  of  nitrification.  Also,  nitrite  is  rapidly  converted  to 
nitrate  in  the  presence  of  oxygen. 

Only  the  Hangingstone,  Clearwater,  and  Surmount  sites  had 
sufficient  numbers  of  samples  to  denote  seasonal  trends.  No  detectable 
seasonal  trends  in  TKN  were  noted  among  the  sites.  Maximum  ammonia  and 
nitrate/nitrite  concentrations,  however,  occurred  in  winter  and  spring. 
The  highest  maximum  value  of  ammonia  was  recorded  in  the  Clearwater 
River  (0.73  mg/L).  The  highest  maximum  value  of  TKN  (9.15  mg/L)  and 
nitrate/nitrite  nitrogen  (0.66  mg/L)  was  recorded  in  the  Hangingstone 
River.  The  maximum  concentration  of  ammonia  recorded  in  the 
Hangingstone  River  was  0.45  mg/L. 

7. 2. 2. 3 Metals.  The  data  for  16  metal  parameters  are  summarized  in 
Appendix  10.6  (Table  70).  Median  values  of  nine  metals  (cadmium, 
chromium,  cobalt,  lead,  mercury,  nickel,  selenium,  silver,  vanadium) 
were  below  the  detection  limit. 

Copper  concentrations  exceeded  surface  water  objectives  in 
four  rivers  (Horse  River  at  Abasands,  Hangingstone  River,  Clearwater 
River  and  Surmont  Creek).  Maximum  concentrations  in  the  Hangingstone 
(43  ug/L)  and  Clearwater  (32  ug/L)  rivers  occurred  in  summer.  Maximum 
values  of  copper  in  Surmount  Creek  (48  ug/L)  occurred  in  spring. 
Maximum  values  of  copper  in  the  Horse  River  at  Abasands  (21  ug/L) 
occurred  in  winter. 

Median  values  for  manganese  and  iron  equalled  or  surpassed  the 
surface  water  objectives.  Guidelines  for  these  two  parameters  were 
established  for  aesthetic  reasons.  Both  substances,  which  are  naturally 
derived  from  river  sediments  and  weathered  bedrock,  may  impart  a taste 
to  water,  and  stain  fabrics  and  porcelain.  Manganese  may  be  taken  up  by 
plants  and  therefore  concentrations  vary  seasonally. 

7.2.3  Benthic  Invertebrates 

Tripp  and  Tsui  (1980)  reported  that  80%  of  benthic 
invertebrates  in  the  southern  tributary  streams  was  dominated  by  aquatic 
insects.  Both  species  diversity  (Shannon  1948)  and  average  number  of 
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taxa  were  greater  in  erosional  habitats.  Cleanwater  taxa  were  abundant 
in  these  habitats.  Representative  taxa  included  Ephemeroptera  (Baetis, 
Heptagenia,  Rhithrogena) , Trichoptera  (Brachycentrus , Glossosoma, 
Hydropsyche)  and  Chironomidae  (Cricotopus , Cryptochi ronomus , 
Eukiefferiel la,  Micro tend ipes , Rheotanytarsus , Thienemanniella) . 

Tripp  and  Tsui  (1980)  suggested  that  the  mean  density  of 
invertebrates  was  related  more  to  longitudinal  river  changes  than  to 
changes  in  substrate.  However,  diversity  and  average  number  of  taxa 
were  typically  higher  in  coarse  substrates  and  lower  in  muskeg  or  sand 
substrates.  Diversity  and  standing  crop  were  highest  in  summer  (August 
1978). 

Data  presented  by  Tripp  and  Tsui  (1980)  were  analysed  to 
determine  differences  in  aquatic  invertebrates  collected  in  August  1978 
from  river  sites  characterized  by  erosional  habitats: 

Stream  Number  of  sites 


Algar  1 
Cameron  1 
Christina  1 
Gregoire  1 
Hangingstone  8 
Horse  4 
Saline  2 
Saprae  3 
Surmount  2 


Emphasis  was  placed  on  Ephemeroptera  (mayflies)  because  they 
were  the  most  common  cleanwater  forms.  Because  the  mayfly  genus, 
Baeti s , occurs  in  both  erosional  and  depositional  habitats,  a correction 
factor  (mayfly  ratio)  was  included. 

A one-way  analysis  of  variance  (ANOVA)  was  used  to  determine 
differences  among  sites  for  various  biological  components  including: 

1.  % mayflies  (mayfly  density/ total  invertebrate  density 

x 100); 
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2.  mayfly  ratio  (mayfly  density  - Baetis  density/' total 
density) ; 

3.  total  density; 

4.  diversity  (Shannon  and  Weaver  1949);  and 

5.  total  number  of  taxa. 

To  satisfy  the  assumptions  of  the  ANOVA,  the  data  were 
transformed  according  to  Sokal  and  Rohlf  (1969)  as  follows:  ratios  (#1 
and  2),  arcsin  (p)5  and  total  density,  In  (x+1).  No  transformations 
were  performed  on  diversity  or  total  number  of  taxa. 

Results  indicated  that  there  were  no  differences  in  any  of  the 
components  among  the  nine  rivers  (Algar,  Cameron,  Christina,  Gregoire, 
Hangingstone , Horse,  Saline,  Saprae,  and  Surmount).  However,  if  one 
outlier  (Hangingstone  River,  site  5)  is  excluded  from  the  analysis, 
there  is  a significant  difference  (p<0.05)  in  the  number  of  taxa  among 
rivers.  Significantly  more  taxa  were  present  in  the  Hangingstone, 
Horse,  and  Christina  rivers;  fewer  taxa  were  present  in  Saline  Creek. 

7.3  RELATIONSHIPS  AMONG  PARAMETERS 

7.3.1  Water  Quality  and  Discharge 

Relationships  were  established  between  water  quality 
parameters  and  discharge  levels  for  the  Hangingstone  (CD0040)  and 
Clearwater  (CD2300)  rivers.  Results  of  the  correlation  analysis 
indicated  that,  of  all  water  quality  parameters  (median  values),  none 
was  significantly  correlated  with  conductivity  at  the  0.01  level. 
Therefore,  only  relationships  between  conductivity  and  discharge  were 
analysed  for  the  two  rivers. 

The  relationships  were  analysed  using  linear,  log/linear,  and 
log/log  transformations  of  the  data.  The  significance  of  each 
regression  was  tested  by  determining  if  the  slope  differed  significantly 
from  zero.  The  log/log  transformation  yielded  the  highest  R2  and  best 
fit  to  the  data  for  both  Hangingstone  and  Clearwater  rivers  (Table  17). 

Relationships  were  examined  to  determine  if  the  regressions  of 
conductivity  on  discharge  were  the  same  in  the  Hangingstone  and 
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Table  17.  Summary  of  regression  equations  of  conductivity  with  discharge  for 
the  Clearwater  and  Hangingstone  rivers.  Equations  with  the  best 
fit  of  the  data  were  of  the  following  log/log  format: 

Y = AxB 

where:  A = intercept,  B = slope,  Y = dependent 

variable  (conductivity),  X = independent 
variable  (discharge). 


River  df  A B R2  p 

(n-2)  (Ho:  B]  = 0) 


Clearwater  43  1821  -0.445  0.71  L0. 001 

(DC2300) 


Hangingstone  41  602  -0.574  0.71  L0.001 

(CD0040) 


L = less  than. 
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Clearwater  rivers.  The  residual  variances  between  rivers  were 
heterogenous  (F-test  for  two  samples)  (i.e.,  the  variances  of  the  two 
samples  were  not  equal  and  so  comparisons  could  not  be  made). 

7.3.2  Relationships  between  Map  Parameters  and  Biological  Components 
A multiple  regression  analysis  was  used  to  determine  the 
influence  of  land  formations  (physiology)  and  channel  features  (distance 
from  river  mouth,  slope,  width,  elevation)  on  the  benthic  invertebrate 
community  of  nine  rivers  (n=27  sites).  Five  dependent  variables  were 
examined  including  percent  mayflies,  mayfly  ratio  (mayflies  - Baetis/ 
total  density),  total  density,  diversity,  and  number  or  taxa.  Separate 
tests  were  performed  for  land  formations  and  channel  features.  Data 
used  in  both  analyses  are  presented  in  Table  18. 

The  percent  mayflies,  mayfly  ratio,  and  total  density  were 
unrelated  to  any  channel  features.  Invertebrate  diversity  increased 
with  elevation  (R2=0.20,  P <0.05)  but  was  unrelated  to  other  factors. 
The  number  of  taxa  increased  with  distance  from  river  mouth  (R2=0.15, 

P <0.05),  but  also  was  unrelated  to  other  features. 

There  were  no  relationships  between  land  formations  and  the 
percent  mayflies,  mayfly  ratio,  density,  or  the  number  of  taxa. 
Diversity,  however,  was  significantly  lower  (P<0.05)  in  sites  located 
within  the  Methy  Portage  Plain  than  in  other  regions. 

In  summary,  there  were  no  particularly  strong  relationships 
between  the  map  parameters  and  the  biological  community. 


Table  18.  Data  set  used  in  multiple  regression  analysis  relating  map  parameters 
and  biological  components.  Data  obtained  from  Tripp  and  Tsui  (1980) 
and  Turchenek  and  Lindsay  (1982). 
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8.  OVERALL  ANALYSIS  OF  WATER  QUALITY  PARAMETERS 

The  purpose  of  providing  an  overall  analysis  of  the  four 
subregions  (north,  east,  west  and  south)  is  to  use  information  on  water 
quality  parameters  to  determine  patterns  of  variability  in  river  sites 
throughout  the  AOSERP  study  area.  The  analysis  of  such  a large  data  set 
is  facilitated  by  the  use  of  muli variate  statistics.  Principal 
component  and  cluster  analyses  were  used  to  delineate  those  sites  with 
similar  water  quality  characteristics. 

Three  types  of  water  quality  parameters  were  analysed: 

1.  Routine  parameters  (including  major  inorganic  ions  and 
suspended  particulates); 

2.  Nutrients,  and 

3.  Metals. 

The  analysis  was  based  on  the  median  values  of  water  quality 
parameters  for  each  of  43  to  44  sampling  sites  in  the  study  area 
(Table  19). 

8.1  ROUTINE  PARAMETERS 

Median  values  of  18  routine  parameters  representing  44  sites 
were  analysed  using  principal  component  and  cluster  analysis.  The  first 
two  principal  components  (PCI  and  PC2)  accounted  for  68.2%  of  the 
variance  in  the  data,  of  which  the  first  component  accounted  for  48.4% 
(Table  20).  Specific  conductivity  and  silica  had  the  highest  and  lowest 
coefficients,  respectively,  for  the  first  principal  component.  Thus, 
river  sites  with  high  values  of  conductivity  and  low  values  of  silica 
have  high  scores  on  PCI.  In  contrast,  river  sites  with  low  values  of 
conductivity  and  high  levels  of  silica  have  lower  scores  on  PCI. 

Another  source  of  variability  was  the  second  principal 
component  (19.8%  of  the  variation).  Here,  non-fil terable  residue  (NFR) 
had  the  highest  positive  coefficient  and  total  alkalinity  had  the  most 
highly  significant  negative  contribution  (Table  20). 

The  relationships  among  sites  are  displayed  as  a dendrogram 
(Figure  44)  and  as  a bivariate  plot  of  factor  scores  for  each  site 
(Figure  45).  Three  major  groupings  designated  as  A,  B,  and  C emerge 
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Table  19.  Site  numbers  with  corresponding  NAQUADAT  codes  which  were  used  in  the 

cluster  and  principal  component  analysis  of  water  quality  parameters  (R, 
routine;  N,  nutrient;  M,  metal). 


Site 

R+N 

Number 

M 

NAQADAT 

Code 

Sub-region 

(N.E.W.S.) 

Site  Number 
R+N  M 

NAQUADAT 

Code  Sub-region 

1 

1 

DD0090 

N 

23 

22 

DA0110 

W 

2 

2 

DD0212 

N 

24 

23 

DA0120 

U 

3 

3 

KF0100 

N 

25 

24 

DA0121 

W 

4 

4 

KF0101 

N 

26 

25 

DA0130 

W 

5 

5 

KF0140 

N 

27 

26 

DA0140 

W 

6 

6 

NA0030 

N 

28 

27 

DA0150 

W 

7 

7 

KF0200 

N 

29 

28 

DA0160 

w 

8 

8 

KF0201 

N 

30 

29 

DA0170 

w 

9 

9 

DD0080 

N 

31 

30 

DA0180 

w 

10 

10 

MD2000 

N 

32 

31 

DA0181 

w 

11 

11 

DA0060 

E 

33 

32 

DA0182 

w 

12 

12 

DA0080 

E 

34 

33 

DA0183 

w 

13 

13 

DA0081 

E 

35 

34 

DA0190 

w 

14 

14 

DA0083 

E 

36 

35 

DB0011 

w 

15 

15 

DA0084 

E 

37 

36 

DB0020 

w 

16 

16 

DA0085 

E 

38 

37 

DB0030 

w 

17 

17 

DA0090 

E 

39 

38 

DB0040 

w 

18 

18 

DC0010 

E 

40 

39 

CC0300 

' s 

19 

19 

DD0020 

E 

41 

40 

CD0040 

s 

20 

20 

DA0070 

W 

42 

41 

CD2300 

s 

21 

- 

DA0098 

W 

43 

42 

CE2000 

s 

22 

21 

DA0100 

W 

44 

43 

CE2121 

s 
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Table  20.  Principal  component  coefficient  (loadings)  for  the  first  two 
components  (PCI  and  PC2):  routine  parameters. 


Parameter 

PCI 

PC2 

pH 

.183 

.034 

Total  alkalinity 

.252 

-.332 

Hardness 

.315 

-.071 

Conductivity 

.329 

.073 

Turbidity 

.052 

.368 

Silica 

.039 

-.166 

Calcium 

.305 

-.066 

Magnesium 

.320 

.095 

Sodium 

.266 

.110 

Potassium 

.266 

.172 

FI uoride 

.224 

.199 

Bicarbonate 

.253 

-.331 

Chloride 

.151 

.233 

Sulphate 

.196 

.365 

Total  organic  carbon 

.165 

-.188 

Total  inorganic  carbon 

.192 

-.325 

Filterable  residue 

.327 

.075 

Non-fil terable  residue 

.093 

.421 

Eigenvalues 

8.71 

3.56 

% total  variance 

48.42 

19.79 

Cumulative  % t.v. 

48.42 

68.21 
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Figure  44.  Cluster  analysis  of  routine  parameters 

at  sample  sites  within  the  four  subregions  . 
Similarity  is  measured  by  error  sum  of 
squares  . 


PC2 
( NFR  ) 


Figure  45.  Principal  component  ordination  of  sampling 
sites  using  routine  parameters  from  the 
four  subregions  . Dashed  lines  enclose  sub- 
groups of  river  sites  . Chemical  labels  on 
axes  are  from  Table  20  . 
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from  the  cluster  analysis  of  the  routine  parameters.  Sites  that  have 
similar  values  of  water  quality  parameters  will  form  a cluster  on  the 
bivariate  plot.  Sites  within  clusters  have  been  delineated  on  a map  of 
the  study  area  (Figure  46). 

The  mean  ± standard  error  of  the  routine  parameters  with  the 
highest  and  lowest  coefficients  (i.e.,  conductivity,  silica,  NFR,  and 
alkalinity)  for  PCI  and  PC2  are  presented  in  Table  21.  The  range  of 
values  for  these  parameters  encountered  within  the  AOSERP  study  area  are 
typical  of  lotic  conditions  throughout  Alberta. 

Sites  within  group  A have  characteristically  low  values  of 
conductivity,  silica,  and  alkalinity,  with  moderate  levels  of  suspended 
solids  and  turbidity  (Table  21).  Suspended  sediments  and  turbidity 
reduce  light  levels  available  for  algal  productivity.  Representative 
group  A sites  included  those  in  the  northeast  (Lake  Athabasca  Delta), 
west  (Ells  River,  Thickwood  Creek),  and  south  (Horse,  Clearwater,  and 
Surmount  rivers)  regions.  Many  of  these  sites  occurred  on  archean 
(granite)  or  devonian  bedrock  and  are  less  productive  than  rivers  which 
flow  over  cretaceous  rock. 

Group  B sites  had  the  lowest  NFR  and  highest  silica  values; 
conductivity  levels  were  moderate  (Table  21).  Substrates  at  these  sites 
receive  maximum  available  light  for  algal  growth  and  abundant  levels  of 
silica  for  diatom  production.  These  sites,  therefore,  should  represent 
areas  of  maximum  productivity  within  the  region.  Representative 
waterways  of  the  eastern  region  include  Muskeg,  Steepbank,  Firebag,  and 
Tar  rivers.  The  Hangingstone  and  MacKay  mainstem  also  are  included 
within  this  group. 

Group  C sites  are  categorized  as  having  maximum  values  of 
conductivity,  NFR,  and  alkalinity;  silica  concentrations  are  moderate 
(Table  21).  These  river  sites  are  often  highly  turbid  and  maintain  a 
high  sediment  load,  factors  which  reduce  light  available  for 
photosynthesis.  Furthermore,  substrate  interstices  are  filled  with 
sediment  which  reduce  the  effective  habitat  available  to  benthic 
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Figure  46. 


Site  groupings  based  on  principal  component 
analysis  using  median  values  of  routine 
parameters  . 
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Table  21.  Mean  (±  standard  error)  of  routine  parameters  indicated  by  the 

highest  and  lowest  loading  coefficients  for  the  first  two  principal 
components.  Groups  are  indicated  in  Figures  46  and  47. 


Parameters 


Group 

(No.  of  Sites) 

Conductivity 

(uS/cm) 

Silica 

(mg/L) 

NFR 
(mg/L ) 

Total  Alkalinity 
(mg/L) 

A (n  = 15) 

169.2  ± 15.76 

5.6  ± 0.83  a 
(4.9  ± 0.49)a 

14.5  ± 3.16 

73.0  ± 6.41 

B (n  = 15) 

276.7  ± 14.43 

8.6  ± 0.71 

8.9  ± 1.50 

143.9  ± 7.70 

C (n  = 15) 

449.9  ± 15.94 

6.9  ± 0.74 

20.9  ± 3.91 

142.6  ± 14.69 

dSite  19  (Richardson  River)  detailed  from  calculation. 
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invertebrates.  Representative  sites  are  lakes  and  channels  of  the 
northwest  region  and  river  sites  draining  the  east  slopes  of  the  Birch 
Mountains. 

8.2  NUTRIENTS 

Clustering  techniques  were  applied  to  44  sites  using  median 
values  of  nutrient  parameters,  including  total  phosphorus  (TP), 
orthophosphate,  nitrate/nitrite  nitrogen,  total  kjeldahl  nitrogen  (TKN) , 
and  ammonia.  The  cumulative  per  cent  variability  removed  by  the  first 
two  principal  components  was  64.0%,  representing  PCI  (34.4%)  and  PC2 
(29.6%)  (Table  22). 

Total  phosphorus  and  nitrate/nitrite  nitrogen  had  the  highest 
and  lowest  coefficients,  respectively,  for  the  first  principal 
component.  Thus,  sites  with  characteristically  high  concentrations  of 
TP  and  low  concentrations  of  nitrate/nitrite  nitrogen  have  high  factor 
scores  on  PCI.  For  the  second  principal  component,  nitrate/nitrite 
nitrogen  had  the  highest  positive  coefficient  and  TKN  had  the  most 
negative  coefficient.  Sites  with  high  nitrate/nitrite  nitrogen  values 
have  high  factor  scores  on  PC2. 

The  relationships  among  sites,  each  of  which  is  represented  by 
a factor  score,  are  presented  as  a bivariate  plot  of  PCI  and  PC2 
(Figure  48).  Sites  with  similar  values  of  TP,  TKN,  and  nitrate/nitrites 
form  clusters  on  the  bivariate  plot.  Most  sites  were  clustered  at  the 
intersection  of  the  two  components.  High  median  values  of 
nitrate/nitrite  nitrogen  were  noted  for  the  Horse  River  (site  40, 

0.2  mg/L)  and  Surmount  Creek  (site  43,  0.32  mg/L).  Thus,  these  two 
sites  are  distinct  from  other  sites  on  the  bivariate  plot. 

Results  from  a cluster  analysis  of  the  data  are  displayed  as  a 
dendrogram  in  Figure  47.  Again,  sites  40  and  43  are  distinct.  Site  30 
(lower  Ells  River)  also  is  isolated  from  the  other  sites.  The  median 
value  of  ammonia  at  site  30  was  0.35  mg/L;  the  range  of  ammonia 
concentrations  for  the  other  43  sites  was  0.03  to  0.19  mg/L.  Median 
values  of  other  nutrients  recorded  for  site  30  were  within  the  range  of 
the  concentrations  of  the  other  sampling  sites. 
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Table  22.  Principal  component  coefficients  (loadings)  for  the  first  two 
components  (PCI  and  PC2):  nutrient  parameters. 


Parameter 

PCI 

PC2 

Total  phosphorus 

.610 

-.083 

Orthophosphate 

.501 

.514 

Nitrate/nitrite 

.100 

.691 

Total  k j el dah 1 nitrogen 

.579 

-.375 

Ammonia 

.178  . 

-.333 

Eigenvalues 

1.72 

1.48 

% of  total  variance 

34.44 

29.60 

Cumulative  l t.v. 

34.44 

64.04 
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Figure  47.  Cluster  analysis  of  nutrients  at  sample 
sites  within  the  four  subregions  . 
Similarity  is  measured  by  error  sum 
of  squares  . 


PC  2 

( NO3/NO*  ) 


Figure  48.  Principal  component  ordination  of  sampling 
sites  using  nutrient  parameters  from  the 
four  subregions  . Chemical  labels  on  axes 
are  from  Table  22  . 
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The  mean  ± standard  error  of  the  nutrients  with  highest  and 
lowest  coefficients  (i.e.,  TP,  nitrate/nitrite,  and  TKN)  for  PCI  and  PC2 
are  presented  for  each  of  the  three  major  groups  (A,  B,  and  C) 
(Table  23).  These  groups  are  also  delineated  on  a map  of  the  study  area 
(Figure  49). 

TKN  is  a measure  of  organic  nitrogen  and  ammonia,  both  of 
which  are  important  in  biological  activity.  Inorganic  forms  of  nitrogen 
include  ammonia,  nitrates,  and  nitrites.  With  adequate  oxygen,  ammonia 
is  converted  to  nitrate  by  a nitrification  process.  Nitrite  is  rapidly 
converted  to  nitrate  in  the  presence  of  oxygen. 

Total  phosphorus  (inorganic  and  organic  phophorus)  is  a 
measure  of  both  dissolved  and  particulate  forms.  The  dissolved  form 
which  is  biologically  available  can  be  a limiting  factor  for  plant 
growth  (McNeely  et  al.  1979).  Orthophosphates  represent  the  dissolved 
form  of  phosphorus.  Although  median  values  of  TP  increased 
substantially  from  group  B (62.2  ug/L  ± 3.23)  to  group  C (117.7  ug/L 
± 9.69),  relatively  minor  increases  occurred  in  orthophosphate  (group  B, 
15.7  ± 1.80;  group  C,  19.1  ± 2.45). 

Group  A sites  were  characterized  by  relatively  low 
concentrations  of  TP  and  TKN  compared  to  other  AOSERP  sites. 
Representative  sites  include  the  Athabasca  Delta  in  the  northeast,  the 
Ells,  Upper  Tar,  and  Thickwood  Creeks  in  the  west  and  Firebag,  Steepbank 
and  Clearwater  rivers  in  the  east. 

Group  B sites  had  moderate  concentrations  of  TP  and  high 
levels  of  nitrate/nitrite  nitrogen  and  TKN  (Table  23)  which  suggests 
that  representative  streams  are  productive.  Sites  in  this  category 
include  lakes  and  channels  of  the  northwest,  the  Muskeg  watershed 
(excluding  Stanley  Creek),  Hangingstone  River,  Poplar  Creek,  lower 
MacKay  River,  and  Joslyn  Creek. 

Group  C sites  reflected  conditions  with  the  highest  TP  and  TKN 
values;  nitrate/nitrite  nitrogen  values  are  lowest  (Table  23).  Stanley 
Creek  in  the  eastern  region  and  several  tributaries  in  the  west 
(Asphalt,  Calumet,  lower  Tar,  Dunkirk,  and  Beaver  rivers)  are 
representative  sites. 
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Table  23.  Mean  (±  standard  error)  of  nutrients  indicated  by  the  highest  and 
lowest  loading  coefficients  for  the  first  two  principal  components. 
Groups  are  indicated  in  Figures  48,  49,  and  50. 


Parameters 

Group 

(No.  of  Sites) 

TP 

(ug/L) 

no3no, 
( ul/l) 

TKN 

(ug/L) 

A (n  = 17) 

46.2  ± 3.61 

34.9  ± 5.35 

804.7  ± 

62.59 

B (n  = 15) 

62.2  ± 3.23 

36.8  ± 5.61 

1179.3  ± 

48.26 

C (n  = 9) 


117.8  ± 9.69 


25.6  ± 6.04 


1173.3  ± 168.90 


mnn 
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Figure  49 


Site  groupings  based  on  principal  component 
analysis  using  median  values  of  nutrients  . 


158 


Alberta  Surface  Water  Objectives  for  total  phosphorus  and 
nitrogen  are  50  ug/L  and  1.0  mg/L,  respectively.  No  surface  water 
objectives  are  available  for  nitrogen  components.  The  mean  of  all 
median  values  of  TP  within  clusters  B and  C exceeded  the  surface  water 
objectives.  The  elevated  values  of  TP,  however,  are  not  necessarily 
related  to  increased  eutrophication.  Rather,  the  apatite  (biologically 
unavailable)  form  of  phosphorus  increases  with  runoff.  Thus,  rivers 
draining  the  east  slope  of  the  Birch  Mountains  which  have  high  TP  values 
will  have  correspondingly  high  levels  of  suspended  sediment. 

8.3  METALS 

Cluster  and  principal  component  techniques  were  applied  to  43 
sites  using  median  values  of  seven  metals  including  aluminum,  arsenic, 
boron,  copper,  iron,  manganese,  and  zinc.  One  site  (DA0098,  mouth  of 
the  Ells  River)  was  deleted  because  of  insufficient  data.  The  seven 
metals,  which  were  selected  from  the  full  suite  typically  analysed, 
exhibited  values  larger  than  the  laboratory  detection  limit. 

The  first  (46.6%)  and  second  (20.9%)  principal  components 
accounted  for  a cumulative  per  cent  variability  of  67.5 % (Table  24). 
Zinc  (along  with  iron)  and  boron  had  the  highest  and  lowest 
coefficients,  respectively,  for  PCI.  Boron  and  aluminum  had  the  highest 
positive  and  negative  coefficients,  respectively,  for  PC2. 

The  relationships  among  sites  were  displayed  as  a bivariate 
plot  of  factor  scores  of  each  site  (Figure  51).  The  dendrogram 
illustrates  relationships  among  sites  indicated  by  the  cluster  analysis 
(Figure  50).  Sites  within  clusters  have  been  delineated  on  a map  of  the 
study  area  (Figure  52). 

The  bivariate  plot  and  dendrogram  reveals  that  there  are  three 
major  groups  (A,  B,  and  C).  Group  A consists  of  21  similar  sites. 
Group  C consists  of  16  sites,  14  of  which  occur  together  and  two  (site 
27,  the  lower  Tar;  site  28,  Joslyn  Creek)  which  are  distinct.  Group  B 
consists  of  two  subgroups,  both  quite  distinct  from  either  groups  A 
or  C. 
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Table  24.  Principal  component  coefficients  (loadings)  for  the  first  two 
components  (PCI  and  PC2):  metal  parameters. 


Parameter 

PCI 

PC2 

Aluminum 

.426 

-.381 

Arsenic 

.257 

.328 

Boron 

.183 

.576 

Copper 

.399 

-.363 

Iron 

.463 

.335 

Manganese 

.344 

.294 

Zinc 

.476 

-.292 

Eigenvalue 

3.26 

1.47 

% total  variance 

46.57 

20.95 

Cumulative  % t.v. 

46.57 

67.53 
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Figure  50.  Cluster  analysis  of  metals  at  sample 
sites  within  the  four  subregions  . 
Similarity  is  measured  by  error  sum 
of  squares 
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Figure  51.  Principal  component  ordination  of  sampling 
sites  using  metal  parameters  from  the 
four  subregions  . Chemical  labels  on  axes 
are  from  Table  24  . 
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Figure  52 


Site  groupings  based  on  principal  component 
analysis  using  median  values  of  metals  . 
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The  mean  ± standard  error  of  the  metals  (zinc,  iron,  boron, 
and  aluminum)  with  the  highest  and  lowest  coefficients  for  the  first  two 
principal  components  are  presented  in  Table  25.  Values  for  zinc  and 
boron  are  within  the  Alberta  Surface  Water  Guidelines  of  50  ug/L  and 
500  ug/L,  respectively.  No  guidelines  are  available  for  aluminum.  The 
iron  values  exceed  the  surface  water  objectives  of  300  ug/L  for  all 
groups  (A,  B,  and  C).  There  are,  however,  few  adverse  effects  of  iron 
to  man  (McNeely  et  al . 1979).  Elevated  levels  of  iron  in  surface  waters 
typically  result  from  the  weathering  of  rocks. 

Group  A sites  exhibited  low  levels  of  zinc,  iron,  and  boron 
with  moderate  levels  of  aluminum  (Table  25).  Representative  sites 
include  the  Athabasca  Delta  of  the  northeast,  as  well  as  the  Firebag, 
Muskeg  (excluding  Kearl  Lake  tributary),  Steepbank,  and  Clearwater 
rivers.  Other  sites  include  the  Ells  River,  upper  Tar,  Thickwood,  and 
Poplar  creeks  in  the  west,  and  Surmount  Creek  in  the  south. 

Group  B sites  exhibited  the  highest  levels  of  zinc  and  iron 
compared  to  group  A and  C (Table  25).  Levels  of  boron  were  moderate, 
whereas  concentrations  of  aluminum  were  relatively  low.  Sites  23 
(Asphalt  Creek)  and  24  (Eymundson  Creek)  within  group  B had  the  two 
highest  values  for  iron  (site  23,  4230  ug/L;  site  24,  4200  ug/L)  and 
zinc  (site  23,  30  ug/L;  site  24,  22  ug/L)  in  the  study  area.  These  two 
sites,  which  represent  the  Asphalt  watershed,  form  a distinct  subgroup 
within  group  B.  Other  sites  within  this  group  include  the  lakes  and 
channels  of  the  northwest. 

Group  C sites  had  the  highest  levels  of  boron,  but  moderate 
levels  of  zinc,  iron,  and  aluminum  compared  to  the  other  groups 
(Table  25).  Representative  sites  occur  predominately  in  the  west 
(Beaver,  MacKay,  Dunkirk,  and  Dover  rivers  and  streams  draining  the  east 
slopes  of  the  Birch  Mountains)  and  south  (the  Horse,  Hangingstone , and 
Christina  rivers).  The  Kearl  Lake  tributary  is  the  only  stream  from  the 
eastern  region  within  group  C.  Elevated  concentrations  of  boron  were 
noted  frequently  in  winter  when  they  are  associated  with  groundwater 
inputs.  Many  of  the  sites  in  group  C correspond  to  groupings  B and  C of 


163 

Table  25.  Mean  (±  standard  error)  values  of  metals  for  sites  within  groups 
determined  by  principal  component  analysis.  The  metals  listed 
represented  the  highest  and  lowest  loading  coefficients  for  the 
first  two  principal  components.  Groups  are  indicated  in 
Figures  50,  51 , and  52. 


Group 

Zn 

(ug/L) 

Fe 

(ug/L) 

Bo 

(ug/L) 

A1 

(ug/L) 

A (n  = 

21) 

6.8  ± 0.49 

827.9  ± 68.29 

88.1  ± 10.52 

129.2  ± 

20.74 

B (n  = 

6) 

21.7  ± 2.88 

2192.5  ± 653.51 

116.7  ± 22.46 

600.8  ± 

118.56 

B'(n  = 

4) 

19.5  ± 3.66 

1181.2  ± 212.25 

82.5  ± 7.50 

575.0  ± 

42.87 

C (n  = 

16) 

9.5  ± 0.93 

1619.7  ± 210.12 

192.5  ± 17.69 

157.5  ± 

21.48 
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the  routine  parameters  ( i . e . , sites  which  exhibited  elevated  saline 
inputs  from  groundwater  sources) . 
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10. 

10.1 


APPENDICES 


LIST  OF  CHEMICAL  PARAMETERS  ANALYSED  AND  THE  ASSOCIATED  CODE 
FOR  WHICH  STANDARD  METHODS  ARE  DESCRIBED  (ENVIRONMENT  CANADA 
1974)  (TABLE  26) 
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Table  26. 

List  of  chemical  parameters  analysed  and  the  associated  code 
for  which  standard  methods  are  described  (Environment  Canada 
1974).  Unless  otherwise  specified,  all  units  are  in  mg/L. 

Parameter 

Code 

pH 

10301L 

Total  alkalinity  as  CaCO^ 

10101L 

Phenol phthalein  alkalinity  as  CaCO^ 

10151L 

Total  hardness  as  CaCO^ 

10602L 

Specific  conductivity  (uS/cm) 

02041L 

Turbidity  (JTU) 

02073L 

True  colour  (relative  units) 

02021L 

Si  1 ica 

14101 L 

Calcium 

20103L 

Magnesium 

12102L 

Sodium 

11102L 

Potassium 

19102L 

Bicarbonate 

06201L 

Carbonate 

06301L 

Chloride 

17203L 

Sulphate 

16306L 

Dissolved  oxygen 

08101L 

Total  organic  carbon 

06001L 

Filterable  residue 

10451L 

Nonfilterable  residue 

10401L 

Total  inorganic  carbon 

06051L 

Phenolic  material 

06532L 

Oil  and  grease  (petroleum  ether  extraction) 

06521L 

Oil  and  grease  (freon  extraction) 

06524L 

Total  phosphorus 

15406L 

Orthophosphate  phosphorus 

15256L 

Total  k j el dah 1 nitrogen 

07013L 

continued  . . 
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Table  26.  Concluded. 


Parameter 

Code 

Ammonia 

07555L 

Nitrate  and  nitrite  nitrogen 

07110L 

Aluminum,  extractable 

13302L 

Arsenic,  dissolved 

33104L 

Boron,  dissolved 

05105L 

Cadmium,  extractable 

48302L 

Chromium  hexavalent 

24 10 1L 

Cobalt,  extractable 

27302L 

Copper,  extractable 

29305L 

Iron,  extractable 

26304L 

Lead,  extractable 

82302L 

Manganese,  extractable 

25304L 

Mercury,  total 

80015L 

Nickel,  extractable 

28302L 

Selenium,  dissolved 

34102L 

Silver,  extractable 

47302L 

Vanadium,  total 

23002L 

Zinc,  extractable 

30305L 
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10.2  SEASONAL  SUMMARY  OF  ROUTINE  WATER  QUALITY  PARAMETERS 

(TABLES  27  TO  38),  NUTRIENTS  (TABLES  39  TO  44),  COLIFORMS 
(TABLES  45  TO  46),  AND  METALS  ( TABLES  47  TO  58),  FOR  12  SITES 
ALONG  THE  ATHABASCA  RIVER. 


Table  27.  Seasonal  summary  of  routine  water  quality  parameters  for  site  CC0012  (centre)  on  the  Athabasca 
River  (1976  to  1983).  Winter  (November  1 to  April  15),  spring  (April  16  to  June  30), 
summer  (July  1 to  August  31),  fall  (September  1 to  October  31). 
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Table  29.  Seasonal  summary  of  routine  water  quality  parameters  for  site  DA0203  (right  bank)  on  the 
Athabasca  River  (1976  to  1983).  Winter  (November  1 to  April  15),  spring  (April  16 
to  June  30),  summer  (July  1 to  August  31),  fall  (September  1 to  October  31). 
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Table  30.  Seasonal  summary  of  routine  water  quality  parameters  for  site  DA0205  (1976  to  1983). 
Winter  (November  1 to  April  15),  spring  (April  16  to  June  30),  summer  (July  1 to 
August  31),  fall  (September  1 to  October  31). 
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Table  31.  Seasonal  summary  of  routine  water  quality  parameters  for  site  DA0205  (right  bank)  on  the 

Athabasca  River  (1976  to  1983).  Winter  (November  1 to  April  15),  spring  (April  16  to  June  30), 
summer  (July  1 to  August  31),  fall  (September  1 to  October  31). 
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Table  32.  Seasonal  summary  of  routine  water  quality  parameters  for  site  DA0206  (left  bank)  on  the 
Athabasca  River  (1976  to  1983).  Winter  (November  1 to  April  15),  spring  (April  16  to 
June  30),  summer  (July  1 to  August  31),  fall  (September  1 to  October  31). 
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Table  33.  Seasonal  summary  of  routine  water  quality  parameters  for  site  DA0206  (right  bank)  on  the 
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Table  34.  Seasonal  summary  of  routine  water  quality  parameters  for  site  DA0207  (left  bank)  on  the 
Athabasca  River  (1976  to  1983).  Winter  (November  1 to  April  15),  spring  (April  16  to 
June  30),  summer  (July  1 to  August  31),  fall  (September  1 to  October  31). 
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Table  35.  Seasonal  summary  of  routine  water  quality  parameters  for  site  DA0207  (right  bank)  on  the 
Athabasca  River  (1976  to  1983).  Winter  (November  1 to  April  15),  spring  (April  16  to 
June  30),  summer  (July  1 to  August  31),  fall  (September  1 to  October  31). 
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Table  37.  Seasonal  summary  of  routine  water  quality  parameters  for  site  DA0208  (right  bank)  on  the 
Athabasca  River  (1976  to  1983).  Winter  (November  1 to  April  15) s spring  (April  16  to 
June  30),  summer  (July  1 to  August  31),  fall  (September  1 to  October  31). 


196 


E 

X 

o o o o 

X 

rv  o «-! 

O X 

StOCDO 

tO 

•-I  03  CO  CD 

It3 

JJ  to 

>1  E 

in  cm  cm  ro 

E 

CM  * — 1 r-l  • — t 

t-  E 

to 

C_J 

> 

E 

3 - — 

o — 

E 

•<-  —1 

•r-  _l 

O U "3 

CO  in  00  ro 

t/1  -v  TD 

O •— 1 "O’  IX 

E 'v'O 

CO  Ch  O to 

3 "'v  CD 

CM  C CM  00 

to  oi  at 

to  31  QJ 

"OWE 

*3-  CM  CM  CM 

to  E E 

t-H  —1  O O 

31  E E 

E 3 



O « — • 

o 

O 

C_> 

CL 

IO 

e 

imncDcn 

E 

in  to  co  to 

■U  E 

O 
t— 

in  to  co  in 

X 

CO  CM  LO  00 

X 

inoooi 

X 

CM  PI  V ft 

to 

00  O O to 

to 

ro  >— t •— i t— t 

to 

tt  CMHP1 

E 

_ ^ 

E 

E 

t/i 

at 

i/i  - — ' 

+->  -v 

CD  _J 

E —l 

IO  __l 

e "v-o 

NOIMCO 

3 -v-3 

voajtt 

o 

•3  CD  QJ 

inocoo 

•r-  CO  QJ 

CM  * — < r— 1 

cl  3>  at 

ro  — t— i cm 

o 

s-  E E 

»— i , — i . — i 

“3  E E 

• — EE 

in 

fO  v— • 

O'— 

S'— 

DC 

1/1 

in 

e 

to  to  CO  CO 

E 

to  to  co  to 

E 

to  ^3-  ro  xo- 

X 

CMOUOOO 

X 

IO  OO  Ol  CM 

X 

ro  O co  co 

>>  n 

NOiOfD 

to 

to 

ro  — — 

4->  E 

r— 1 1 — 1 * — < 

E 

E 

E 

E 

CD 

3 ' 

X3  ■ — • 

tO  _J 

•r-  _J 

•t-  _J 

-x  \ xi 

<tCMCO>-t 

toxxj 

ntvcoat 

S-  -X.-3 

<— t rv  co  «-t 

o 

r-  CM) 

■>o-  01  cn  o 

CD  31  QJ 

O 31  CD 

ro  r-t 

in 

CEE 

r— t 

E E E 

' — EE 

CM 



-E  v-’ 

to 

c_> 

<0 

s: 

+-> 

o 

(—  c 

to  to  co  to 

E 

to  in  co  to 

E 

to  in  co  in 

X 

rv  01  o O 

to  to 

X 

Ol  31  Ol 

X 

cr>  ro  co  cm 

JO 

to 

<3-CMCCHt 

to 

O r-t  CM  co 

E 

tv  tv  oo  CO 

00  00 

E 

E 

ai 

4-> 

CM  HHH 

o o 

to 

E *-v 

E ' — ' 

3 —I 

O —1 

E 

'tcorvct- 

to  to 

-r-  ~'v'3 

CM  00  CO  O 

JJ  \-3 

in  cm  in  in 

U (Ml 

*3-  CM  CM  CO 

S-  31  CD 

IV.  r-t  r-t  CM 

E 

tv  tv  rv.  tv. 

CO  CO 

«—  E E 

to  E E 

t— 1 1 — t • — t ?— i 

to  - — 

U — - 

C_J 

JO 

CO 

CL 

T3 

to  tv  CO  CO 

OJ 

E 

to  to  CO  to 

E 

in  "O"  m 

E 

IVIVtVIV 

E 

IO  *3-  «d-  «3- 

1-  3>  J- 

S-  31  S- 

S-  31  s- 

QJ  E CD 

CD  E 3 

CD  E CL) 

4->  — E — 

O CD 

+->  -r-  £ — 

O CD 

-t-  E r— 

O CD 

E S-  E 1 — 

3 3 

E S-  I — 

3 3 

E J-  Ir. 

2 2 

a 3 id 

CO  o 

■i-  a 3 to 

in  a 

•t-  Q.  3 to 

in  o 

2 W 1/1  U. 

c£  <_) 

2 WWIjl 

«=C  cj 

2 in  in  ll 

<C  o 

continued 


197 


CO  C vo  o 


fo  oo  co  n 
o o o o 


i-unuirx 
— O O O 
o o o o 

o o'  o o" 


VO  CO  CM  C\J  LD 

o o o o o 

o o o o o 

o o o o o 


in  LD  VO  IT) 


VO  CO  O CM 
VO  CM  C CO 
CSJ  CO 


in  in  vd  in 


X 

E 


o *3-  in  co 

OlVOlOr-l 

PJ-IrHM 


"O 

OJ 

~o 

13 

o 

c 

o 

o 


CU  - 
O'  _J 

'-»»  ~o 
Ol  DMU 
— EE 


O C C co 

VO  C LO  VO 
CVJ  »“•  1— I f"H 


c 


in  in  in  in 


ro 


-Q 

f0 


S-  cn  i- 

QJ  C (D 

E — O O 

c i-  E — 3 3 

•I-  a 3 (O  WD 

3 i/o  oo  U-  co 


Table  38.  Seasonal  summary  of  routine  water  quality  parameters  for  site  DD0010  (centre)  on  the  Athabasca 
River  (1976  to  1983).  Winter  (November  1 to  April  15),  spring  (April  16  to  June  30), 
summer  (July  1 to  August  31),  fall  (September  1 to  October  31). 
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Table  40.  Seasonal  summary  of  nutrients  for  sites  DA0203  (left  and  right  banks)  on  the  Athabasca  River 


201 


>4 

C5 

r"0 


en-o 


o o o o 


CM  o o o 

o o o o 


o o o o 


NOOO 

o o’  o*  o’ 


o 

CO 


«3-  CTi  O OO 


O O O O 


O O O O 


O O O O 


o o o o 


s- 

CL 

c 


<toioco 


CD 


s- 

Q. 

CO 


PO 

s- 

s-  a) 

CL  jD 

c o 

4-> 

o u 
+->  o 


<D  i— 
_D 

E S- 
CU  CU 
> .Q 
O E 
S QJ 
4-> 
CL 
5-  CU 
QJ  CO 
+->  ' 
SC 


■st^MO 


•JOiON 


^tOiOCO 


fO 


ro^ 
CO  I— 
cr>  ro 


o co 

-P  C5 
CD 
CD  3 
C 
CT> 

r-  O 
— ' -P 


^-o 
cn  a> 
3 E 


O O O r- 
cn  ro 
rH  mo 


CM  O O 

'd-  00  ro  oo 
ro  to 


ri-aiOCO 


~o 

CD 


O QJ 
<C  r— 

O • 


aj  jz. 


O 4— 


J- 

CL3  (O 
l/H/Ui. 


•r-  E r-  O O 


less  than. 


Table  41.  Seasonal  summary  of  nutrients  for  sites  DA0205  (left  and  right  banks)  on  the  Athabasca  River 
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Table  42.  Seasonal  summary  of  nutrients  for  sites  DA0206  (left  and  right  banks)  on  the  Athabasca  River 
(1976  to  1983).  Winter  (November  1 to  April  15),  spring  (April  16  to  June  30),  summer  (July  1 
to  August  31),  fall  (September  1 to  October  31). 
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Table  43.  Seasonal  summary  of  nutrients  for  sites  DA0207  (left 
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Table  44.  Seasonal  summary  of  nutrients  for  sites  DA0208  (left  and  right  banks)  on  the  Athabasca  River 
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median  value  = minimum  value. 


Table  45.  Summary  of  total  coliforms  (NAQUADAT  parameter  code:  36001L)  for  12  sites 
along  the  Athabasca  mainstem.  Units  are  recorded  as  most  probable  number 
(MPN)  per  100  mL.  Winter  (November  1 to  April  15),  spring  (April  16  to 
June  30),  summer  (July  1 to  August  31),  fall  (September  1 to  October  31). 
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If  value  = 320  or  MPN/mL  is  greater  than  320  or  2400,  respecti vely. 

No  data,  but  for  purpose  of  analysis  DA0208  (right  bank)  values  were  assigned  to 
DD0010  (centre). 


Table  46.  Summary  of  fecal  conforms  (NAQUADAT  parameter  code:  3601 1L)  for  12  sites 
along  the  Athabasca  mainstem.  Units  are  recorded  as  most  probable  number 
(MPN)  per  100  mL.  Winter  (November  1 to  April  15),  spring  (April  16  to 
June  30),  summer  (July  1 to  August  31),  falk  (September  1 to  October  31). 


207 


i 
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No  data,  but  for  purpose  of  analysis  DA0208  (right  bank)  values  were  assigned  to 
DD0010  (centre). 


Table  47.  Seasonal  summary  of  metals  for  site  CC0012  (centre)  on  the  Athabasca 

River  (1976  to  1983).  Winter  (November  1 to  April  15),  spring  (April  16 
to  June  30),  summer  (July  1 to  August  31),  fall  (September  1 to  October  31). 
All  values  in  ug/L. 
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Extractable. 

Dissol ved. 

ASWO  = Alberta  Surface  Water  Objectives,  CDWG  = Canadian  Drinking  Water  Guidelines. 
If  n = 2,  median  value  = minimum  value. 


Table  48.  Seasonal  summary  of  metals  for  site  DA0203  (left  bank)  on  the  Athabasca 
River  (1976  to  1983).  Winter  (November  1 to  April  15),  spring  (April  16 
to  June  30),  summer  (July  1 to  August  31),  fall  (September  1 to  October  31). 
All  values  in  ug/L. 
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median  value  = minimum  value. 


Table  49.  Seasonal  summary  of  metals  for  site  DA0203  (right  bank)  on  the  Athabasca 
River  (1976  to  1983).  Winter  (November  1 to  April  15),  spring  (April  16 
to  June  30),  summer  (July  1 to  August  31),  fall  (September  1 to  October  31). 
All  values  in  ug/L. 
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Table  51.  Seasonal  summary  of  metals  for  site  DA0205  (right  bank)  on  the  Athabasca 
River  (1976  to  1983).  Winter  (November  1 to  April  15) s spring  (April  16 
to  June  30),  summer  (July  1 to  August  31),  fall  (September  1 to  October  31). 
All  values  in  ug/L. 
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Table  52.  Seasonal  summary  of  metals  for  site  DA0206  (left  bank)  on  the  Athabasca 
River  (1976  to  1983).  Winter  (November  1 to  April  15),  spring  (April  16 
to  June  30),  summer  (July  1 to  August  31),  fall  (September  1 to  October  31). 
All  values  in  ug/L. 
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median  value  = minimum  value. 


Table  53.  Seasonal  summary  of  metals  for  site  DA0206  (right  bank)  on  the  Athabasca 
River  (1976  to  1983).  Winter  (November  1 to  April  15),  spring  (April  16 
to  June  30),  summer  (July  1 to  August  31),  fall  (September  1 to  October  31). 
All  values  in  ug/L. 
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median  value  = minimum  value. 


Table  54.  Seasonal  summary  of  metals  for  site  DA0207  (left  bank)  on  the  Athabasca 
River  (1976  to  1983).  Winter  (November  1 to  April  15),  spring  (April  16 
to  June  30),  summer  (July  1 to  August  31),  fall  (September  1 to  October  31). 
All  values  in  ug/L. 
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Table  55.  Seasonal  summary  of  metals  for  site  DA0207  (right  bank)  on  the  Athabasca 
River  (1976  to  1983).  Winter  (November  1 to  April  15),  spring  (April  16 
to  June  30) , summer  (July  1 to  August  31),  fall  (September  1 to  October  31). 
All  values  in  ug/L. 
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Table  56.  Seasonal  summary  of  metals  for  site  DA0208  (left  bank)  on  the 

Athabasca  River  (1976  to  1983).  Winter  (November  1 to  April  15), 
spring  (April  16  to  June  30),  summer  (July  1 to  August  31), 
fall  (September  1 to  October  31).  All  values  in  ug/L. 


217 


X 

o © 

x 

HOinn 

X 

m to  cn  id 

03 

O CM 

03 

CMOCMlf) 

03 

E 

E 

CM  CM 

E 

l/l 

4-> 

X 

■M 

o 

LU 

X 

LU 

© o 

T3 

o o o 

-O 

O CO  o o 

■O 

«— * cn  co  cn  in  o 

e 

CD 

o O 

CD 

a) 

--)rsO<3  O O 

03 

o 

o 

E 

LD  o 

to 

E 

r— 1 If)  If) 

(J 

E 

in 

s_ 

LO 

QJ 

E 

o 

C 

m 

(V 

rvi 

cn 

03 

c 

co  jo  co 

e 

e 

tt  in  in  cn 

e 

cn  cd  in 

X 

in  o w k 

X 

o o o o 

03 

03 

CO  N O CO 

X 

CM  CM  •-« 

E 

O CO  CM  O 

E 

ncMCis 

03 

*3-  in 

E 

in 

+j 

1 k 

o 

o 

4-> 

1— 

X 

LU 

T3 

O in  in  in 

u 

-O 

*3-  T-I  o CD 

03 

id  in  o o 

E 

-O 

1 — 1 v—H  i — H r— | 

a) 

....  O o 

E 

CcJCDNlD  O O 

3 

03 

e 

E 

o ■— < <— l O -HIT) 

C 

r— l r~l  OO  CO 

E 

a> 

o 

T3 

& 

i- 

03 

i- 

E 

«=£ 

03 

> 

e 

*3-  *3-  "3-  -53- 

e 

*=}-  «3-  *3- 

E 

<3-  cn  co  in 

X 

o o o o 

03 

n Oin  ^ 

X 

• — i » — 1 CO  i 

X 

HOCNO 

E 

<30CM 

03 

o> 

«— t 

E 

E 

+-> 

•M 

X 

X 

X 

LU 

LU 

LU 

s — ' 

— ' 

* — ' 

E 

T3 

O O LO  O 

s- 

03 

1-lCOtt  rH 

•o 

rHCMtJ-rH 

CD 

O 1^-  r-H 

a> 

0) 

o o 

CD 

e 

E 

•— l <3-  <— < 

CL 

E 

CM  O 

-X 

E 

»— ( 

Du 

o 

CJ 

E 

O 

3 

O 

3 

C 

E 

co  >a-  rj- 

e 

<3  cn  cd  cn 

E 

CO  ^3-  *3-  *53- 

E 

S-  DlL 

e 

S-  CD  S- 

E 

S-  C7>  S- 

o 

CD  E CL) 

o 

0)  E 03 

O 

03  E 0) 

to 

-M  •»-  E < — O CJ5 

to 

+->•<-  E 1 — O O 

to 

+->•<-  E ■ — O O 

03 

E S-  E — 3 3 

03 

C S-  £ r-  3 3 

03 

E L.  E 3 3 

03 

•r-  CL  3 03  WQ 

03 

•f—  Q.  D3  03  (/)□ 

03 

■<-  CL  Z3  03  in  Q 

in 

3 LO  1/0  U-  < U 

LO 

3 CO  W Ll  cC  O 

cn 

3iotnCL  ec  c_) 

• C\J 
CO 

-»->  II 

CO 

“O  c 
O 4— 


CO  -Q 


median  value  = minimum  value. 


Table  57.  Seasonal  summary  of  metals  for  site  DA0208  (right  bank)  on  the  Athabasca 
River  (1976  to  1983).  Winter  (November  1 to  April  15),  spring  (April  16 
to  June  30),  summer  (July  1 to  August  31),  fall  (September  1 to  October  31). 
All  values  in  ug/L. 
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median  value  = minimum  value. 


Table  58.  Seasonal  summary  of  metals  for  site  DD0010  (Centre)  on  the  Athabasca 

River  (1976  to  1983).  Winter  (November  1 to  April  15),  spring  (April  16 
to  June  30),  summer  (July  1 to  August  31),  fall  (September  1 to  October  31). 
All  values  in  ug/L. 
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median  value  = minimum  value. 
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10.3  WATER  QUALITY  DATA  SUMMARY  OF  ROUTINE  PARAMETERS  (TABLE  59), 
NUTRIENTS  (TABLE  60),  AND  METALS  (TABLE  61)  FOR  WATER  BODIES 
IN  THE  NORTHERN  REGION. 


Table  59.  Northern  region  water  quality  data  summary  (routine  parameters),  1976  to  1980. 
See  Figure  4 for  site  location. 
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Table  60. 

Narthern  region  water  quality  data  summary  ( 

nitrogen  and  phosphorus) 

Site 

Total  Phosphorus  (ug/L) 

Or 

■thophosphate  (uq/L) 

Code3 

Description 

n 

med  max 

n 

med  max 

DD0090 

Jackfish  Creek 

12 

40  93 

12 

9 

18 

DD0212 

Big  Point  Channel 

47 

55  930 

47 

14 

235 

KF0100 

Mamawi  Lake  Channel 

13 

51  390 

12 

6 

280 

KF0101 

Chenal  des  Quatre  Fourches 

33 

65  470 

32 

14 

360 

KF0140 

Prairie  River 

12 

70  310 

12 

9 

31 

NA0030 

Riviere  des  Rochers 

48 

48  370 

48 

12 

70 

KF0200 

Lake  Claire  (near  Birch  R.) 

12 

71  370 

12 

13 

30 

KF0201 

Lake  Claire  (west  of  Willow 

Pt.) 

11 

80  230 

11 

10 

180 

DD0080 

Richardson  Lake  (center) 

9 

42  138 

9 

6 

19 

MD2000 

Lake  Athabasca  - Sandy  Pt. 

9 

12  40 

9 

4 

9 

ASWO 

50 

CDWG 

Total  Kjeldahl  N (mq/L) 

Ammonia  (mq/L) 

Nitrate  & 

Nitrite 

Code 

n med  max 

n 

med 

max 

n 

med 

max 

DD0090 

12  0.86  1.90 

11 

0.03 

0.17 

12 

0.03 

0.27 

DD0212 

33  0.60  1.40 

47 

0.04 

1.10 

47 

0.04 

0.66 

KF0100 

13  1.06  3.26 

12 

0.03 

0.10 

12 

0.02 

0.34 

KF0101 

31  0.59  9.15 

32 

0.05 

0.16 

32 

0.06 

0.36 

KF0140 

12  1.62  6.04 

11 

0.06 

0.34 

12 

0.03 

0.27 

NA0030 

47  0.50  6.00 

47 

0.02 

0.11 

49 

0.04 

0.80 

KF0200 

12  1.12  3.33 

11 

0.04 

0.15 

12 

0.03 

0.24 

KF0201 

11  1.16  3.35 

11 

0.05 

0.08 

11 

0.07 

0.20 

DD0080 

9 1.40  1.78 

8 

0.05 

0.31 

9 

L0.01 

0.17 

MD2000 

9 0.37  1.20 

8 

0.03 

0.14 

9 

0.01 

0.08 

ASWO 

CDWG 

9 Rivers,  00AT07-;  Lakes,  01AT07-. 


Table  61.  Northern  region  water  quality  data  summary  (metals).  All  values  in  ug/L. 
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10.4  WATER  QUALITY  DATA  SUMMARY  OF  ROUTINE  PARAMETERS  (TABLE  62), 

NUTRIENTS  ( TABLE  63),  AND  METALS  ( TABLE  64)  FOR  WATER  BODIES 
IN  THE  EASTERN  REGION 


Table  62.  Eastern  region  water  quality  data  summary  (routine  parameters),  1976  to  1982. 
See  Figure  4 for  site  locations. 
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Chloride  Sulphate  Dissolved  Oxygen 

Code  (mq/L)  (mq/L)  (mg/L) 

(00AT07-)  n med  max  n med  max  n med  max 
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Total  Organic  Carbon  Filterable  Residue  Nonfil terable  Residue 

Code  (mq/L)  (mq/L ) (mq/L) 

(00AT07-)  n med  max  n med  max  n med  max 
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Total  Inorganic  Carbon  Phenolic  Material  Oil  and  Grease 

Code  (mq/L)  (mq/L)  (mq/L) 

(00AT07-)  n med  max  n med  max  n med  max 
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Table  63.  Eastern  region  water  quality  data  summary  (nitrogen  and  phosphorus),  1976  to  1982. 
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Total  Phosphorus  Orthophosphate 

Code  (uq/L)  (uq/L) 

(00AT07-)  n med  max  n med  max 
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Table  64.  Eastern  region  water  quality  data  summary  (metals),  1976  to  1982.  All  values  in  ug/L 
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10.5  WATER  QUALITY  DATA  SUMMARY  OF  ROUTINE  PARAMETERS  ( TABLE  65), 
NUTRIENTS  (TABLE  66),  AND  METALS  (TABLE  67)  FOR  WATER  BODIES 
IN  THE  WESTERN  REGION 
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Table  66.  Western  region  water  quality  data  summary  (nitrogen  and  phosphorus). 


Code 

(00AT07-) 

Site 

Description 

Total  Phosphorus  (ug/L) 
n med  max 

Orthophosphate  (uq/L) 
n med  max 

DA0070 

Poplar  Ck.  at  Hwy.  63 

62 

47 

240 

62 

13 

66 

DA0098 

El  1 s R.  near  mouth 

11 

37 

170 

11 

10 

34 

DA0100 

Upper  Ells  R. 

27 

37 

150 

27 

13 

40 

DAO 110 

Unnamed  Ck.,  N of  Ft.  MacKay 

17 

17 

130 

17 

15 

26 

DAO 120 

Asphalt  Ck.  near  Ft.  MacKay 

10 

113 

630 

10 

10 

40 

DA0121 

Eymundson  Ck.  upstr.  Asphalt 

8 

150 

360 

8 

17 

80 

DA0130 

Pierre  R. , N of  Ft.  MacKay 

11 

66 

230 

11 

11 

20 

DAO 140 

Calumet  near  Ft.  MacKay 

13 

124 

720 

13 

37 

72 

DAO 150 

Lower  Tar  R. , NW  of  Ft.  MacKay 

15 

120 

1400 

15 

19 

50 

DA0160 

Joslyn  Ck.  above  Ells 

22 

145 

720 

22 

15 

110 

DA0170 

Lower  Ells  R.  above  Joslyn 

18 

32 

340 

18 

8 

60 

DA0180 

Beaver  R.  above  Syncrude 

18 

120 

490 

18 

20 

235 

DA0181 

Beaver  R.  at  Hwy.  63 

21 

50 

600 

21 

L10a 

46 

DA0182 

Bridge  Ck.  diversion  at  Hwy  63 

32 

57 

220 

32 

18 

80 

DAO 183 

Bridge  Ck.  Beaver  R. 

17 

41 

168 

17 

17 

90 

DAO 190 

Upper  Tar  R. , NW  of  Ft.  MacKay 

12 

70 

243 

12 

14 

50 

D80011 

MacKay  R.  at  Hwy.  63 

37 

65 

260 

37 

24 

ICO 

DB0020 

Dover  R.  above  MacKay  R. 

17 

69 

300 

17 

13 

39 

DB0030 

Dunkirk  R.  near  Ft.  MacKay 

27 

91 

500 

28 

20 

180 

DB0040 

Thickwood  Ck.  above  MacKay  R. 

15 

50 

240 

15 

13 

40 

ASWO 

CDWG 

50 

Code  Total  Inorganic  Carbon  (mg/L)  Phenolic  Material  (mq/L)  Oil  and  Grease  (mg/L) 

(00AT07-)  n med  max  n med  max  n med  max 


DD0070 

29 

29 

93 

59 

0.003 

0.020 

60 

0.8 

4.3 

DA0098 

13 

19 

22 

13 

0.001 

0.015 

12 

0.9 

1.6 

DAO 100 

27 

11 

49 

17 

L0.001 

0.020 

18 

0.4 

1.6 

DA0110 

17 

30 

51 

15 

L0.001 

0.022 

14 

0.5 

2.2 

DAO 120 

10 

13 

61 

10 

L0.001 

0.020 

10 

0.4 

1.7 

DA0121 

11 

16 

150 

11 

L0.001 

0.007 

9 

0.5 

2.0 

DA0130 

11 

15 

28 

10 

L0.001 

0.029 

10 

0.4 

1.0 

DA0140 

13 

43 

112 

11 

L0.001 

0.026 

11 

0.5 

1.4 

DAO 150 

15 

33 

65 

13 

0.001 

0.C24 

13 

0.1 

1.6 

DA0160 

22 

33  , 

81 

20 

L0.001 

0.023 

20 

0.7 

1.5 

DA0170 

18 

19 

25 

14 

L0.001 

0.024 

14 

0.5 

8.2 

DAO 180 

17 

30 

105 

14 

L0.001 

0.023 

14 

0.3 

2.2 

DA0181 

21 

36 

95 

18 

L0.001 

0.010 

19 

0.7 

1.8 

DA0182 

27 

27 

57 

29 

0.001 

0.031 

30 

0.7 

4.9 

DAO 183 

- 

- 

- 

17 

0.006 

0.016 

17 

1.4 

3.1 

DA0190 

12 

24 

70 

11 

L0.001 

0.019 

11 

0.2 

2.2 

DB0011 

34 

26 

90 

32 

L0.001 

0.010 

34 

0.8 

2.5 

DB0020 

17 

56 

92 

13 

L0.001 

0.023 

13 

0.4 

1.7 

DB0030 

29 

25 

62 

20 

L0.001 

0.014 

20 

0.2 

1.2 

DB0040 

15 

15 

36 

9 

L0.001 

0.028 

10 

L0.1 

4.4 

ASWO 

CDWG 


a 


L = less  than 


Table  67.  Western  renion  water  quality  data  summary  (metals).  All  values  in  ug/L. 
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Table  67.  Continued  . 
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10.6  WATER  QUALITY  DATA  SUMMARY  OF  ROUTINE  PARAMETERS  (TABLE  68), 
NUTRIENTS  (TABLE  69),  AND  METALS  (TABLE  70)  FOR  WATER  BODIES 
IN  THE  SOUTHERN  REGION 
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Table  69.  Southern  region  water  quality  data  summary  (nitrogen  and  phosphorus). 
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Table  70.  Southern  region  water  quality  data  summary  (metals).  All  values  in  ug/L. 
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